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Inhaltsübersicht 


Von gewissen Reihenentwicklungen ausgehend werden elementare Appro- 
ximationen in geschlossener Form erhalten für die Zustandsgleichung eines 
idealen Fermi-Dirac-Gases bei beliebigem Grad der Entartung, ohne von 


irgendwelchen numerisch berechneten Werten Gebrauch zu machen. Hierbei mye: 


wird verlangt, daß jede zulässige Approximation in den beiden Extremfällen, 
nämlich der fast vollkommenen Entartung bzw. Nichtentartung, nicht nur 
den Druck des Gases, sondern auch seine spezifische Wärme richtig wieder- 
geben soll. Gewisse Ergebnisse werden tabellarisch mit bekannten Werten 
verglichen. 


In einer früheren Arbeit!) habe ich derThomas- mi- 
Funktion y (x x), (0 < x< oo), in geschlossener Form betrachtet. Die Unter- 
suchung ging von den bekannten Eigenschaften a 


aus. Weiterhin wurden auch auf rein numerischen Wege berechnete Funk- 
tionswerte hinzugezogen. Nun ist es möglich, die damals verwendete Methode 
auch auf andere der in der theoretischen Physik vorkommenden transzen- 
denten Funktionen erfolgreich auszudehnen. In anderen Worten, ist f(x) eine 
solche Funktion, so versucht man einfache algebraische Funktionen g (x) zu 
finden, welche f (x) über den ganzen Definitionsbereich z,< x < 2, genügend 
genau approximieren, wobei von etwa bekannten Reihenentwicklungen von 
f(x) in der Nähe von x, und x, — mögen sie absolut oder nur asymptotisch 
konvergieren — Gebrauch zu machen ist. Unter Umständen kommt man dann 
sogar ganz ohne irgendwelche numerische Behandlungen aus. 

In der vorliegenden Arbeit möchte ich nun gerade ein solches Beispiel 
betrachten, und zwar handelt es sich um die Zustandsgleichung eines idealen 
nicht-relativistischen Fermi-Dirac-Gases für beliebigen Grad der Entartung. 
Es wird angenommen, daß dasselbe aus N gleichen Teilchen besteht, von denen 
jedes die Masse m und den Spin % besitzt. Es sei P der Druck des Gases, 
T seine absolute Temperatur und V sein Volumen, also n = N/V die Teilchen- 
dichte. Dann kann man immer die Zustandsgleichung in der Form 


1) H. A. Buchdahl, Ann. Physik 17, 238 (1956). A) 
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346 Annalen der Physik. 7. Folge. Band 3. 1959 I 
anschreiben (k = Boltzmannsche Konstante), wobei x zu n T” pro- 
portional ist, wie aus dem Virialsatz und den allgemeinen Prinzipien der 
Thermodynamik leicht folgt. Es wird nun eine elementare Approximation 
9 (x), für die Funktion f (x) gesucht. Die Anforderungen, die dabei an g (x) 
zu stellen sind, werden weiter unten im vierten Abschnitt erörtert. Zuerst 
sollen die genauen Gleichungen, welche die Grundlage der gegenwärtigen 
Untersuchung bilden, kurz zusammengestellt werden. 

Es soll noch bemerkt werden, daß die spezifische Wärme bei konstantem 
Volumen (pro Teilchen) allgemein durch e Formel 


3 


2. Die Grundgleichungen 


Man definiere e* durch die Aussage, daß die Anzahl der Quantenzustände 
mit Energie <e* gerade N ist. Dann zeigen z. B. Fowler u. Guggenheim?), 
daß P durch das Integral 


e 
ö 
gegeben ist, wobei sich der Parameter y aus der Beziehung 
dt 
berechnet. Ferner ist die Energie des Gases _ Fah wes yee 
| 
U=ZPV, ‘> (6) 
während &* mit » durch 
verbunden ist, wobei } die Plancksche Konstante ist. Schreibt man E 


(und das sind gerade die von McDougall und EN auf numerischem 
Wege berechneten Funktionen), so kann man auch schreiben 


2) R. H. Fowler u. E. A. Guggenheim, Statistical Thermodynamics, Cambridge 
University Press, 1939, Kap. 11. 

8) J. McDougall u. E. C. Stoner, Philos. Trans. Roy. Soc. London, (A) 287, 67 
(1938). 
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gegeben ist. Schließlich erhält man aus (6), (9), und (10) die der Gl. (3) ent- yes — 


sprechende Formel 
mit F’ — OF /oy. <j 


8 Reihenentwicklungen 


Die folgenden Entwicklungen sind bekannt*). Für y< 0 hat man de 
zu. a, welche für y < 0 absolut konvergieren : et 


G(y) = eu (1 — + + (13) 


Andererseits hat man für genügend große positive y (y > co) die asympto- 
tischen Reihen 


(1 — 2s + 4 0 (68 (12) 


N) 


2 Int 
Fir große T, also y— —oo, ergibt (9) dann natürlich sofort die klassische 
Zustandsgleichung P=nkT. Schreibt man jetzt 
F(y)=x, oder umgekehrt y = q(z), (16) 


so kann man oben überall y durch x ausdrücken, und erhält dann schließlich : 
auf diese Weise die gewiinschte Funktion 


f (2) = @ (q (2) (4 (x) = 216 (q (2). 


2 


Explizit ausgeschrieben hat man also wegen (12) und (14) für genügend kleine ER, Pr 
positive x 
Vx n 3y3 
und für oo 
= (2) (5 =) 1-5 -5(5) +]. (19) 
Endlich ist also 
2 (yas) ~ (sya ~*) (yes) 
2 bn2 (3 2\-‘/s /3 
[1+ (=) + (2+ 00). ey 
Hierbei ist x natürlich durch (10) d.h. 
gi 


mk)! 


Ähnliche Reihenentwicklungen sind schon z. B. von Stoner‘) aufgestellt 
vorden. Er hat aber dann nur U und c direkt ausihnen berechnet. 


4) E.C. Stoner, Philos. Mag. 21, 145 (1936). 
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4. Über Forderungen an die Güte der Approximation 


Be; Die Funktion f (x), deren Reihenentwicklungen in der Nähe der unteren 
J _ und oberen Grenzen des betrachteten Gebietes durch (20) und (21) gegeben 
ee sind, ist diejenige, deren Approximation g (x) nun gesucht wird. Es entsteht 
5 jetzt die Frage, welche Forderungen man jetzt im Prinzip an die Güte einer 
solchen Approximation stellen soll. Betrachtet man nämlich etwa die Funktion 


so ist zwar g(0)=1 und g(z)w +) gi wenn £—00; d.h. (23) verletzt 


nicht die Gesetze der klassischen Thermodynamik, und steht im Einklang mit 
den Extremfällen, in welchen das Gas überhaupt nicht bzw. vollkommen ent- 
artet ist. Es ist aber weder die Zustandsgleichung des schwach entarteten 
Gases noch die spezifische Wärme des fast vollkommen entarteten Gases 
durch (23) richtig wiedergegeben. Es soll also verlangt werden, daß g (x) 
sowohl das schwach entartete Gas als auch das stark entartete Gas samt 
seiner spezifischen Wärme in erster Näherung genau richtig beschreiben soll, 
während es gleichzeitig für beliebigen Grad der Entartung den Gasdruck mit 
nur geringem Fehler darstellen soll. Um diesen Anforderungen gerecht werden 


ay gu können, muß man also schon Funktionen betrachten, die mindestens jeweils 
u die zwei ersten Glieder von (20) und (21) genau wiedergeben. Eine solche 

Funktion g (x) muß also mindestens drei verfügbare Konstanten enthalten, 
= ri wenn schon vorausgesetzt ist, daß g (0) =1 ist. Dabei soll hier versucht 
werden mit den obigen Entwicklungen allein auszukommen, d. h. ohne irgend- 


welche Zahlentabellen zu benutzen. 


+ 


Am naheliegendsten ist wohl eine Funktion der Form 


g(t) aa +b 22-213 Le g2n/3y1n (n > >): (24) 
wobei n einen möglichst ungekünstelten Wert besitzen soll. Hat man n einmal 
gewählt, dann kann man dann die drei Konstanten a, b, c so bestimmen, daß 
jeweils die gewünschten Glieder von (20) und (21) genau wiedergegeben werden. 
Faßt man nun jeweils die dritten Glieder ins Auge, ohne sich dabei um deren 
konstante Faktoren zu kümmern, und verlangt ihre richtige Abhängigkeit 
von x, so müßte einmal » = 5 und das andere Mal n = 11/2 sein. Diese 
Werte sind sich glücklicherweise verhältnismäßig ähnlich. Nimmt man nun 
n = 5, so haben die dritten Glieder der Reihen, die aus (24) entspringen, das 
falsche Vorzeichen, d.h. man könnte erwarten, daß überall g (x) > f (x) ist. 
Andererseits, mit n = 11/2 findet man, daß die in Frage kommenden Glieder 
zwar das richtige Vorzeichen haben, die Absolutbeträge ihrer konstanten 
Koeifizienten aber zu groß sind. In diesem Fall kann man also wohl überall 
f (x) > g (x) erwarten. Ein recht günstiger Wert von n wäre also vielleicht 
der Durchschnittswert von 5 und 11/2, also rn = 21/4; und es ist dann zu er- 
warten, daß g(x) — f(x) teilweise positiv und teilweise negativ ausfallen 
wird. Für n = 21/4 ist die Approximation natürlich etwas „gekünstelter“ 
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5. Die Approximation n = 5 


Für f (x) hat man nun die Approximation 
g(a) jh, (25) 


Werden jetzt die Koeffizienten a, b, c auf die oben auseinandergesetzte Weise 


bestimmt, so erhalt man aa 


= 0,0395618, 
Mit (25) wird (3) zu ic A 
3,21 3 
Tabelle 1 

y x g (x) f(x) 5%) (% 
— ‚016128 | 1,0032 | 1,0032 0 | 1,498 1,497 ee oe 
—3 ‚043366 | 1,0087 | 1,0087 , 0 1,493 1,494 Be 
‚11456 | 1,0230 | 1,0228 | ‚02 | 
‚29050 | 1,0588 | 1,0576 | 11 | 1,456 1,458 | —,2 

0 ‚62809 | 1,1376 | 1,1834 | ,87 | 1,398 10 | —5 

1 1,3964 | 1,2801 | 1,2708 | 73 | 1,310 | 1314 | 

2 2,5025 | 1,4884 | 1,4752 | ,90 1,202 1,193 ee: 

3 3,9770 | 1,2498 | 1,7356 | ,82 | 1,084 1,061 

4 5,7707 2,0496 | 2,0365 | ,66 | 9639 9341 | 320 

5 7,8380 | 2,3762 | 2,3648 | | ‚8522 8210 | 38 

6 10,144 | 2,7214 | 2,7116 | ,36 | ‚725 

7 | 12,665 | 3,0795 | 3,0711 27 6711 ‚644 2 

8 | 15,380 | 3,4466 | 3,4395 | ‚21 ‚6013 576 430 

9 | 18278 | 3,8203 | 3.8143 | 16 | 521 4,2 

10 21,344 | 4,1989 | 4,1938 | 12 ‚4933 ‚474 4,0 

11 24,572 | 4,5812 | 4,5767 | 10 ‚4515 ‚434 4,1 

12 | 27,952 | 4,9663 | 4,9624 | 08 ‚4157 ‚400 4,0 

13 | 31,477 | 5,3536 | 5,3602 | ‚06 | ,3848 ‚371 3,8 

14 35,143 | 5,7426 5,7397 | 05 ‚3580 ‚345 3,8 

15 38,943 | 6,1331 61305 | ,04 ‚3346 ‚323 3,7 

16 42,873 | 6,5248 | 6,5225 | ,04 ‚3139 ‚303 3,6 

17 46,929 | 6,9174 | 6,9154 | ,08 ‚2955 ‚286 3,5 

18 51,106 | 7,3109 7,3090 | ‚2791 271 3 

19 55,402 | 7,7050 | 7,2084 | ‚02 2644 ‚257 3 

20 59,813 | 8,0997 | 8,0983 | ‚02 | 2612 ‚244 3 


In der Tab. 1 sind die aus (25) und (26) berechneten Werte von g und c 
zusammengestellt für den Intervall 0,016 << «< 60, d.h. —4< y< + 20, 
und mit aus den Tafeln von McDougall und Stoner’) berechneten Werten 
von f(x) und c(=c®) verglichen®). Man sieht, daß wie erwartet wurde 
9 — f überall positiv ist. Die prozentuelle Differenz 6 hat einen Höchstwert, 


5) Für y > 5 ist c aus ihnen nur mit dreistelliger Genauigkeit berechenbar. 


q 
u 
= 
| | 


der unter eins liegt. 
und wächst dann bis auf ungefähr + 4 an. 
 dangsam ab. 


2 


Es ist jetzt 


6. Die Approximation n = 21/4 
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63 x? 18 "he 

1250 \125 


— 2125 


Für c ist die prozentuelle Differenz 6’ zuerst negativ, 
Leider fällt danach 6” nur recht 


3 15 13 aie 
Tabelle 2 
y x g(x) | f(z) 6 (%) 5’ (%) 
ool 0,16128 | 1,0032 | 1,0032 ee 1,498 1,497 ine 
onl 0,43366 | 1,0086 | 1,0087 Be 1,494 1,494 peas 
‚11456 | 1,0225 1,0228 | —,03 1,484 1,484 
=n ‚29050 | 1,0564 1,0576 | —,11 1,459 1,458 an 
0 67809 | 1,1802 | 1,1834 | —,28 1,406 1,404 to 
1 1,3964 1,2655 | 1,2708 | —,41 1,310 1,314 | — 3 
2 2,5025 1,4695 | 1,4752 | —,39 1,183 1,198 | — ‚8 
3 3,9770 1,7308 | 1,7356 | —,28 1,047 1,061 
4 5,7707 2,0329 | 2,0365 | —,17 ‚9199 9341 | —15 
5 7,8380 2.3624 | 2,3648 | —,10 ‚8088 8210 | —1,5 
6 10,144 2,7102 2,7116 | —,05 ‚150 | ‚725 —1,4 
7 12,665 3,0704 3,0711 | —,02 ‚369 644. | —1,1 
8 15,380 3,4393 3,4395 de ‚5720 676 | —1,0 
9 18,278 3,8143 3,8143 ate ‚5177 ‚521 — ,06 
10 21,344 4,1939 4,1938 vn ‚4720 ‚474 ne 
11 24,572 4,5770 4,5767 an ‚4332 ‚434 As 
12 27,952 4,9627 | 4,9624 ves ‚3999 400 ‘hin 
13 31,477 5,3505 | 5,3502 ‚3712 ‚371 
14 35,143 | 5,7400 | 5,7397 die ‚3461 ‚345 = 
15 38,943 | 6,1308 | 6,1305 u ‚3240 ‚323 er 
16 42,873 6,5228 | 6,5235 ‚3046 ‚303 
17 46,929 6,9157 6,9154 > ‚2873 ‚286 Sal 
18 51,106 7,3093 7,3090 a ‚2718 ‚271 > 
19 55,042 7,7036 7,2034 ‚2579 ‚257 
20 59,813 8,0985 | 8,0983 en ‚2452 ‚244 si 


In Tab. 2 werden genau wie in Tab. 1 die aus (27) und (28) berechneten 
Werte mit den aus den Tabellen entnommenen verglichen. Differenzen, welche 
durch Unsicherheiten in der letzten Stelle mehr oder weniger ihren Sinn ver- 
lieren, sind durch Striche ersetzt worden. 


Wie erwartet wurde, ist diese Approximation bedeutend besser als die durch 


die Wahl x = 5 bedingte. 


ö wechselt jetzt tatsächlich das Vorzeichen; be- 
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sonders bemerkenswert ist aber die Verbesserung von g(x) fir x = 15, 
— der relative Fehler bleibt unterhalb einem Teil in zehntausend!, — und von 
c für => 20. Überhaupt sind die Ungenauigkeiten außerhalb des Intervalls 
01S 2 20 wohl für alle praktischen Zwecke vollkommen vernachlässigbar, 
wobei auch noch im Auge zu behalten ist, daß sich ja sowohl g als auch c 
asymptotisch genau richtig verhalten. Man sieht hier also, wie man durch 
Zusammenstückeln von Reihenentwicklungen auf obigem Wege zu sehr brauch- 
baren elementaren Approximationen in geschlossener Form gelangen kann. 


Tasmania (Australia), Physics Department, University of Tasmania, 
Hobart. 


Bei der Redaktion eingegangen am 20. November 1958. = 
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der Linien des Röntgenspektrums 
Von T.T. Kakuschadse ERW 


Inhaltsübersicht 


Mit Hilfe der Störungstheorie gelang es F. Bloch!), das Entstehen des 
Kxz-Satelliten der Ka,-Linie zu deuten. Dieser Theorie von Bloch wider- 
spricht jedoch das Fehlen eines analogen Satelliten bei der Linie Kß,, sowie 
das Vorhandensein eines sehr intensiven Satelliten Kp’. 

Eine Berücksichtigung der effektiven Wechselwirkung der Elektronen 
mit den Wärmeschwingungen des Gitters?) erweitert die Anwendbarkeit der 
Blochschen Theorie und erlaubt von einem einheitlichen Standpunkt aus 
nicht nur die Bildung der Satelliten Ka, und Kf’, sondern auch das Erschei- 
nen einer Asymmetrie der Kx,- und Kxa,-Linien zu erklären. 


Zur Untersuchung der gestellten Aufgabe benutzen wir folgende fünf 
Voraussetzungen. 


1. Die relative Wahrscheinlichkeit des Entstehens eines Satelliten durch 


doppelten Übergang ist durch die Formel Bee b 
’ W = W,- W. 1 
MT gegeben; wobei W,- die relative Wahrscheinlichkeit der 
M Anregung des Atoms zur Ausstrahlung des Satelliten 
und W,- die relative Wahrscheinlichkeit des gleichzeitigen 
“tes | doppelten Ubergangs der Elektronen, die an der Bildung 

he T T des Satelliten teilnehmen, ist (Abb. 1). 

_ I | 2. Nach der Theorie des doppelten Überganges ent- 
+ rrr? steht ein kurzwelliger Satellit, wenn gleichzeitig mit dem 
ir Kd, Ka’ Übergang, der der erzeugenden Linie entspricht, ein 


Abb.1.Schematische Übergang eines äußeren Elektrons auf ein tieferes Niveau 
Darstellung der Ent- stattfindet. Nach der Quantenstrahlungstheorie®) erhielt 
— Pn es in diesem Falle Bloch?) fiir die relative Wahrschein- 
Pfeile = rn lichkeit W, des gleichzeitigen doppelten Elektronen- 
bareSatelliten,punk- überganges den angenäherten Ausdruck: 


tierte Pfeile — nicht 5 
beobachtbare Satel- e (2) 
liten W, (2) 


» F. Bloch, Physic. Rev. 48, 187 (1935). 

2, T.I. Kakuschadse, Ubergangsmetalle und Legierungen, 1957; N. § S. 
Akulow u. T. I. Kakuschadse, DAH CCCP 77,593 (1951). 

3) L.L. Schiff, Quantum Mechanic, New York 1955. 
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Hier sind r, und r, die Radien der inneren und äußeren Bahnen im zwei- 
fach ionisierten Atom des Elementes mit der Ordnungszahl Z + 1; E,— E, 
ist die Energie der Verschiebung des Satelliten von seiner erzeugenden 
Linie. 

3. Nach der Theorie des doppelten Uberganges entsteht ein langwelliger 
Satellit, wenn gleichzeitig mit dem Übergang, der der erzeugenden Linie ent- 
spricht, ein Sprung eines äußeren Elektrons auf ein noch höheres Niveau 
stattfindet. In diesem Falle ist die relative Wahrscheinlichkeit W, dafür, 
daß beim Übergang des ersten Elektrons aus dem Zustand mit der Wellen. 
funktion y’(1) in den Zustand w’’(1) (erzeugende Linie), ein zweites Elektron 
aus dem Zustand y’ (2) in den Zustand y” (2) angeregt wird, proportional 
dem Quadrat des Matrixelementes?) 


P= (1)y (1), (2) y (2), (3) 


wobei dr, und dr, die Elementarvolumina im Konfigurationsraum des ersten 
und des zweiten Elektrons darstellen. 1,5 ist der Abstand zwischen den 
Elektronen. 

4. a) Die effektive Wechselwirkung der Elektronen mit den Wärme- 
schwingungen des Gitters?) führt zum Kreislauf der 4s—3d-Elektronen 
in den quasistetigen Bändern der Übergangsmetalle vom Eisentyp. In Abb. 3 
ist dieser Kreislauf durch die Pfeile angezeigt. Die Lebensdauern der Elek- 
tronen im d-Band, 6-Band und A-Band sind mit r,, 75 und r,, bezeichnet. 
Nach 2) wird das 3d-Elektron bei Zimmertemperatur in der Zeit tT, ~ 10-1! sec 
in das 6-Band gehoben (Abb. 3). In der Zeit t3 ~ 10-25 sec fällt das Elektron 
aus dem 6-Band in die Fermi-Ebene des 4s-Bandes unter Bildung des /- 
Bandes (Abb. 3). Aus dem A-Band fällt das Elektron in t ~ 10-12 see zurück 
in die Fermi-Ebene des 3d-Bandes°), d.h. 


~10- sec, ~ 102 sec, ~ 10-" sec. (4) 

Die Dauer K -> L,,, des Uberganges des Atoms eines Elementes der Eisen- 
gruppe®) ist 
Tx Z 1012 sec. (5) 


7 
Aus der relativen Intensität ve ) folgt, daß die Übergangsdauer K > My, 44; 


Ts S ~ 10-" sec (6) 
beträgt. Die Formeln (4), (5) und (6) führen zu der Ungleichheit 


b) Die Zahl der Elektronen sei proportional der Lebensdauer t in den 
entsprechenden Bändern. Wir beschränken uns auf einmalige Anregungen. 
In diesem Falle ist die Zahl der 3d-Elektronen, die in das 4s-Band überführt 


4) D. Coster u. R. Kronig, Physica 2, 13 (1935); G. Wentzel, Z. Physik 48, 524 
(1927); E. Fues, Z. Physik 48, 726 (1927). 

5) Bei den Elementen Cu—Mn finden die zu den in Abb. 3 gezeigten Übergängen ent- 
gegengesetzten Übergängen seltener statt?). 

6) W. Heitler, Quantum Theory of Radiation, 1954. 

7) H.T.Meyer, Wissenschaftliche Veröffentlichungen aus dem Siemens-Konzern 
7, 108 (1929). 
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werden können, gleich der Zahl N der Atome im betrachteten Metallstück 
zu setzen. Folglich beträgt die Zahl der Elektronen bei Zimmertemperatur 


im d-Band N 
4 im A-Band AN ~ 10-1 N 


5. Im Ubergangsmetall der Eisengruppe unterscheiden wir zwei Arten 
von d — s-Übergängen. 

a) Die Übergänge finden zwischen den 3d- und 4s-Termen eines und des. 
selben Ato.nes statt. Dabei beträgt die Energiedifferenz zwischen den 3d. 
und 4s-Termen eines Übergangselementes der Eisengruppe AZ w2Ry. 
Der 4s — 3d-Übergang findet dabei in der 


Zeit 7,, 10-22 sec statt. In Abb. 2a ist 
eine schematische Darstellung der Energie- 
no. --— differenzen zwischen den 4s- und 3d-Termen 
qT der einzelnen Ionen im Gitter gezeigt. 
b) Die Wechselwirkung von Elektronen 
Jd mit Wärmeschwingungen des Gitters (2) be- 


DE: wirkt Übergänge zwischen den 3d- und 4s. 

Termen zweier verschiedener Ionen. Die 
Abb. 2. Schematische Darstellung Termen beträgt dabei AE~1Ry. In 
der s—d-Übergänge im Metall: Abb. 2b ist schematisch die Differenz 
a — in ein und demselben Atom, zwischen den 3d- und 4s-Termen zweier ver- 
b — in verschiedenen Atomen schiedener Ionen im Gitter dargestellt. In 
éh ey Abb. 3 ist eine schematische Darstellung 
von Übergängen in quasistetigen 4s- und 

8d-Bandern für das Metall im ganzen ge- 
Aw zeigt. Der Satz 5a wird für die Erklärung 
der Entstehung des Satelliten Ka,, und der 


Satz 5b wird für die Erklärung des Satel- 
liten Kf’ benutzt. Die angeführten Sätze 
erlauben die Voraussage einer Reihe von 
Satelliten der Ka,- und Kß,-Linien. 
aie. Wenn gleichzeitig mit K — L,,,- oder 
K > (die die Linien Ka, 
ONE) bzw. Kf, erzeugen*) ein Sprung eines 3d- 
dé Elektrons auf den 4s-Term desselben Ions 
Abb. 3. Schematische Darstellung Oder umgekehrt stattfindet (Abb.2a, Satz 5a), 
des Elektronenkreislaufes in den so sind die Satelliten 1, 2, 7, 8 zu erwarten 
quasistetigen 4s- und 3d-Bändern (Abb. 4). Für Elemente der Eisengruppe 
der u 2 — werden die Energien dieser Satelliten gegen- 
™ über den entsprechenden erzeugenden Linien 
in die Größe AE ~ 2 R y verschoben sein, die ungefähr gleich der Differenz 
der 3d- und 4s-Terme in einem und demselben Ion ist. Vier Satelliten 3, 4, 9 
und 10 (Abb. 4) entstehen, wenn gleichzeitig mit den K>L,,,- oder K> My, yy 
Übergängen ein Sprung des Elektrons aus den angefüllten Teilen des 4s- 
oder 3d-Bandes in das 6-Band oder umgekehrt stattfindet (Abb. 3, Satz 5b). 


b Energiedifferenz zwischen den 3d- und 4s.| 
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Für Elemente der Eisengruppe werden Energien dieser Satelliten gegenüber 
den Energien der entsprechenden erzeugenden Linien um die Größe AE ~ 1 Ry 
verschoben sein, die gleicher Ordnung mit der energetischen Breite des quasi- 
stetigen 4s-Bandes ist. 

Vier Satelliten 5, 6, 11, 12 (Abb. 4) entstehen, wenn gleichzeitig mit den 
K > Lyy- oder K > Myı.111-Übergängen ein Sprung eines 4s-Elektrons aus 
dem A-Band auf die Fermi-Ebene des 3d-Bandes oder umgekehrt stattfindet 
(Abb. 3, Satz 4b). Bei Zimmertemperatur werden die Energien dieser Sa- 
telliten gegenüber den Energien der entsprechenden sie erzeugenden Linien 
um die Größe AH ~0,1 Ry verschoben sein, die größenordnungsmäßig 
gleich der Höhe der Aufschichtung 4 ist?). 

Wir haben noch die Wahrscheinlichkeit der Entstehung der Satelliten 
abzuschätzen, da der entsprechende Satellit nicht beobachtet werden kann, 
wenn diese Wahrscheinlichkeit zu klein ist. 

Nach Bloch!) ist die Wahrscheinlichkeit Wı,z dafür, daß im Moment der 
K-Ionisation des Atoms sein 3d-Elektron auf das 4s-Niveau desselben Atoms 


übergeht (Satz 5a), ungefähr gleich eins, d.h. 


Die Wahrscheinlichkeit, W,,,, daß das K-ionisierte Atom ein 3. Elektron 
besitzt, das in das A-Band geworfen ist (Abb. 3), beträgt nach (8) Satz 4b: 
AN 

Wi,4 0,1. (10) 
Diese Bedingungen (9) und (10) ermöglichen die Entstehung der kurzwelligen 
Satelliten 6, 12, 2,8 (Abb. 4). Jedoch befinden sich nach Formel (7), (Satz 4a), 
vor dem K > M,ı.111-Übergang (Kß’-Ausstrahlung) alle Elektronen, die auf 
die 4s-Terme übergehen, praktisch je 


auf Fermi-Ebene des 3d-Bandes Ka, 
(Ta S tg.) Sie können entweder |? 3| ER [+ 2| 
geichzeitig mit dem K — L,,,-Uber- 
gang übergehen, indem sie den Ka,- Kp’ KB, y 
Satelliten bilden, oder sie bilden ein- az 170 4] 
fache s->d-Ubergiinge, bevor der 

K — M,ı.111-Übergang stattfindet. Das Ee 
K-ionisierte Atom mit einem angeregten Abb. 4 


3d-Elektron auf dem 4s-Term desselben 
oder eines anderen Ions ist bereit zum Aussenden von kurzwelligen Satelliten 
der Linie Kx,, und nicht bereit zum Aussenden von kurzwelligen Satelliten 
der Kß,-Linie. Infolgedessen werden die Satelliten 12 und 8 (Abb. 4) nicht 
beobachtet. 

Es bleibt die relative Wahrscheinlichkeit W, des gleichzeitigen doppelten 
Überganges für die Satelliten 6 und 2 (Abb. 4) nach Formel (2) (Satz 2) abzu- 
schätzen. Für beide Satelliten der Elemente der Gruppe Fe (1) gilt 


a 


2 
(Hier ist Ay der Bohrsche Radius.) 
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Die Verschiebungsenergie der Satelliten beträgt 

AE, ~ 2,25 Ry, AE, ~0,1 Ry. (12) 

An den Formeln (1), (2), (9), (10) (11) und (12) erhalten wir für die relative 

Wahrscheinlichkeit W der Satelliten 6 und 2 (Abb. 4) der Elemente der 

1 

330 

Die Wahrscheinlichkeit W, dafür, daß im Moment des Uberganges in den 

Zustand der K-Ionisation das Atom ein in das 6-Band geworfenes 3d-Elektron 


hat (8) (Abb. 3, Satz 4b) ist An 10-4. Ein entsprechender Satellit fehlt 


sowohl für die Ka,-Linie als auch für die Kf,-Linie. Folglich werden die 
Satelliten 4 und 10 (Abb. 4) nicht beobachtet. 

Das Entstehen der Satelliten 5 und 11 (Abb. 4) ist durch das Pauli- 
prinzip verboten: Das A-Band wird als vollkommen besetzt angenommen. 
Auf diese Weise haben von den acht möglichen Satelliten 5, 6, 4, 2, 11, 12, 
10, 8 (Abb. 4) nur zwei genügende Entstehungswahrscheinlichkeiten, nämlich 
2 und 6 (Abb. 5) 

Wir gehen nun zu den langwelligen Satelliten 1, 3, 7, 9 über (Abb. 4). 
Die Größen p,, P3, Pz, Pg sind nach Formel (3) (Satz 3) bestimmt. Der Ver- 
gleich dieser Größen läßt den Schluß zu, daß py bedeutend kleiner als alle 
übrigen drei p,, pz und p, ist. Daraus folgt, daß die Entstehung des Satelliten 
9 (Abb. 4) sehr viel unwahrscheinlicher ist®), als die Entstehung der Satelliten 
1, 3 und 7 (Abb. 4). 

Das Experiment ergibt das Gegenteil: und zwar werden die Satelliten 1, 
3 und 7 (Abb. 4) überhaupt nicht beobachtet, während die Röntgenspektren 
aller Übergangsmetalle der Eisengruppe und ihrer Legierungen äußerst 
helle Satelliten 9 (Abb. 4) aufweisen, die als Satelliten Kf’ bekannt sind. 
Wir nehmen an, daß diese Kf’-Satelliten ihre Entstehung dem aktiven Ein- 
fluß energiereicher Wärmequanten verdanken. 

Bei endlicher Temperatur finden unter den 3d-Elektronen beliebiger 
Atome von Übergangsmetallen oder Legierungen äußerst häufige d > s- 
Übergänge statt?) (Satz 4 und 5b). Es versteht sich, daß diese Übergänge 
gleich häufig stattfinden, unabhängig davon, ob das Atom sich im Zustand 
der K-Ionisation befindet oder nicht. 

Das K-ionisierte Atom besitzt eine hohe Ionisationsenergie. Im Ion, das 
sich im Zustand derK-Ionisation befindet, kann gleichzeitig mit demX>M 43;- 
Übergang ein Sprung eines 3d-Elektrons auf das 4s-Band stattfinden (auf 
den 4s-Term eines anderen Atoms). Dieser Sprung wird jedoch nicht durch 
Wärmequanten verursacht, wie bei nichtionisierten Atomen, sondern durch 
die Energie des K > M,ı.ı1-Überganges, was zum Entstehen des Kß’-Satelliten 
führt. Tatsächlich ist es wahrscheinlicher anzunehmen, daß ein 3d—4s-Über- 
gang eher auf Grund der Energie des K— My, ,11:- Überganges stattfindet, 
als durch die Wechselwirkung mit Warmequanten, da die Wechselwirkung 
zwischen M,, 4;;- und 3d-Elektronen des betrachteten Atoms sehr viel größer 
ist, als die Wechselwirkung zwischen 3d- und 4s-Elektronen verschiedener 
Atome. 


Eisengruppe die Größen ~0,1 bzw. ~ a5 


8) Die Satelliten 1, 3, 7, 9, (Abb. 4) werden wegen ihrer geringen Entstehungswahr- 
scheinlichkeit nicht beobachtet. 
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Energiereiche Wärmequanten bestimmen durch ihren induzierenden Ein- 
fluß die Energie des d— s-Uberganges. Wir haben also nach Formel (2) 
(Satz 2) für die relativen Wahrscheinlichkeiten W, des gleichzeitigen Über- 
ganges von Elektronen im Falle des Kf’-Satelliten (9 in Abb. 4): 


~ = 0,69R y) (12) 
und im Falle des Satelliten 3 (Abb. 4) te ; 
1 a a 
We na, AE ~ 0,69 Ry). (13) 


Auf diese Weise können also infolge der induzierenden Wirkung der Wärme- 
quanten die Satelliten Kf’ und 3 (Abb. 4) beobachtet werden. Die Satelliten 
1 oder 7 (Abb. 4) werden jedoch nicht beobachtet, da die Wärmequanten im, 
Falle der Übergangsmetalle keine Energien der Größenordnung 2 Ry be- 
sitzen, die für den erzwungenen Übergang in diesem Atom notwendig ist?) 
(Satz 5); deshalb fehlen diese auch in Abb. 5. 

Wir bestimmen nun die relative Intensität W der Ausstrahlung von lang- 
welligen Satelliten Kf’ und 3 (Abb. 4). 

Die Wahrscheinlichkeit W, dafür, daß das Atom, das in den Zustand der 
K-Ionisation übergegangen ist, ein angeregtes 3d-Elektron enthält, beträgt 


177 N, 
(14) 
Hier bedeuten N, und N, die Zahl der Atome mit nicht angeregten bzw. 
angeregten 3d-Elektronen. In der Zeit r,,< tg gelingt es im Durchschnitt 


N, Atomen, in den Normalzustand über- 


zugehen (t,, ist die Lebensdauer des Ka eset 
Elektrons im A-Band (Abb. 3). Im : 2 he 
Gleichgewichtszustand werden in dem- IE 2| 

selben Zeitintervall 7,, von N, nor- -2 -7 0 7 2—eky 
malen Atomen bei N, Atomen 3 d-Elek- Kp’ KB, 

tronen angeregt; folglich ist in einem 9] 

bestimmten Atom (von N,„ Atomen) 7 2 Ry 


die Wahrscheinlichkeit der Anregung Ann 6 


eines 3d-Elektrons proportional 2 


und die Wahrscheinlichkeit dafür, daß die 3d-Elektronen des Atoms nicht 
angeregt bleiben, ist proportional 

.Im Zeitintervall t,, gelingt es einer bestimmten Zahl » von Atomen, die 
sich im Zustand der K-Ionisation befinden, einen K — M,, 41,-Ubergang 
auszuführen. Im Ergebnis solcher Übergänge kann z. B. die Kß,-Linie oder 
der Kß’-Satellit entstehen, je nachdem, ob der Übergang K>M 1.111 ohne einen 
d — s-Sprung oder gleichzeitig mit einem d — s-Sprung stattfindet. Die Über- 
gänge K — M,, 4, hängen nicht davon ab, ob das äußere 3d-Elektron im 
Atom erregt ist oder nicht. Es werden also bei vÜbergängen der FormK >My, 44 


N„—N 
* ° Quanten der erzeugenden Linie und W, ıT- 
n 


Quanten des Satelliten ausgestrahlt, wobei W, die Wahrscheinlichkeit des 


(von N,, Atomen) » 
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Doppelüberganges ist; (z.B. K > M,,ın und 3d — 48). Infolgedessen er- 

halten wir für die relative Intensität des Satelliten zur erzeugenden Linie 
(von N, Atomen) 

N, 

W W, (15) 

und für die relative Intensität des langwelligen Satelliten zur erzeugenden 

Linie (unter Berücksichtigung der Gesamtzahl von (N,,-+ N,) Atomen) 
nach (14) und (15) 

N,N, 


W=W, (16) 
Für die Satelliten Kf’ und 3 (Abb. 4) erhält man nach Formel (16): 
Wee = N2 — N}? 


W sats = (17) 


32 N2— N?° 


_ Der Satellit 3 (Abb. 4), der an der Basis der hellen Ka,-Linie liegt, kann 


unter Umständen infolge starker Diffusität nicht beobachtet werden. 


Folgerungen und experimentelle Angaben 


Die relative Intensität der Ausstrahlung des Kß’-Satelliten des Eisens 
hängt von den Untersuchungsbedingungen ab. Wenn man annimmt, daß die 
Temperatur der Probe bei der Sekundärerregungsmethode ~ 400° K und 
bei primärer Erregung ~ 700° K beträgt, dann ergibt sich nach (17) die rela- 
tive Intensität von Kg, in völliger Übereinstimmung mit dem Experiment 
zu ~ 20% bzw. ~ 50%"). (Auf Grund des experimentell ermittelten 
Temperaturganges der Elektronenwärmekapazitäten bei tiefen Temperaturen, 
sowie aus Messungen der magnetischen Suszeptibilität entnehmen wir, daß 

an ” die Zahl der Eisenatome mit 

Ka, angeregten 3d-Elektronen in 
5 erster Näherung proportional 

der absoluten Temperatur ist 
(N, und sie bei T~ 


S 1200°K Werte von N, 
: —> ky erreicht?)). 

2 . 1%, Die Linie Kx, und der 
Ss Satellit o (Abb. 6) sind nicht 


trennbar. Das Zusammen- 


Abb. 6 mit der ihn erzeugenden Linie 


moon Ka, (punktierte Linien) er- 
gibt eine stark asymmetrische Ka,-Linie (Abb. 6). Vollkommen analog 
entsteht die Asymmetrie der Ka,-Linie in völliger Übereinstimmung mit 


dem Versuch. 


9) Ja. M. Fogel, J. exp. theoret, Physik 10, 455 (1940). 3 kei 
10) G. Ortner, Nature 117, 828 (1926). 
u) M.T. Druywesteyn, Z. Physik 48, 707 (1928). 
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Von den zwölf in Abb. 4 gezeigten Satelliten werden in voller Uberein- % 
timmung mit dem Experiment nur drei Satelliten beobachtet: Satellit 2 
(Abb. 5 oder Abb. 6), der als Ka,-Satellit bekannt ist; Satellit 5 (Abb. 5 oder 
Abb. 6), der als Kf’ bekannt ist und Satellit 6 (Abb. 5 oder Abb. 6), der die 
Asymmetrie der Ka,-Linie bewirkt. ae 
Die beobachtbare starke Diffusität des KA’-Satelliten, die Abhängigkeit 
der Wellenlänge vom Typ der chemischen Bindung und von der Gitterstruktur, 
wwie die ungleiche Verschiebung der Ka, und Kß’-Satelliten von der sie er- R a 
magenden Linien erhärtet vollständig die Rolle der Übergänge zwischen den 2 : 
Bändern?) bei der Satellitenbildung. 2 


4- 


Bei der Redaktion eingegangen am 9. Dezember 1958. 
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a Frage der Existenz von singularitatsfreien Lésungen 
der allgemein-relativistischen Feldgleichungen, 
ate die Teilchenmodelle darstellen könnten 


Von A. Papapetrou und H. Treder al, 


Mit 4 Abbildungen 


Inhaltsübersicht 


7 Die von 1 Einstein gestellte Frage nach der Existenz von singularitats- 
_ freien Lösungen der Gravitationsgleichungen, die die Grenzbedingung g,, Nps 
für r-— oo erfüllen (und die als Modelle von beständigen Teilchen 
verwendbar sein könnten), wurde für den Fall der stationären Lösungen durch 
die Untersuchungen von Einstein und Pauli und von Lichnerowiez 
negativ beantwortet. In der vorliegenden Arbeit wird die Existenz von zeit- 
lich periodischen singularitätsfreien Lösungen der Gravitationsgleichungen 
diskutiert. Aus den Ergebnissen von Papapetrou über das asymptotische 
Verhalten von periodischen Feldern folgt, daß eine solche Lösung nur in einem 
endlichen Innenbereich des dreidimensionalen Raumes, in dem das Gravi- 
_ tationsfeld stark ist, zeitlich periodisch sein könnte und mittels einer Sprung- 

fläche an ein stationäres Außenfeld anzuschließen wäre. Aus Stellmachers 
Untersuchungen über die mit den Feldgleichungen verträglichen Sprünge 
folgt nun, daß diese Sprungfläche eine glatte Nullhyperfläche sein muß. Es 
wird dann ein geometrischer Satz über glatte Hyperflächen hergeleitet, aus 
dem sich ergibt, daß eine derartige Sprungfläche nicht existieren kann. Dem- 
nach erlaubt Einsteins Theorie des reinen Gravitationsfeldes auch keine 
zeitlich periodischen singularitätsfreien Lösungen und somit keine Modelle 
_ bestiindiger Teilchen. — Auf dasselbe negative Ergebnis führt auch die Dis- 
_ kussion der notwendigen Eigenschaften von zeitlich periodischen singularitäts- 
freien Lösungen der Einstein-Maxwellschen Feldgleichungen für das 
kombinierte Gravitations- und elektromagnetische Feld. 


. Einleitung 


Einsteins Programm einer einheitlichen Feldtheorie bedingt die Vor- 
stellung, daß die Elementarteilchen als singularitätsfreie Lösungen der Feld- 
gleichungen darstellbar sind. Nach Einsteins Idee sind die Elementar- 
teilchen als Feldkonzentrationen aufzufassen, deren Existenz durch die 
_ Wirkung des Feldes auf sich selbst gewährleistet wird und somit eine Folge 
der Nichtlinearität der Feldgleichungen ist. 

Einstein!) hat sehr früh die Frage aufgeworfen, ob vielleicht schon die 
allgemein-relativistische Theorie des reinen Gravitationsfeldes derartige 
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Modelle von beständigen Teilchen zuläßt. Derartige Teilchenmodelle würden 
einem Gravitationsfeld g,, entsprechen, das die Einsteinschen Gravita- 

tionsgleichungen 

= 0 (1) 
im ganzen Raum erfüllt. Damit ein solches Gravitationsfeld ein beständiges 
Teilchen darstellt, muß bei passender Wahl des Koordinatensystems ferner 
einer der beiden folgenden Fälle vorliegen: Be 
a) Ju» ist im ganzen Raum zeitunabhängig, d.h. ee 


= Gur (a, 22, 23) (2a) 


oder 

b) gu» hängt mit einer Periode 7’ periodisch von der Zeit ab, d.h. es ist 
Jur Jur 2, 2, x°) Jur (2, x, + T). (2b) 
Im Fall b) werden wir voraussetzen, daß g,, wesentlich periodisch von der 
ritäts- | Zeit abhängt, schließen also jedes g,, aus, das nur scheinbar von der Form 
—>Nu» | (2b) ist und durch eine passende Koordinatentransformation auf die stationäre 

ilchen | Form (2a) gebracht werden kann?). 
durch Die Existenz und die Art solcher singularitätsfreien Lösungen hängt 
»wiez | natürlich von den Randbedingungen ab, die für die g,, gefordert werden. 
1 zeit- | Wir werden unseren Betrachtungen die von Einstein?) schon bei seiner 
ungen | ersten Untersuchung des zeitunabhängigen singularitätsfreien Gravitations- 


tische | feldes benutzte Grenzbedingung 7 


aravi- 
rung- | zugrunde legen, also fordern, daß für r > oo die Metrik asymptotisch in die 
:hers | pseudo-euklidische Minkowski-Metrik übergeht. Aus (3) folgt, daß für 
rünge | genügend große r das Gravitationsfeld noiwendig schwach ist. Die Bedingung 
3. Es | (3) ist im wesentlichen der Annahme äquivalent, daß es physikalisch sinnvoll 
t, aus ist, ein isoliertes Teilchen zu betrachten und anzunehmen, daß die Wechsel- 


Dem- wirkung von zwei sehr weit voneinander entfernten Teilchen sehr schwach 
keine | ist®). 

odelle Für das erste Problem der zeitunabhängigen singularitätsfreien Lösungen 
> Dis- | haben nun zunächst Einstein sowie Einstein und Pauli®) gezeigt, daß aus 


itäts- (1) und (3) das Verschwinden der Gesamtenergie J, des Feldes folgt. Denn 
- das für derartige Felder gilt die Beziehung 


x’ = const x’ = const et 
2) Vgl. A. Papapetrou, Ann. Physik (6) 20, 399 (1957). ni A 
Vor- 8) A. Einstein, Revista (Univ. Nac. de Tucuman) A 2, 11 (1941). 
Feld- 4) A priori darf man wohl nicht den Fall ausschließen, daß die Grenzbedingung (3) 


t physikalisch unzutreffend ist und man vielmehr von der Hypothese eines räumlich end- 
— lichen Weltalls etwa im Sinne der Betrachtungen Eddingtons (vgl. Fundamental 
1 die Theory, Cambridge 1948) ausgehen muß. In diesem Fall wäre die Grenzbedingung (3) 
Folge durch einfache Eindeutigkeitsbedingungen zu ersetzen. Es sei jedoch bemerkt, daß bei 
Zugrundelegung eines räumlich endlichen Weltalls die Behandlung des Einsteinschen 
Teilchenproblems sehr viel komplizierter wäre als bei Verwendbarkeitder Grenzbedin 
: (3). Denn in einem solchen kosmologischen Raum müßten im Prinzip alle Teilchen des 
ırtige Universums gleichzeitig berücksichtigt werden. 

5) A. Einstein, l.c. Anm. 3; A. Einstein u. W. Pauli, Annals of Math. 44, 131 
(1943). 
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und die rechte Seite von (4) verschwindet wegen (1). Liehnerowicz®) hat 

_ dann weiter bewiesen, daß es bei Erfüllung der Grenzbedingung (3) überhaupt 
keine nichttrivialen zeitunabhängigen singularitätsfreien Lösungen der 
Einsteinschen Gravitationsgleichungen (1) gibt: Die vierdimensionale 
Riemannsche Raum-Zeit-Mannigfaltigkeit ist bei Gültigkeit von (1), (2a) und 
(3) mit der ebenen Minkowski-Welt identisch. 


Nach diesem negativen Ergebnis ist nur noch die zweite Möglichkeit, 
nämlich die Frage der Existenz von zeitlich periodischen singularitätsfreien 
Lösungen, zu diskutieren. Nun hat Papapetrou’) für solche Lösungen be- 

= daß sie in dem gemäß (3) existierenden räumlich unendlich ausge- 
_ dehnten Bereich r > a, in dem das Feld schwach ist, notwendig zeitunab- 
hängig sein müssen. Ist also das Gravitationsfeld im 
ganzen Raum schwach, so gibt es überhaupt keine 
zeitlich periodischen singularitätsfreien g,,. Hieraus 
folgt als eine notwendige Bedingung für die Existenz 
eines zeitlich periodischen singularitätsfreien Gravi- 
tationsfeldes der durch Abb. 1 veranschaulichte Feld- 
verlauf: In dem räumlich unendlichen Bereich r > a 
ist das Gravitationsfeld überall schwach und stationär. 
In dem räumlich endlichen Bereich r< a ist das 
Gravitationsfeld hingegen teilweise stark und zeitlich 
periodisch. Dieser Innenbereich muß also in zwei 
durch eine Hyperfliche & voneinander getrennte 
Teilbereiche zerfallen. Innerhalb des von 2 um- 
Abb. 1 schlossenen Gebiets B ist das Gravitationsfeld stark 
und zeitlich periodisch. Im AuBenbereich A, zu dem 
Bars der Bereich r>.a des schwachen Feldes gehört, ist das Gravi- 
tationsfeld hingegen stationär. 


\ Ein derartiger Verlauf des Gravitationsfeldes ist nur dann möglich, wenn 

auf der Hyperfläche & gewisse Ableitungen der g,, unstetig sind. Wären 

' nämlich auf X die Ableitungen jeder Ordnung stetig, so wäre die Fortsetzung 

des zeitunabhängigen Feldes aus dem Bereich A in den Bereich B nach den 

En von Lichnerowicz®) über das Cauchysche Problem der all- 

gemeinen Relativitätstheorie ebenfalls zeitunabhänig. In der Tat wäre dann 

die Funktion g,, im Bereich B einfach die analytische Fortsetzung der Funk- 

tion g„, im Bereich A. Die Hyperfläche Y muß demnach eine Sprungfläche 
des betrachteten Problems sein. 


In dem folgenden $ 2 werden wir die Eigenschaften der Sprungfläche £ 
näher untersuchen. Im $ 3 wird ein geometrischer Satz über glatte Hyper- 
flächen bewiesen. Im $ 4 werden wir dann mit Hilfe dieses Satzes beweisen 
können, daß es in der allgemeinen Relativitätstheorie keine zeitlich periodi- 

schen, die Grenzbedingung (3) befriedigenden, singularitätsfreien Gravitations- 
5 felder gibt. Und wir werden ferner zeigen, daß derselbe Satz auch für die 


= 


u 6) A. Lichnerowicz, C. R. hebd. Séances Acad. Sci. 222, 432 (1946); Theories 
ER de lag gravitation et de l’Electromagnétisme, Paris 1955, S. 198146. 

I 4. 7) A. Papapetrou, Ann. Physik (6) 20, 399 (1957); (7) 1, 186 (1958); 2, 87 (1958). 
! 8) A. Lichnerowicz, ]l. c. Anm. 6, S. 27—33. 
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Einstein-Maxwellsche Theorie des kombinierten Gravitations- und elek- 


tromagnetischen Feldes gilt. 


2. Die Sprungfliche 


2 Ex, 

Der von uns zu betrachtende Raum ist eine reguläre Riemannsche 
Raum-Zeit-Mannigfaltigkeit V,, deren Metrik den Einsteinschen Gravita- 
tionsgleichungen (1) und der Grenzbedingung (3) genügt. Die Signatur der 
Metrik muß überall dieselbe wie die der Minkowskischen Metrik n,„, sein, 
die g,, müssen also lokal immer auf die Minkowskische Form zurückführbar 
sein. Entsprechend der Grenzbedingung (3) ist der V, in den drei unab- 
hängigen raumartigen Richtungen offen. Ferner muß der V, aus physikali- 
schen Gründen auch in Richtung der Zeitkoordinate offen sein. Der betrach- 
tete V, ist somit topologisch homolog zur pseudoeuklidischen Minkowski- 
Welt. 


Für die Beweisführung im $ 3 genügt es zu bemerken, daß es in einem 
solchen V, eine Schar von sich nicht schneidenden raumartigen Hyperflä- 
chen gibt derart, daß jeder Punkt des V, auf genau einer dieser Flächen liegt. 
Durch diese Schar von raumartigen Hyperflächen wird eine dazu orthogo- 
nale Schar von zeitartigen Weltlinien 
definiert®). Dabei kann man das Koor- xe 
dinatensystem so wählen, daß die raum- 
artigen Hyperflächen der betrachteten 
Schar mit den Hyperflächen x = const 
zusammenfalien. Dann sind die zu 
dieser Schar orthogonalen zeitartigen 
Weltlinien die Parameterkurven von x° 
(vgl. Abb. 2). 

Bei dem von uns behandelten Pro- 
blem des zeitlich periodischen singulari- 
tätsfreien Feldes ist es zweckmäßig, die Abb. 2 
Schar der raumartigen Hyperflächen so 
zu wählen, daß sie im Bereich r > a des schwachen Feldes (Abb. 1) mit den 
Hyperflächen x° = const desjenigen Koordinatensystems zusammenfallen, 
in dem dieses Feld zeitunabhängig ist. Es sei noch die für das Folgende 
wesentliche Tatsache hervorgehoben, daß die durch die Gleichung r= a 
bestimmte Hyperfläche glatt und zeitartig ist. 


Wir betrachten nun die auf Abb. 1 gezeigte Sprungfläche 2. Es muß zu- 
nächst die für die folgende Beweisführung wichtige Tatsache hervorgehoben 
werden, daß die F.äche Y ganz im Innern des durch die zeitartige Fläche r = a 
begrenzten Bereichs liegt, diese Fläche also nirgends schneidet oder berührt. 
Wegen der Periodizität des gesamten Gravitationsfeldes wird auch die Sprung- 
fläche Y mit derselben Periode 7 wie g,, zeitlich periodisch sein müssen. 


Da ferner das Feld im Bereich B als wesentlich und nicht nur scheinbar 
periodisch vorausgesetzt wurde, im Bereich A dagegen zeitunabhängig ist, 
muß es auf & nichtforttransformierbare Sprünge gewisser Ableitungen der 
Jur geben. Nach den Ergebnissen von Stellmachers Diskussion des Sprung- 


®) L. P. Eisenhart, Riemannian geometry, Princeton 1949, S. 41—43. 
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problems in der allgemeinen Relativitatstheorie!) muß daher die Fläche F 
eine Nullfläche des betrachteten regulären Riemannschen V, sein. Ist 
also die Sprungfläche XY etwa durch die Gleichung 


2=2 2, 2, 2°) = (5) 


ia gegeben, so ist ihr Normalvektor 


g"” Pu Pv = PP” Py = 0. (6) 


Aus den Feldgleichungen (1) folgt noch eine weitere wichtige Eigenschaft 
‘g der Sprungfläche: Die StoBwellen des Gravitationsfeldes breiten sich längs 
der Bicharakteristiken der Einsteinschen Feldgleichungen aus, d.h. längs 
der Trajektorien des Normalvektors p, der Sprungfläche. Diese Trajektorien 
_ sind die auf der Nullfläche liegenden geodätischen Nullinien mit der Gleichung 


pe = 0. (7) 


_ Ist nun y,, der Sprungtensor der niedrigsten unstetigen Normalableitung 
von gy»), so folgt aus den Feldgleichungen (1) für die ,,Sprungerregung“ 


nach Stellmacher der skalare Erhaltung) 

V-9Ep).=0. (9) 


Aus (9) folgt insbesondere, daß, wenn die Sprungerregung E an irgendeinem 
Punkt einer auf XY liegenden geodätischen Nullinie nicht verschwindet, E dann 
längs der ganzen geodätischen Nullinie von Null verschieden sein wird. Der 
Skalar E ist aber positiv definit in dem Sinne, daß aus E=0 das Ver- 
schwinden des nichtforttransformierbaren Teils von y,, folgt. Daher ergibt 
sich aus (9), daß der nichtforttransformierbare Teil von y,, auf der ganzen 
Nullgeodäten von Null verschieden ist, wenn er an einem ihrer Punkte nicht 
verschwindet. 

Diese aus den Feldgleichungen folgende Eigenschaft der Sprünge hat zur 
Folge, daß die Sprungfläche X notwendig eine einheitliche Nullfläche sein 
muß, welche jede auf ihr liegende Nullgeodäte als Ganzes enthält. Eine solche 
Nullfläche ist notwendig frei von Kanten, kann also nicht aus Teilen von sich 
schneidenden Nullflächen zusammengesetzt sein, da dann die Trajektorien 
des Normalvektors der Sprungfläche & entweder Knicke oder Endpunkte be- 


10) K. Stellmacher, Math. Ann. 115, 740 (1938). In dieser Arbeit ist der Fall be- 
handelt, daß die niedrigsten unstetigen Ableitungen von gy» die zweiten oder höheren 


r Ableitungen sind. Die Stellmacherschen Ergebnisse gelten auch für die Sprünge der 


ersten Ableitungen wie von A. Papapetrou u. H. Treder, Mathematische Nachrichten, 


0 (im Druck) bewiesen wurde; vgl. auch St. O’Brien u. J. L. Synge, Comm. of the 


Dublin Inst. for Advanced Studies A No. 9 (1952). Die gu» selbst können bei Forderung 
einer r n Metrik nicht springen, wie in einer späteren Arbeit gezeigt werden soll. 
11) Vgl. K. Stellmacher, I. c. Anm. 10. 

12) Nach A. Papapetrou u. H.Treder, 1.c. Anm. 10, gilt diese Formel auch für 


die Sprünge erster Ordnung. 
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sitzen würden ®). — Eine einheitliche Nullfläche kann man sich am einfachsten 
durch eine Schar von Nullgeodäten gebildet denken, die eine glatte zwei- 
dimensionale raumartige Fläche V, orthogonal schneiden. Dabei müssen 
nur die Richtungen dieser Nullgeodäten auf der ganzen Fläche V, einheitlich 
gewählt werden, so daß die benachbarten Punkten des V, entsprechenden 
Nullgeodäten kontinuierlich ineinander übergehen“), da sonst keine kon- 
tinuierliche, den vierdimensionalen Raum JV, teilende, Nullfläche 2 ent- 
stehen würde. 

Eine einheitliche Nullfläche hat im allgemeinen an jedem Punkt genau 
einen Normalvektor, ist also eine glatte Hyperfläche. Die einzige Ausnahme 
bilden die sogenannten konischen Punkte: An solchen Punkten treffen sich 
verschiedene auf der Nullfläche liegenden Nullgeodäten mit dem Ergebnis, 
daß die Nullfläche dort keine Normale besitzt. Es ist aber leicht einzusehen, 
daß die für unser Problem brauchbaren Sprungflächen keine konischen Punkte 
besitzen können, sondern überall glatt sein müssen. Dies ergibt sich unmittel- 
bar aus Stellmachers Behandlung des Sprungproblems, in der die Ein- 
deutigkeit des Normalvektors p, eine wesentliche Voraussetzung ist. Formal 
kann man die Unvereinbarkeit von konischen Punkten der Sprungfläche mit 
dem Sprungproblem auf Grund des Stellmacherschen Erhaltungssatzes (9) 
erkennen. Hat nämlich die Sprungfläche XY etwa den konischen Punkt Px, 
so folgt aus (9), daß die Erregung E und somit der Sprungtensor y,, an Px 
unendlich groß werden. Ein regulärer metrischer Tensor g,, läßt sich aber in 
der Umgebung der Sprungfläche in der Form 


Juv = + Yuvh(z) (10) 


entwickeln, wobei g~,, und y,, regulär sein müssen, während h (z) eine geeig- 
nete Sprungfunktion ist"). Wird y,, unendlich groß, so ist diese Entwick- 
lung nicht möglich, d. h. die Metrik g,, wird am Punkt P, singular. 

Für das im folgenden $ 3 abzuleitende geometrische Theorem über das 
Verhalten von einheitlichen Nullhyperflächen ist es jedoch unwesentlich, ob 
die Nullfläche konische Punkte besitzt oder nicht. Daher werden wir im 
§ 3 beliebige einheitliche Nullhyperflächen, also auch solche mit konischen 
Punkten, in die Untersuchung einbeziehen. 


3. Ein geometrischer Satz über glatte Hyperflächen 


Wir legen unseren Betrachtungen wieder einen offenen regulären Rie- 
mannschen V, mit singularitätsfreier Metrik g,, zugrunde, der nach den 
Ausführungen am Anfang des $ 2 topologisch homolog zur pseudo-euklidischen 
Minkowski-Welt ist. Für diesen V, beweisen wir zunächst folgendes Lemma: 


13) Durch das obige Argument wird nicht ausgeschlossen, daß unser Sprungproblem 
mehrere einheitliche Nullflächen Z’ als Sprungflächen besitzt, die sich eventuell schneiden. 
Für die Beweisführung des $ 4 genügt es aber zu bemerken, daß jede dieser Flächen 2 die 
zeitartige Hyperfläche r = a nicht schneiden darf. 

14) Man beachte nämlich, daß es an jedem Punkt des V, zwei zum V, orthogonale 
Nullrichtungen gibt. Die Fläche V, stellt eine momentane Front der sich auf 2 aus- 
breitenden Gravitationsstoßwelle dar. Die obige Bedingung für die Richtung der Null- 
geodäten bedeutet anschaulich, daß bei einem geschlossenen V, entweder eine auslaufende 
oder eine einlaufende Gravitationswelle vorliegt. 

15) H. Treder, Ann. Physik (7) 2, 225 (1958). 
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tae Bar geschlossene glatte Hyperfläche besitzt Punkte, an denen ihr Normal. 
DR ox vektor raumartig, Punkte, an denen er zeitartig, und auch Punkte, an denen 
me er isotrop ist. 


et Der Beweis ergibt sich unmittelbar aus der Abb. 2, auf der neben der Schar 
der zeitartigep Parameterkurven von x und der dazu orthogonalen raum. 
> Ys artigen Hyperflächen x’ = const auch die glatte geschlossene Hyperfläche 8 
i im eingezeichnet ist. Wegen der Geschlossenheit und der Glätte von S wird es 
: ade Er sicher Parameterkurven von x? geben, die die Fläche S an Punkten P, tan- 
m. gieren; an solchen Punkten ist der Normalvektor von S raumartig. Ferner 

oy wird es auch Hyperflächen x = const geben, die die Fläche S an Punkten P, 
_ tangieren, und an solchen Punkten ist der Normalvektor von S zeitartig, 
Wegen der Kontinuität des Normalvektors einer glatten Fläche folgt aber aus 
_ Stetigkeitsgriinden, daß zwischen den Punkten P, und P, auch Punkte liegen 
müssen, an denen der Normalvektor ein Nullvektor ist. 


Aus diesem Lemma folgt sofort, daß in der zugrunde gelegten Riemann- 
schen Raum-Zeit-Mannigfaltigkeit eine glatte Hyperfläche, deren Normal- 
_ vektor überall raumartig oder überall zeitartig oder überall isotrop ist, not- 
wendig eine (in mindestens einer Richtung) offene Hyperfläche ist!%). Das- 
selbe Ergebnis gilt ersichtlich auch für ganze einheitliche Nullflächen mit 
konischen Punkten. 


a ee Wir betrachten im folgenden glatte Hyperflächen, deren Normalvektor 

überall denselben Charakter hat (d.h. an allen Punkten der Hyperfläche ent- 

_ weder raumartig oder zeitartig oder isotrop ist) und werden folgenden Satz 
beweisen: Zwei derartige Hyperflächen mit Normal- 

5,559 sm s' s verschiedenen Charakters müssen sich not- 
wendig schneiden. 

Um konkret argumentieren zu können, wollen 
wir den Fall einer glatten, überall zeitartigen Hyper- 
fläche S und einer glatten Nullfläche XY betrachten 
(Abb. 3). Der folgende Beweis läßt sich dann ohne 
weiteres auf die andern möglichen Kombinationen 
von Hyperflächen übertragen. Wir werden unseren 

p Satz dadurch beweisen, daß wir annehmen, daß 

die Flächen S und 2 der Abb. 3 überhaupt keinen 

\ gemeinsamen Punkt besitzen und dann zeigen, daß 

diese Annahme zu einem Widerspruch führt. 

ais Abb. 3 Von der Fläche S ausgehend, errichten wir ein 

a auf S basierendes GauBsches Koordinatensystem 

in der Weise, daß in ihm die Hyperfläche S die Gleichung x! = 0 hat, 

# uss während die Parameterkurven von x! raumartige geodätische Linien sind, 
die die zeitartige Hyperfläche S orthogonal Die 


xz! = const = b (11) 


; 16) Man überlegt sich leicht an Hand des Minkowskischen Raums, daß jede überall 

Ws, bit raumartige Hyperfläche in allen drei unabhängigen Richtungen offen sein muß, während 

eine zeitartige Hyperfläche in den beiden auf ihr gelegenen unabhängigen raumartigen 
Richtungen geschlossen sein kann. De 


17) L. P. Eisenhart, 1. c. Anm. 9, S. 57—59. eet 
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bestimmt dann eine Hyperfläche S’ (Abb. 3), die ebenfalls zeitartig und 
(mindestens für kleine Werte von b) glatt ist und bei passender Wahl des 
Vorzeichens von 6 näher an der Fläche 2 liegt als die ursprüngliche Fläche 8. 
Wenn das benutzte Gaußsche Koordinatensystem bis zu genügend großen 
Werten der Konstante b verwendet werden kann, werden wir mit seiner Hilfe 
eine Schar von glatten zeitartigen, zu S geodätisch parallelen Flächen kon- 
struieren können, die auch eine Fläche S, enthält, welche die Fläche S schneidet 
(Abb. 3). 

Das Gaußsche Koordinatensystem bricht aber im allgemeinen dort zu- 
sammen, wo sich die zur Ausgangsfläche S senkrechten geodätischen Linien, 
d.h. die Parameterkurven von x!, schneiden. Die Gaußsche Konstruktion 
ührt dann zu einer caustischen Fläche Sx, die die Grenze der Verwendbarkeit 
des Gaußschen Koordinatensystems darstellt. Es kann nun vorkommen, daß 
die Fläche Sx so nahe an der Ausgangsfläche S liegt, daß wir unter Verwendung 
des auf S basierenden Gaußschen Koordinatensystems keine Fläche von 
der Art S, konstruieren können. Liegt nun z. B. die Fläche 8’ so nahe an der 
caustischen Fläche Sx, daß sie zwar noch glatt ist, aber schon Stellen starker 
Krümmung besitzt, so können wir sie zunächst an diesen Stellen glätten. 
Diese Glättung erfolgt dadurch, daß wir die stark gekrümmten Bereiche der 
Hyperfläche 8’’ kontinuierlich durch andere ersetzen, die ebenfalls zeitartig 
sind, aber eine geringere Krümmung aufweisen. Auf der nach dieser Modifi- 
zierung entstandenen glatten zeitartigen Fläche SC) basieren wir dann ein 
neues Gaußsches Koordinatensystem, in dem jetzt $( die Gleichung x! = 0 
bekommt und können so mit der Konstruktion der Schar von glatten zeit- 
artigen Hyperflächen weiter in Richtung von 2 fortschreiten. 

Wenn wir auch mit Hilfe dieses zweiten Gaußschen Koordinatensystems 
nicht bis zu einer Fläche der Art S, gelangen können, werden wir den ge- 
schilderten Prozeß erneut wiederholen. Da die geodätische Entfernung 
zwischen der ursprünglichen zeitartigen Fläche S und der Nullfläche & endlich 
ist, werden wir nach einer endlichen Wiederholung dieses Prozesses zu einer 
geglätteten zeitartigen Hyperfläche $® gelangen, die so nahe an 2 liegt, 
daß das auf ihr basierende Gaußsche Koordinatensystem zwischen S™ und 
2 keine caustische Fläche besitzt, so daß wir von S™ ausgehend mit Hilfe 
des Gaußschen Verfahrens zu einer Fläche der Art S, gelangen können. 


Die Fläche S, schneidet die Fläche &, während 8 mit & voraussetzungs- 
gemäß keinen gemeinsamen Punkt besitzen soll. Demnach muß die Schar 
der zu $(®) geodätisch parallelen glatten zeitartigen Flächen, die zwischen S™ 
und 8, liegt, eine Fläche S, enthalten, die die Fläche & an einem Punkt P 
tangiert!%) (Abb. 3). An diesem Punkt werden also die Flächen S, und & 
— die ja als glatte Flächen vorausgesetzt wurden — denselben Normalvektor 
besitzen müssen. Dieser Vektor muß nun aber als Normalvektor der zeit- 
artigen Fläche S, ein raumartiger, hingegen als Normalvektor von 2 ein 
isotroper Vektor sein. 

Somit führt die Annahme, daß die glaite zeitartige Hyperfläche S mit 
einer glatten Nullhyperfläche 2 keinen gemeinsamen Punkt besitzt, zu einem 
Widerspruch. Wir haben damit unseren Satz bewiesen. 


18) Die weitere Beweisführung gilt auch für den Fall, daß der Punkt Pim Unendlichen 
liegt, d ns also die Fläche S, asymptotisch zu 2 verläuft. 
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Unser Satz erlaubt ohne weiteres noch die beiden folgenden einfachen 
Verallgemeinerungen: 
u 1. Es ist nicht nötig, sich auf überall glatte Nullhyperflächen zu be. 
schränken: Unser Satz gilt auch für (ganze) einheitliche Nullflächen mit 
konischen Punkten. Dies folgt unmittelbar aus der Bemerkung, daß der 
Punkt P auf der Abb. 3 kein konischer Punkt sein 
5 5 S 5m kann (in welchem Fall die Nullfläche an P keinen 
Normalvektor hätte und daher der obige Wider- 
spruch beseitigt sein würde). Denn wäre P ein koni- 
J scher Punkt Pg, so hätte man auch den auf Abb. 4 
gestrichelt gezeichneten Teil von 2 zu berück- 
sichtigen, um die ganze einheitliche Nullfläche 2 
zu bekommen. Dann gäbe es aber notwendig zu 
S™ geodätisch parallele Hyperflächen, die rechts 
von S, liegen und & schneiden, wie z. B. die Fläche 
’ 2 S, auf Abb. 4, im Widerspruch zu der Annahme, 
cu Abb. 4 daß S, die erste Fläche der Schar ist, die 2 
schneidet. 
. An Stelle von glatten raum- oder zeitartigen Hyperflächen können in 
unserem Theorem ersichtlich auch Hyperflächen mit endlich vielen Kanten 
treten, wenn sie nur so beschaffen siud, daß man durch Glättung der Kanten 


bat daraus glatte Flächen erhalten kan , die überall raum- oder überall zeitartig 
sind). 


4. Folgerungen beziiglich der Existenz von zeitlich " 


periodischen singularitätsfreien Feldern 


Wir haben im $ 2 gesehen, daß zeitlich periodische singularitätsfreie Lö- 
sungen der Einsteinschen Gravitationsgleichungen, die die Grenzbedingung 
_ (3) erfüllen, nur dann existieren können, wenn es eine einheitliche Nullhyper- 

fläche 2 gibt, die ganz innerhalb des Bereichs r< a verläuft, also keinen ge- 
 meinsamen Punkt mit der glatten zeitartigen Hyperfläche r =a besitzt. 
u Nach dem im $ 3 hergeleiteten geometrischen Theorem kann es aber keine 
derartige Nullhyperfläche geben. Hieraus folgt unmittelbar, daß es keine 
periodischen singularitätsfreien Lösungen von (1) gibt, welche (3) befriedigen. — 
Da nach Lichnerowicz auch keine nichttrivialen stationären Lösungen 
dieser Art existieren, so folgt, daß die allgemeine Relativitätstheorie kein 
Modell eines beständigen Teilchens, das nur aus dem Gravitationsfeld auf- 
gebaut ist, erlaubt. 


Das im § 3 abgeleitete geometrische Theorem über einheitliche Null- 
= 5 von zeitlich periodischen singularitätsfreien Lösungen der Einstein-Max- 
wellschen Feldgleichungen 
| Ruy = — For Pet, 
19) Eine Nullhyperfläche mit Kanten setzt sich notwendig aus verschiedenen einheit- 


; lichen Nullflächen zusammen. Eine Glättung der Kanten derart, daß dadurch eine glatte 
und somit einheitliche Nullfläche entsteht, ist daher unmöglich. 
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Hierbei ist F,,, = A,,,—A,,, und A, das elektromagnetische Potential. 
Für g,, wird wiederum die Grenzbedingung (3) gefordert. Ferner sollen auch 
hier 9„, und A, wesentlich periodisch sein und nicht nur scheinbar von der 
Zeit abhängen. D. h. es soll nicht möglich sein, durch Koordinaten- und Eich- 
transformationen die Größen g,, und A, auf eine zeitunabhängige Form zu 
bringen. 

Bei diesem Problem wird man ähnlich wie im Fall des reinen Gravitations- 
feldes vorgehen: Nach Papapetrou”) kann auch das kombinierte Einstein- 
Maxwellsche Feld in dem räumlich unendlichen Bereich r > a, in dem das 
Gravitationsfeld als Folge der Grenzbedingung notwendig schwach ist, nicht 
von der Zeit abhängen. Eine periodische Zeitabhängigkeit des kombinierten 
Feldes wäre also nur in demjenigen Gebiet des Innenbereichs r< a möglich, 
in dem das Gravitationsfeld stark ist. Hieraus folgt aber wiederum, daß es 
eine Hyperfläche X geben muß, auf der gewisse Ableitungen von g,, und A, 
Sprünge besitzen. Weil eine wesentliche Periodizität von g,, und A, be- 
stehen soll, dürfen diese Sprünge nicht durch Koordinaten- oder Eichtrans- 
formationen eliminierbar sein. 

Die Forderung, daß die Feldgleichungen (12) im ganzen Raum, also auch 
auf der Sprungfläche 2, gültig sind, ist nun wieder eine stark einschränkende 
Bedingung für die Möglichkeit von nichtforttransformierbaren Sprüngen. 
Für den Fall, daß die niedrigsten unstetigen Ableitungen von g,, und A, 
die zweiten oder höheren Ableitungen sind, hat Stellmacher!®) gezeigt, 
daß auch hier, wie im Fall des reinen Gravitationsfeldes, die Sprungfläche Y 
eine Nullfläche sein muß. Nach Papapetrou und Treder!®) sind nicht- 
forttransformierbare Unstetigkeiten der ersten Ableitungen von g,, und A, 
ebenfalls nur an Nullflächen möglich. Wie beim reinen Gravitationsfeld 
diirften schlieBlich auch im Fall des kombinierten Feldes Unstetigkeiten der 
Potentiale g,, und A, mit den Gleichungen (12) und der Forderung der Regu- 
larität der Metrik nicht vereinbar sein”*). 

Die Sprungfläche X muß auch wieder eine einheitliche Nullfläche sein. 
Dies ergibt sich unmittelbar daraus, daß auch in der kombinierten Einstein- 
Maxwellschen Theorie ein skalarer Erhaltungssatz für die Sprungerregung 
des Gesamtfeldes 


(13) 

er 


In diesen Formeln haben y,, und p, dieselbe Bedeutung wie zuvor, während 
gy, den Sprungvektor von A, bedeutet. E-+4,J ist wieder positiv definit 
in dem Sinne, daß aus E+4J = 0 das Verschwinden der nichtforttrans- 
formierbaren Teile von y,, und , folgt. Aus (14) ergibt sich wieder, genau- 
so wie früher aus (9), daß die Sprungerregung auf einer auf &’ liegenden geodä- 
tischen Nullinie nirgends verschwindet, wenn sie an einem Punkt dieser Linie 


20) A. Papapetrou, Ann. Physik (7) 1, 186 (1958). 

2!) Wie man sich im Rahmen der speziellen Relativitätstheorie leicht überlegt, hat 
eine nicht forttransformierbare Unstetigkeit von A, zur Folge, daß auf der Sprungfläche 2 
nicht nur die elektromagnetische Energiedichte, sondern auch die Gesamtenergie des 
elektromagnetischen Feldes unendlich groß wird. Eine unendliche Gesamtenergie wird 
aber in der allgemeinen Relativitätstheorie notwendig zu nichtregulären g.» führen. 


as 
5 | | 
| 
4 
®B 


= 
Ef 


ungleich Null ist. Die Sprungfläche & muß also jede auf ihr liegende geodä- 
tische Nullinie als Ganzes enthalten und somit eine einheitliche Nullfläche sein, 

Da diese einheitliche Nullfläche 2 wieder ganz innerhalb des Bereichs 
r< a verlaufen muß, darf sie mit der glatten und gemäß der Grenzbedingung 
(3) zeitartigen Hyperfläche r = a keinen gemeinsamen Punkt besitzen. Nach 
dem Theorem des $ 3 ist dies aber nicht möglich. Es gibt somit keine zeitlich 
periodischen singularitätsfreien Lösungen der Einstein-Maxwellschen 
Feldgleichungen (12), die die Grenzbedingung (3) erfüllen. 


Es bleibt noch die Frage, ob es stationäre singularitätsfreie Lösungen der 
Einstein-Maxwellschen Feldgleichungen gibt, die die Grenzbedingung (3) 
erfüllen. Diese Frage wurde aber schon durch die Arbeiten von Thiry und 
Lichnerowicz™) negativ beantwortet. Wie in der Einsteinschen Theorie 
des reinen Gravitationsfeldes existieren demnach auch in der Einstein- 
Maxwellschen Theorie des kombinierten Gravitations- und elektromagneti- 
schen Feldes keine Lösungen, die Modelle von beständigen Teilchen darstellen 
können. 

Das Theorem des $ 3 entscheidet übrigens auch dann über die Existenz 
von zeitlich periodischen Lösungen der Einsteinschen Gravitationsglei- 
chungen oder der Einstein-Maxwellschen Gleichungen des kombinierten 
Feldes, wenn im Innenbereich r< «a periodisch bewegte Quellen des Feldes 
(in Form von Singularitäten oder von genügend kleinen Körpern) existieren. 
Auch in diesem Fall gilt nämlich nach Papapetrou”’), daß das Feld im 
Außengebiet r > a zur Erfüllung der Grenzbedingung (3) notwendig stationär 
sein muß. Es muß also eine im Endlichen des gewöhnlichen dreidimensionalen 
Raumes verlaufende Sprungfläche existieren, die das periodisch zeitabhängige 
innere Feld vom stationären singularitätsfreien Außenfeld trennt. 


Da der V, außerhalb des von der Sprungfläche & umschlossenen räumlich 
endlichen Gebiets eine reguläre Metrik besitzt, kann nach dem Theorem des 
$ 4 die Fläche &' — da sie die zeitartige Hyperfläche r = a nicht trifft — keine 
einheitliche Nullfläche sein. Daher wird es bei einem solchen Feld außer den 
sich im Innern von 2’ bewegenden Singularitäten (bzw. kleinen Körpern) 
auch auf 2 selbst flächenhaft verteilte Quellen des Feldes geben müssen. 
Es liegen dann ähnliche Verhältnisse vor wie bei dem aus der Theorie des 
reinen Maxwellschen Feldes bekannten Fall, wo außerhalb einer geschlosse- 
nen dünnen leitenden (metallischen) Schicht ein zeitunabhängiges Feld exi- 
stieren kann, während sich im Innern geladene Körper beliebig bewegen 
dürfen. Für das reine Gravitationsfeld ist ein solches Feld von Marder 
angegeben worden?®). Allerdings handelt es sich dabei um ein axialsymme- 
trisches Feld, welches in der Richtung der Symmetrie-Achse unendlich aus- 
gedehnt ist, die Grenzbedingung (3) also nicht erfüllt. 

Schließt man dagegen jede flächenhafte Quellenverteilung aus, so folgt 
aus der vorigen Überlegung, daß es kein die Grenzbedingung (3) erfüllendes, 
von periodisch bewegten räumlich verteilten Quellen (Singularitäten oder 
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22) Y. Thiry, J. Math. pures et appliquées (9) 30, 275 (1951); A. Lichnerowicz, 
l. c. Anm. 6, S. 227—231. Die in diesem Beweis verwendete zusätzliche Forderung, daß 
die elektromagnetische Feldstärke F,., im Unendlichen mindestens wie r? verschwindet, 
ist als unabhängige Annahme nicht notwendig, da sie aus den Gravitationsgleichungen 
(12a) und der Grenzbedingung (3) folgt. 

23) L. Marder, Proc. Roy. Soc. London (A) 246, 133 (1958). nike es 
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kleinen Körpern) erzeugtes reines Gravitations- oder kombiniertes Einstein- 
Maxwellsches Feld geben kann. Dieses Ergebnis scheint uns die reale 
Existenz von energietragenden Gravitationswellen überzeugend zu beweisen). 

Unsere Resultate schließen natürlich nicht die Existenz von Modellen 
beständiger Teilchen in anderen Feldtheorien aus, in denen z.B. das Ein- 
steinsche Gravitationsfeld (mit symmetrischem g,,) mit einem Feld gekoppelt 
ist, das nicht das Maxwellsche Feld ist. Insbesondere erscheint es wichtig 
zu prüfen, ob in Einsteins Feldtheorie mit unsymmetrischem g,, derartige 
Teilchenmodelle möglich sind. Dies erfordert die Untersuchung der Frage, 


ob Einsteins unsymmetrische Feldgleichungen®) 
| (15) 


zeitunabhängige oder zeitlich Br singularitätsfreie Lösungen g,, 
besitzen. Dabei wäre man geneigt, folgende Verallgemeinerung der Grenz- 
bedingung (3) zu fordern): 


Juv > Nu» und Jur > 0 fir r>o. i (16) 


Es ist nun zwar nicht geklärt, ob es nichttriviale stationäre singularitäts- 
freie Lösungen von (15) gibt, die die Bedingung (16) erfüllen?”). Jedoch 
würden auch diese Lösungen — wenn sie überhaupt existieren — keine physi- 
kalisch verwendbaren Teilchenmodelle ergeben, da für die Gesamtenergie 
eines derartigen Feldes wiederum das Theorem von Einstein und Pauli 
gilt®). Auch in der unsymmetrischen Feldtheorie ist die Gesamtenergie 
eines stationären Feldes nämlich durch 


gegeben, wobei jetzt T), durch Schrödingers?®) Beziehung 
definiert ist und 7',, die Bedeutung 
Pus Bae + $ ger Rox + Fan Jur gv 


hat. Aus den Feldgleichungen (15) folgt dann das Verschwinden des Ma- 
terietensors 7',, und somit aus (17) das Verschwinden der Gesamtenergie 
eines stationären singularitätsfreien Feldes auch in der unsymmetrischen 
Feldtheorie®). 


24) Vgl. auch D. Geissler, A. Papapetrou u. H.Treder, Ann. Physik (7) 2 
344 (1959). 

25) Vgl. A. Einstein, The Meaning of Relativity, London 1951, S. 127—157. 

Pe Vgl. hierzu auch M. Wyman, Canad. J. Math. 2, 427 (1950). 

27) Vgl. M. Lenoir, C. R. hebd. Séances Acad. Sci. 287, 424 (1953). 

22) A. Papapetrou, Physic. Rev. 78, 1105 (1948), vgl. auch A. Einstein, The 
Meaning of Relativity, Princeton 1955, S. 163. 

29) E. Schrödinger, Proc. Roy. Irish Acad. (A) 52 (1948). 

3) Für Teilchen mit verschwindender Ruhmasse gibt es kein Ruhsystem. Sie werden 
daher nicht durch Lösungen der Feldgleichungen darstellbar sein, die die Grenzbedingung 
(3) erfüllen. 
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Dagegen ist es nicht a priori ausgeschlossen, daß die unsymmetrische 
Theorie Einsteins zeitlich periodische singularitätsfreie Lösungen zuläßt, 
die die Grenzbedingung (16) erfüllen. Eine dafür notwendige Bedingung 
‚läßt sich aus unseren obigen Resultaten sofort herleiten: Die Feldgleichungen 
(15) sollten Lösungen zulassen, die schon im Bereich des schwachen Feldes 
zeitlich periodisch sind, im ragen zu der Theorie des reinen symmetrischen 


Fall war. Denn sonst würde uch in der he Theorie eine 
Sprungfläche 2 benötigt, die das Gebiet des zeitunabhängigen von dem des 
periodischen Feldes trennt. Diese Fläche 2 soilte nach den Ergebnissen von 
Lichnerowicz*!) über das Cauchysche Problem in der unsymmetrischen 


Feldtheorie eine Nullhyperfläche der durch 
« 
definierten Riemannschen Metrik 


(18) 


sein. Auf Grund des Theorems des $ 3 würden wir dann aber wieder auf die 

Unméglichkeit eines solchen Feldes schlieBen. 

Wenn zeitlich periodische singularitätsfreie Lösungen der Feldgleichungen 
_ der unsymmetrischen Feldtheorie existieren, dürfte man wohl erwarten, 
daß diese Lösungen für die Weiterentwicklung des Einsteinschen Teilchen- 
problems von Bedeutung sein werden. 


B Herrn Dipl.-Mathematiker W. Prieß danken wir für seine Hilfe bei der 
Anfertigung des Manuskriptes. 


31) A. Lichnerowicz, l. c. Anm. 6, S. 275—288. 


Berlin, Institut für reine Mathematik der Deutschen Akademie der 
Wissenschaften. 


Bei der Redaktion eingegangen am 16. Januar 159. 
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Über die Elektrophorese _ 
in Argon—Quecksilber-Niederdruckentladungen') 


Von. H. Morgenroth 
Mit 16 Abbildungen 


Inhaltsübersicht 
In einer mit Gleichspannung betriebenen Argon— sili: a 
druckentladung hängen die Hg-Konzentration und die Elektronentemperatur 
vom Abstand von der Kathode ab. Die entsprechenden Beziehungen konnten 
formuliert und die gewonnenen Gleichungen durch Intensitäts- und Elek- 
tronentemperaturmessungen im wesentlichen bestätigt werden. In einem ge- 
sonderten Abschnitt wurde als Anwendung der erzielten Ergebnisse die Elek- 
tronenbeweglichkeit für verschiedene Elektronentemperaturen und Queck- 
silberkonzentrationen berechnet. 


Stamm! gi 
Zusammenstellung der in der Arbeit benutzten Zeichen: 00... 
Der obere Index (i) bezieht sich stets auf das i-te Gas. Speziell ist 3 
i=1 Argon, = 2: Quecksilber = 
b,,b, Elektronen-, Ionenbeweglichkeit a 
D® _ Ambipolare Diffusionskoeffizienten subs 
E, Längsfeldstärke 
fu» = FU, Un», Umuv) = 1 fir U,= 0 und weicht im allgemeinen nicht 
viel von diesem Wert ab = = 
i Entladungsstrom 
Intensität einer Linie entsprechend dem Übergang >u 
k,, kg definiert durch (2.3) rt 
n Konzentration der Neutralgasatome 
N,,%, Konzentration der Ladungsträger; Zusatzindex 0: Konzentration in 
der Achse des Rohres 
Qmu» Wirkungsquerschnitt an der Stelle 
R,r  Rohrradius, Abstand von der Achse 2 
= S Qu, (U,)-v-F(v) - dv: Stoßintegral 


U, Elektronentemperatur in eV 

Potentialdifferenz zwischen den Niveaus uw und v. Für das Ionisations- 

potential steht auch U® an Stelle von Ufo. 

Umu = Uy, — Ug», wo U, die Anregungsspannung für das Maximum des 
Wirkungsquerschnittes in der Approximationsfunktion Q,, (U,) ist 


aus der Diplomarbeit des Verf. 
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(0) = n,/n. = Trägerverteilung. Ist Y (e) die 


so ist Y, = Y (0)-20:-do ein Faktor der Ordnung 4 
& definiert Such (1.10) 


y Korrekturfaktor des StoBintegrals bei Berücksichtigung von Stufen- 
E prozessen in der Ionisation, siehe (4.5) 

Proportionalitatsfaktor I,, und n 

definiert durch (1.11) 

1. Wurzel der Besselfunktion = 2,405 


r/R b 
mittlere Trägerlebensdauer 
% mittlere Lebensdauer des v-ten Anregungszustandes. Pu 
I. Einleitung 


In einer mit Gleichstrom betriebenen Argon— Quecksilber-Niederdruck- 
entladung wird nach Ablauf einer gewissen Zeitspanne nach dem Zünden der 
Entladung in Kathodennähe eine Stauung von Quecksilber beobachtet. Dazu 
u muß die Hg-Konzentration klein gegenüber der Argonkonzentration sein. 
_ Der stationäre Zustand stellt ein Gleichgewicht zwischen dem in Richtung 
ur Kathode fließenden Quecksilberionenstrom und einem entgegengesetzt 
gerichteten, gleichgroßen Diffusionsstrom neutraler Quecksilberatome dar. 
Zu einer qualitativen Erklärung dieser Erscheinung gelangt man unter Beach- 
tung des Konzentrationsunterschiedes und des Unterschiedes in den Ioni- 
_ sationspotentialen für Quecksilber und Argon. Danach ergibt sich für die 
Quecksilberkomponente ein hoher Ionisationsgrad im Vergleich zu dem des 
& und durch die Driftbewegung der Hg-Ionen im Längsfeld wird ein 


wesentlich größerer Anteil des gesamten vorhandenen Quecksilbers zur 

R stammt von Druyvesteyn?). Diese Arbeit wurde in der Folgezeit nicht 
temperatur in Gasgemischen, jedoch ohne Berücksichtigung des Entmischungs- 

ee der Anregungsverhältnisse in Gasgemischen. Durch Wechselstrombetrieb 
Experimentell wurde die Entmischung in einer Arbeit von Penning‘) 

_ Eine Anwendung der Elektrophorese wurde von Riesz und Dieke®) angege- 


x Kathode hin bewegt wie entsprechend beim Argon. 
Die erste eingehendere theoretische Untersuchung der Elektrophorese 
wesentlich weiter entwickelt. Dorgelo, Alting und Boers®) berechnen 
- durch Erw eiterung der Schottkyschen Diffusionstheorie die Elektronen- 

u : effektes. Unter der gleichen Voraussetzung beschäftigen sich einige Arbeiten 
a mit theoretischen und experimentellen Untersuchungen des Gradienten und 
oder anderweitige VorsichtsmaBregeln wird stets dafiir gesorgt, daB keine 

Elektrophorese auftritt. 
untersucht, in der der von Druyvesteyn für einen Bereich nahe der Kathode 
berechnete lineare Konzentrationsabfall spektroskopisch gemessen wurde. 
ben. Die Autoren nutzten die bei Gleichstrombetrieb auftretende Entmischung 
zur Erzielung eines hohen Reinheitsgrades von Gasen aus. 


5 
2) M. J. Druyvesteyn, nage hy , 255 (1935). ea 
3) H. B. Dorgelo, H. Alting, C 3 Boers, Physica 2, 959 (1935). Pe os 
4) F. H. Penning, Physica 1, 763 (1934). i ng 


5) R. Riesz, G. H. Dieke, J. appl. Physics 25, I, 196 (1954). 
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In der theoretischen Behandlung des Effektes wurde von Druyvesteyn 
nicht beachtet, daB die Elektronentemperatur vom Quecksilbergehalt ab- 
hängt, d.h. nach vollzogener Entmischung eine Funktion des Abstandes 
von der Kathode sein muß. In der vorliegenden Arbeit wird daher versucht, 
die theoretischen Grundlagen des Effektes in dieser Richtung genauer zu 
formulieren. 


II. Theorie der Entmischung 

Unter der daß die Quecksilberkonzentration unab- 


hängig vom Abstand von der Rohrachse ist), ist we Diffusionsstrom durch 
den Rohrquerschnitt 


(2) 
Die Dichte des positiven Quecksilberionenstromes ist 


Nach Integration über den Querschnitt und Division durch —e, ergibt sich 
daraus der Teilchenstrom zur Kathode, der dem Diffusionsstrom (1. 1) gleich- 


zusetzen ist: 


Zur Integration der bereits von Druyvesteyn re rg 
gleichung (1.2) ist die Konzentration der positiven Hg-Ionen ns als Funk- 


tion von n'” darzustellen. 
Unter Vernachlässigung von Volumenrekombination und Stufenprozessen 


ist 


(1.3) 


(i) (i) 
s — A ¢,,VU,e %, 


(1.4) 
4) — 1.99. 108. Om ur 
muv 
als StoBintegral. as 


Um einen expliziten Ausdruck fiir t zu gewinnen, wird | angueemeinen, 
daß die in der Trägerbilanz auftretenden Ionenströme der Argon- und Queck- 
silberionen voneinander unabhängig sind und sich analog zum Fall des ein- 
heitlichen Gases mit Hilfe gewisser ambipolarer Diffusionskoeffizienten 
D“ darstellen lassen. Dann lautet die Trägerbilanz 


(4) po D = > (1.5) 
6) Derebenfallsin der Arbeit von Druyveste yn behandelte Radialeffekt der Elektro- 
phorese wird in der vorliegenden Arbeit nicht beriicksichtigt. 
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‘ 
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wo z die Zahl der von einem Elektron in einer sec bewirkten Ionisationen 


des i-ten Gases ist. Mit 2” =n®"s$%, und mit (1.3) ergibt sich daraus 


1 2 u) 
D . (1.6) 


Von den ambipolaren Diffusionskoeffizienten läßt sich zeigen, daß sie in 
Analogie zum Fall des einheitlichen Gases in der Form 


D® = U, ; (1.7) 

geschrieben werden können. b, stellt dabei die Beweglichkeit der Ionen des 
Fer i-ten Gases unter den gegebenen Bedingungen dar, d.h., daß bei kleiner 


J Konzentration der Zusatzkompouente, im vorliegenden Falle Quecksilber, 
für die Quecksilberionen deren Beweglichkeit in Argon einzusetzen ist. Bei 
geringen Konzentrationsunterschieden zwischen beiden Gasen müssen die 
 riehtigen b,-Werte jedoch experimentell ermittelt werden. 
” Aus den Gln. (1.3), (1.6) und (1.7) folgt fir die Konzentration der Queck- 


gilberionen 
(1.8) 


=N No 80 co 
e \Hı 


 Drückt man in dieser Gleichung n,, gemäß 
i = b, Ex Yır FR? 

- durch den Entladungsstrom ¢ aus und geht damit in die Diffusionsgleichung 

(1.2) ein, so erhält man 


4 Wye fore) \ 19 


Die Integration von (1.9) läßt sich formal in einfacher Weise durchführen, 

wenn man die Elektronenbeweglichkeit b, als Konstante behandelt. Der 
dadurch auftretende Fehler ist aber nur für kleine Konzentrationen der Zu- 
 satzkomponente klein. Die Abhängigkeit von 6, von der Elektronentemperatur 
and der Quecksilberkonzentration wird im letzten Abschnitt noch ausführ- 
licher diskutiert. 

Mit Abkürzungen 


A(x) —A(0) = dx (1.11) 


liefert die Integration von (1.9) 
£. Elektronentemperatur 

Die explizite Berechnung von (1.12) verlangt die Kenntnis der Elektronen- 
temperatur als Funktion vom Abstand von der Kathode. Nimmt man ein- 
schränkend an, daß die nee die Feldstärke und die 
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Elektronenkonzentration nicht von x abhängen und vernachlässigt mar 
weiterhin den Einfluß der metastabilen Atome beider Gase, so läßt sich für 
U,(x) eine Differentialgleichung angeben, die bei Kenntnis eines Randwertes 
graphisch gelöst werden kann. Die Einschränkungen können nur bei kleinen 
Konzentrationen der Zusatzkomponente vertreten werden. Andernfalls müßten 
die Ausgangsgleichungen durch eine entsprechende Anzahl von Gleichungen 
ergänzt werden, wodurch die mathematische Behandlung außerordentlich 
erschwert wird. 


Bestimmungsgleichungen für die Elektronentemperatur sind die Träger- 
bilanzausdrücke der Form (1.8) für die Ionenkonzentrationen beider Gase, 
die durch die für den Fall zweier Gase geltenden Quasineutralitätsbedingung 


miteinander verknüpft sind. Man erhält nach Kürzung von n, und mit (1.12) 
“ 


Hı b, U,b, 


i 


Schreibt man die StoBwahrscheinlichkeiten nach (1.4) aus und setzt ee 2 


a) 40) 44) pe (2) 2) 42) pe 


so führt (2.2) zu 


kzexp|- [A (x) —A(0)] U. | e exp 
Durch Logarithmierung und nachfolgender Differentiation nach x “ak mn 
sich daraus 
ve )ewt- | U} (x) 


(2) + U' (x) = 


und schließlich mit 


— 
~ 
| 
- 
® 
— 
~ 


2 (2 2) 
U, U 
der Ausdruck 
1— & 
at N vol uv? & v2 exp] 
- 3, . | U, l e (x) U, (x) ky N vv) 
1— | 
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Die Funktion g [U,(x)] ist nach einigen Umformungen darstellbar in der 
(1) \2 zu) 
— — odeı 
Ai B: x Mit (2.6) lautet die gesuchte Differentialgleichung für U,(x) endgültig 
at (8) 
kon 


Aus (2.6) geht hervor, daß die Funktion g* für kleine U, nicht viel von 1 
abweicht. Für größere U, nähert sich g* der Null, um von dort sehr rasch zu 


= ee stets positiv. Fiir Fin 


die 
2... geht g* > — oo, daher U, (x) > 0. Aus (2.8) geht hervor, daß es sich dabei 
um die Stelle handelt, an der nur noch Atome des Gases (1), also des Grund- 
gases vorhanden sind: 
ky eco U, | e Mit 


Dies ist bis auf den Faktor / U, genau die Trägerbilanz des Gases (1). 
Eine entsprechende Abgrenzung auf der Kathodenseite findet man nicht, 
da sich U, so lange ändert, wie die Zusatzkomponente n“) sich ändert. Dies 
ist zunächst bis zur Begrenzung der positiven Säule durch die Kathode der } wo 
Fall. Steigt n allerdings über eine gewisse kritische Konzentration hinaus f ma 
an, so wird die Entladung vollkommen von den Atomen der Zusatzkompo- | Be 
nente getragen. Die Voraussetzungen für das Bestehen der Diffusionsgleichung f die 
sind nicht mehr erfüllt und der Entmischungsvorgang bricht ab. Damit f nüy 
_ trifft dann auch die Gl. (2.7) nicht mehr zu, aus der hervorgeht, daß U; =0 | licl 
erst für U, = 0 ist. geg 
Die Gl. (2.7) läßt sich graphisch auswerten, indem man das durch die | ter 
rechte Seite gegebene Richtungsfeld zeichnet. In dieses Richtungsfeld ist 
die gesuchte Kurve derart einzupassen, daß sie durch einen Anknüpfungspunkt } gle 
verläuft, der in Form einer Randbedingung ermittelt werden muß. Dazu muß | Ni 
an einer Stelle des Rohres die Konzentration der Zusatzkomponente bekannt J spe 
sein, was experinentell leicht verwirklicht werden kann. Mittels der so ge- | dr 
 wonnenen Elektronentemperaturkurven findet man nach (1.12) den Konzen- 
trationsverlauf der Zusatzkomponente in der Längsachse des Entladungs- 
-rohres. 
3. Anregung und Intensitätsverlauf 


Unter der Voraussetzung, daß nur Direktanregung erfolgt, gilt für die 
r Intensität einer Linie, die dem Übergang » > u eines Atoms des Gases (i) 


hohen negativen Werten überzugehen. Damit bleibt der Nenner v 7 
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Diese Gleichung stellt die theoretische Grundlage für die Bestimmung _ 
des Konzentrationsverlaufes aus Intensitäts- und Elektronentemperatur- 
messungen dar. Mit einem Proportionalitätsfaktor ö,, folgt aus (3.2) für 
konstantes n, 


(i) | 


7) 


Für kleine Beimengungen dürfte (3.3) trotz der Näherung n, = const gut sein. 
Fiir die Uberlegungen beziiglich der Elektronenbeweglichkeit wird noch 
die Änderung von n) benötigt. Durch Differentiation nach x folgt aus (3.! 


(2) (2) 7(2) 
on (2) 0 (2) U jv owl 7 


an einer Stelle des Rohres die gesuchte Konzentration bereits bekannt ist. 
Mit J,,, und U, an dieser Stelle ergibt sich ö, „. na 


4. Stufenprozesse 


In den bisherigen Uberlegungen sind Stufenprozesse nicht beriicksichtigt 
worden. Bei genügend hohen Gasdichten und Elektronenkonzentrationen muß 
man jedoch mit dem Auftreten solcher Prozesse rechnen. Quantitative 
Berechnungen stoßen auf große Schwierigkeiten, da die Lebensdauern und 
die Anregungs- bzw. Ionisationsquerschnitte der Zwischenniveaus nur unge- 
nügend oder gar nicht bekannt sind. Die folgenden Rechnungen geben ledig- 
lich die Möglichkeit zu entscheiden, ob in einem gegebenen Fall an einer 
gegebenen Stelle des Rohres zur Erklärung der dort gemessenen Elektronen- 
temperaturen Stufenprozesse herangezogen werden müssen. 

Gesucht ist die Erweiterung des Ausdruckes (1.3) unter der Annahme 
gleichzeitig auftretender Direktionisation und Ionisation über ein bestimmtes 
Niveau mw. Volumenrekombination wird nicht berücksichtigt. In diesem 
speziellen Fall reduziert sich der von Ludwig”) angegebene allgemeine Aus- 
druck für die u des v-ten Zustandes im realen Plasma auf den Ausdruck 


N,=NnN, 
b 

n, MT, 

n 1 + Bou 


sd E. H. Ludwig, Z. angew. Physik 5, 377, 421 (1953). 
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Durch den Nenner in (4.2) ist den Stößen 2. Art im w-ten angeregten Zustand 
ee getragen. Setzt man : 


und beachtet 
| U, U ool | U, oo} 

80» Suco = - € exp) U, = beexp | A 
hy = NN, T 1 + 


1+n, T, B, 


a, b und ¢ hängen nur schwach von U, ab. Damit ist auch B,, = er, wo 
u 


das Verhältnis der statistischen Gewichte für die Niveaus 0 und u ist, 


mie konstant und der Klammerausdruck hängt merklich nur von der 
ab: 


Ny» = NN, T Y (N,). (4.4) 


_ Fir kleine Änderungen von n, darf y als Korrekturfaktor behandelt werden, 

der insbesondere an den Betrachtungen in den Abschnitten 1 und 2 nichts 

ändert. y >1 weist auf Stufenprozesse der vorausgesetzten Form hin. 
Entsprechendes gilt auch für Stufenanregungen. Eine nähere Betrachtung 


5. Zusammenfassung 


Die in den Abschnitten 1 und 2 gewonnenen Gleichungen (1.12) und (2.7) 
stellen ein simultanes Gleichungssystem für den Verlauf der Konzentration 
der Zusatzkomponente als Funktion vom Abstand von der Kathode des Ent- 
ladungsrohres dar. Die Lösung der Gleichungen läßt sich mit Hilfe der be- 
ae kannten Ionisationsquerse hnitte der beteiligten Gase sowie der Ionenbeweg- 
ze lichkeiten für einen speziellen Fall (Strom, Sattdampfdruck und Druck des 
Ee _ Grundgases gegeben) graphisch ermitteln. Aus den gewonnenen Abhängig- 
keiten n®9 = n® (x) und U, =U, (x) läßt sich ferner der Gang der Intensität 
Be _ einer Linie entsprechend dem Ü bergang v — u für beide Gase errechnen. Die 
I = entsprechenden Kurven erscheinen zusammen mit den Meßergebnissen im 
experimentellen Teil der Arbeit. 

II. Experimentelle Ergebnisse und Vergleiche 
Die Daten des Entladungsrohres sind der Abb. 1 zu entnehmen. W ährend 
_ des Betriebes wurde in der Umgebung des Rohres eine Temperatur von 75° C 
_ aufrechterhalten. Dadurch wurden Quecksilberniederschläge im Inneren des 
Rohres, insbesondere in den Totvolumina, vermieden. Das Rohr e stellt die 
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k zeigt, dab man den dabei auftretenden Korrekturfaktor mit der in Abschnitt 3 
. os . . 
bereits angeführten Konstanten 6,, zusammenfassen kann. Da es sich in der 
7 a _ vorliegenden Arbeit nur um relative Intensitätsmessungen handelt, interessiert m 
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Verbindung des Quecksilbervorrates mit dem Entladungsrohr her. Der 

Quecksilbervorrat diente zur Abtrennung des Entladungsrohres von der 

Pumpanlage und gleichzeitig zur Erzeugung des gewünschten Quecksilber- 


x 


Abb. 1. Entladungsrohr. A: Anode (Molybdänblechzylinder), K: Kathode (Barium- — 
oxydkathode, Heizstrom etwa 8 A), a: gegabeltes Heizrohr, b: Verbindungsrohr zur 
Pumpe bzw. zum Quecksilbervorrat, I bis V: unbewegliche Zylindersonden 


- 
Abb. 2. Anordnung des*Spektrographen vor dem Rohr. a: Spektrograph, b: Spaltko 
ce: Tubus von 3 cm Durchmesser zur Ausblendung eines kurzen Abschnittes der positiven 
Säule, d: Abbildungslinse, e: Führungsschiene 


Hg 5790/70 5461 4916 4858/47/44 4078 4046 3663/54/50 J 


198 
Ar 4091 4158 


Abb. 3. Die Spektralaufnahme läßt deutlich die Abnahme der Quecksilberintensität und 

die Zunahme der Argonintensität längs des Rohres in Richtung Kathode— Anode erkennen. 

Die Aufnahmen entsprechen Abständen von je 6 cm. Auswertung der Platten am Schnell- 
photometer 


dampfdruckes im Entladungsgefäß. Die Sonden I bis V waren unbewegliche 
Zylindersonden, deren von der Glasisolation freien Teile sich in der Achse 
des Entladungsrohres befanden. 

Die Intensitätsmessungen wurden entweder mittels eines ee er 
oder direkt mit GN und Photozelle durchgeführt. Aus der posi 
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: ane tiven Säule wurde ein kurzes Stück ausgeblendet und auf den Spalt des 
‘Spektrographen bzw. der Photozelle abgebildet. Spektrograph bzw. Photo- 
_ zelle konnten in fester Justierung parallel zum Rohr verschoben werden. 
_ Langs der positiven Säule wurden von einem eingestellten Entladungs- 
zustand jeweils 14 Aufnahmen im Abstand von je 6 cm gemacht. 
en Sondencharakteristiken bestimmt. Andererseits wurde versucht, unter der 
Es Voraussetzung konstanter Argonkonzentration zusätzlich Elektronentem- 
ee ee peraturwerte aus Argonintensitätsmessungen zu gewinnen. Dazu muß aller- 


ela Er dings nach (3.3) die Elektronentemperatur wenigstens an einer Stelle des 
Rohres bekannt sein, d. h., mindestens eine Sondenmessung ist unvermeidlich. 


Die Elektronentemperaturen wurden einerseits nach der Methode der 


Im Laufe der vorbereitenden Messungen zeigte es sich, daß die Füllung 
des Rohres mit Quecksilber beträchtliche Zeit in Anspruch nahm. Brannte 
im Rohr bei zunächst ausgefrorenem Quecksilbervorrat eine reine Argon- 
entladung und wurde zu einer bestimmten Zeit das Kühlgefäß durch das 
_ Wärmebad ersetzt, so konnte bei fortgesetzter Zunahme des Quecksilber- 
gehaltes im Rohr erst nach 4 bis 5 Stunden die erwartete Stationarität nach- 
_ gewiesen werden. Dabei wurden zum Nachweis des stationären Zustandes die 
Spannung, die Grenzstromstärke für den Einsatz von Plasmaschwingungen 
und der Gang der Intensität einer Quecksilberlinie in Abhängigkeit von der 
Zeit gemessen. Die Vorgänge, die zu der genannten unerwartet großen Zeit- 

_ differenz führen, wurden nicht näher untersucht. Eine eingehende Unter- 
7 _ suchung des nichtstationären Zustandes erscheint jedoch lohnend, und zwar 
auch für den Fall, wo nur der Entmischungsv organg selbst beobachtet wird, 
d. h., wo ünschte im Rohr schon ist. 


= Für die a; des TEEN wurden nur die Werte 7, = 30 
En. 50°C mit den dazugehörigen Sattdampfdrucken p, = 2,79 - 103 und 
28.10 mm Hg gewählt. Bei Temperaturen des Quecksilbervorrates 
oberhalb 50° wird die Entladung vollständig vom Quecksilber getragen. Es 
er findet keine Entmischung mehr statt. Die Untersuchung des Effektes bei 
Temperaturen unterhalb von 30° erfordert einen größeren technischen Auf- 
Ea wand und führt qualitativ zu keinerlei neuen Ergebnissen. Eine Unter- 
teilung des Temperaturbereiches zwischen 30 und 50° ist nicht möglich, da 
sich die Elektronentemperaturen bei Änderung von T,= 30 auf T, = 40° 
etwa um 10%, ändern und damit noch in dem bei Sondenmessungen auf- 
_ tretenden Fehlerbereich liegen. 


7. Zusammenstellung der in den Rechnungen benutzten Konstanten 


. Der Quecksilberpartialdruck an der Stelle «= x, (Mündung des Zu- 
 Jeitungsrohres in das Entladungsrohr), entspricht dem Sattdampfdruck des 
 — Quecksilbervorrates bei dessen Temperatur. Da die Temperatur im Rohr 
x an der Stelle x, von der Temperatur des eserves vr ist, 


feat ther wobei für T die Umgebungstemperatur des Rohres einzusetzen 
ist wad be Daraus ergeben sich bei den gewählten Temperaturen 7, des 
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Quecksilbervorrates am Fixpunkt x, die Konzentrationen 
T, = 30° C:n'® (z,) = 7,85 - 10% cm? 
T, = 50° C:n® (z,) = 3,6 - 1014 
Die Argonkonzentration entspricht dem Fülldruck p = 1 mm Hg. Sie 
wurde nicht variiert. Die Füllung erfolgte bei Zimmertemperatur: 
nd) — 3,3 - 1016 em-3, 


Für den konstanten Faktor des Stoßintegrals (siehe 1.4) ergaben Mes- 
sungen an reinem Argon 


= 1,286 - 10-8. 


Der entsprechende Wert für Quecksilber wurde nach Zahlenangaben aus 
der Arbeit von Ludwig berechnet. Danach wird die Stoßquerschnitts- 
funktion Quo (U,) durch Überlagerung zweier Ausdrücke für Quo (U,) mit 
passenden Q,,90 und Umax = U,,90 approximiert. Für den Faktor des 
StoBintegrals ergibt sich daraus 


AS, = 6,95 - 10-8. 


Bei Argon führten ähnliche hai mit Zahlenwerten aus Granowski®) 
zu geringfügigen Diskrepanzen mit den Messungen. Daher wurde der oben 
genannte experimentelle Wert benutzt. 

Die in der Literatur vorhandenen Angaben über die Beweglichkeiten von 
Argon- und Quecksilberionen in ihren eigenen Gasen sind im vorliegenden 
Falle nicht ohne weiteres verwertbar. Dies gilt insbesondere für die Queck- 
silberionen, deren Beweglichkeit wegen des großen Konzentrationsunter- 
schiedes hauptsächlich durch Stöße mit den Argon-Atomen bestimmt wird. 
Daher wurde die Beweglichkeit der Quecksilberionen aus den Messungen selbst 
errechnet, obwohl dies den konsequenten Vergleich der dargestellten Theorie 
mit den Intensitäts- und Elektronentemperaturmessungen stört. Zur Er- 
mittlung aus Messungen genügt die Kenntnis der Elektronentemperatur und 
der Konzentrationen n“) und n‘) an der Stelle x, des Rohres. Aus (2.2) ergab 
sich 

(2) R cm? 

= 1,083 "Vol 
der Literatur, die sich auf die Beweglichkeit der Argonionen in Argon selbst 
beziehen, hinreichend sein, da die Argonkonzentration um mindestens drei 
Größenordnungen höher liegt, als die Quecksilberkonzentration. Holstein?) 
gibt Beweglichkeiten für Atom- und Molekülionen an. Aus zwei neueren 
Arbeiten von Boyd und Morris’) und Boyd") geht hervor, daß in 
dem hier zugrunde liegenden Strom- und Feldstärkebereich nur die Beweg- 
lichkeit der Atomionen in Frage kommt: 


9 3 cm? 
bp — 1,22 10 V 1 
olt - sec 


8) W. L. Granowski, Der elektrische Strom im Gas. Ubersetzung: Berlin 1955. 
%) T. Holstein, Physic. Rev. 82, 567 (1951). 

10) A. L. F. Boyd, D. Morris, Proc. physic. Soc. A 68, 1 (1955). 
1) A. L. F. Boyd, Proc. physic. Soc. A 68, 11 (1955). 


x H. Morgenroth: Uber die Elektrophorese in Argon—Quecksilber-N iederdruckentladungen 383 7 
» 
. 
a 
| 
F 
f 
5 
' 
j 
Ä 
4 
. : 2 


Annalen der Physik. 7. Folge. Band 3. 1959 4 


Da die Elektronenbeweglichkeit von U, und n2 abhängt, darf sie eigentlich 
nicht als Konstante behandelt werden. Die Theorie des Elektronentem- 
u peratur- und Quecksilberkonzentrationsverlaufes wurde jedoch unter der 
Voraussetzung b, = const entwickelt. Diese Voraussetzung trifft zu für kleine 
Konzentrationen der Zusatzkomponente n) und für geringe Schwankungen 
der Elektronentemperatur. Da zunächst nicht sicher war, unterhalb welcher 
Konzentration der Einfluß des Quecksilbers auf die Elektronenbeweglichkeit 
zu vernachlässigen ist, wurden den Berechnungen des folgenden Abschnittes 
(T, = 30°) versuchsweise zwei Werte zugrunde gelegt. Der eine dieser 
Werte ergibt an der Stelle x = x, die gemessene Neigung der Elektronen- 
temperaturkurven: | 


Er bezieht sich auf n@(x,) = 7,85- 10% und U, = 14000° K. Der zweite 
tee bezieht sich auf die gleiche Konzentration und auf U, = 15000° K: 


m? 
b,, = 2,9 106 


Volt - 
Er liefert in einem mittleren Bereich des Rohres die richtige Neigung der 
Quecksilberintensitätskurve für i= 3A. Die Werte b, und 6,, reichen 
aus, um die Theorie durch Vergleich mit den Messungen prinzipiell zu be- 
stätigen. Der Fall 7, = 50° kann damit jedoch nicht mehr erfaßt werden. 
In dem entsprechenden Abschnitt wird daher nur qualitativ diskutiert. 
Dagegen eröffnet sich hier die Möglichkeit einer Ermittlung der Funktion 
b, = b, (n™, U,). Die Ergebnisse dieser Betrachtungen bestätigen die An- 
nahmen 6,, und b,.. 
Der Diffusionskoeffizient für die Diffusion der Hg-Atome in Ar ergibt 
nach einer Formel von Stefan-Maxwell [Eucken 12), S. 344] für den 


oe vorliegenden Fall zu 


D = 75 em?/sec. 


8. Resultate bei T, 
81 Verlauf der Elektronentemperatur für verschiedene Strom- 
stärken 
Um in bezug auf die Messungen eine Diskussionsgrundlage zu haben, 


sollen an dieser Stelle zunächst die sich aus (2.7) ergebenden Abhängigkeiten 
 U,(x) fir i= 2, 3,4 und 5A erörtert werden. Die Elektronentemperatur 
am Fixpunkt x, ist zur Bestimmung der Beweglichkeit der Quecksilberionen 


bereits experimentell zu U,(x,) = 14000° K ermittelt worden. Dieser Wert 


Nach (2.7) ist der Anstieg dem Strom proportional. Steilere Kurven 
_ entsprechen also höheren Stromstärken. Es zeigt sich, daß alle Kurven auf 


12) A. Eucken, Lehrbuch der chem. Physik II, 1 (1950). IE 


he: 


, 


der 
pera 
einn 
Tray 
Elel 
höh 
trat 
aus 
der 
Str 
mis 
wal 
ein 
pul 
mit 
len 
bei 
ni? 
ste 
rat 
tra 
| ne! 
St 
dr ret 
ra’ 
lac 
ell 
de 
E 
ei 
T 
m 
hi 
hi 
li 
infolge der konstanten Quecksilberkonzentration an dieser Stelle unab- 
.  hängig vom Entladungsstrom. Alle aus (2.7) zu ermittelnden Kurven müssen u 
daher den Punkt (14000°, x,) durchlaufen. Man gewinnt sie, indem man für h 
einige U,-Werte mit Hilfe der besprochenen Konstanten den jeweiligen An- 
stieg von U, berechnet. Das Resultat ist in Abb. 4 wiedergegeben. r 
¢ 
t 


ntlich 
ntem- 
r der 
kleine 
ungen 
elcher 
chkeit 
nittes 
dieser 
‘onen- 


g der 
ichen 
u be- 
rden. 
tiert. 
ktion 


H. Morgenroth: Uber die Elektrophorese in Argon—Quecksilber-N iederdruckentladungen 385 


der Anodenseite in die Gerade U, = const = 15000° K (Elektronentem- 
peratur für reines Argon bei p = 1 mm Hg, unabhängig von der Stromstärke) 
einmünden, und zwar um so eher, je größer der Entladungsstrom ist. Da laut 
Trägerbilanz der niedrigeren 

tration entspricht, wird, wie 
aus Abb. 4 hervorgeht, von 
der Theorie für höheren 
Strom die stärkere Ent- 
mischung verlangt. Man er- 
wartet also experimentell an KER® Xs 30 36 42 48 Su 60 66 72 BA 


‘ [cm): 
Abb. 4. Theoretischer Verlauf der Elektronentempe- 
(Kathodenseite) atur. A: Anode, K: Kathode, x,: Mündungsstelle 
mit steigendem Strom fal- der Quecksilberzuführung in das Entladungsrohr 
lende Elektronentemperatur 


bei steigender Konzentration 
steigende Elektronentempe- 77000 
ratur bei fallender Konzen- 76000 + E 
tration der Zusatzkompo- 15000 
nente. 14000 

Da an einer beliebigen 33001 
Stelle des Rohres die theo- ot 
retischen Elektronentempe- 
raturwerte für alle Ent- 
ladungsströme innerhalb K 6 2 ®B xz 30 36 42 48 54 60 65 72" 78 A 
eines 10proz. Bereiches [em] 


liegen, läßt sich bezüglich Abb. 5. Verlauf der Elektronentemperatur nach 
der Stromabhängigkeit der Sondenmessungen. Theoretischer Verlauf für i=3A 


Elektronentemperatur an 

einer Stelle des Rohres die 

Theorie nicht durch Sonden- 17000 + fee) ? pom 
messungen bestätigen, da 1600} 
hierbei alle Meßwerte inner- 15000 + 
halb der Fehlergrenzen 14000 H 
liegen, die fiir Sondenmes- 13000 + 
sungen veranschlagt werden 1 2 

müssen. Die Verhältnisse 

gestalten sich erst bei höhe- K6 26 x 30 36 60 66 72°78 

rem Quecksilbergehalt giin- 


stiger, wie die Messungen {ith cssungen, berechnet nach (3.3) unter der Vor. 
bei T, = 50° zeigen. aussetzung, daß die Elektronentemperatur an der 
Abb. 5 belegt die grund- Stelle x = x, bekannt ist. Theoretischer Verlauf für 
sätzliche Übereinstimmung i=3A gestrichelt 
der Abhängigkeit der Elek- 
tronentemperatur vom Abstand von der Kathode mit dem entsprechenden 
theoretischen Verlauf jedoch auch für 7, = 30°. Grobe Abweichungen zeigen 
sich an Sonde V. Hier steigen die Elektronentemperaturwerte beträchtlich 
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über den Wert hinaus an, der aus der Trägerbilanz für reines Argon gefordert 
wird. Die Ursache für diese offenbar systematischen Abweichungen konnten 
nicht geklärt werden. Der Versuch einer mathematischen Erfassung der 
Ortsabhängigkeit von Stufenprozes- 


sen, Ionen- und Elektronenbeweg- 
_ lichkeiten erforderte einen hohen 
0 Aufwand, der nur zu geringen Korrek- 
8r turen führte. 
Abgesehen von den Abweichungen 
of 4r an Sonde V sind die aus der Rück- 
rechnung aus Argon-Intensitätsmes- 
oh sungen gewonnenen Ergebnisse noch 
je. befriedigender (Abb. 6). Eine Abhän- 
: 4 wt gigkeit vom Strom ist wegen der Ge- 
Pay ringfügigkeit des Effektes allerdings 
; auch hier nicht erkennbar. 


8.2 Hg-Konzentration und In- 


Abb. 7. Die den Konzentrationsberech- tensität der Hg-Linien füri—=3A 


nungen zugrunde liegende Integralfunktion Der Verlauf der Quecksilber- 


A(x) für die durch Sondenmessungen (S) 
und durch Argonintensitätsmessungen (A) konzentration wird grundsätzlich 


gewonnenen Elektronentemperaturkurven durch (1.12) wiedergegeben. Zur Aus- 


K 612% 2% 30 36 42 48 54 60 66 72 78 84 A 
lem] 


fir i= 3A wertung dieser Formel muß zunächst 
a ee der Verlauf der Integralfunktion A (x) 
9 Bu, au diskutiert werden. Neben dem theo- 


retischen Verlauf der Elektronen- 
temperaturkurve fiir ©=2A sind 
dabei wegen der Abweichungen auch 
die gemessenen Elektronentempe- 
raturkurven zu beriicksichtigen. Der 
aus Sondenmessungen gewonnene 
Verlauf von U, (x) wurde hierzu mehr 
oder weniger willkürlich geglättet. Die 
Integration wurde graphisch durch- 
geführt. Die Ergebnisse sind in 
Abb. 7 wiedergegeben. Es zeigt sich, 
daß die Schwankungen der anoden- 

seitig gemessenen Elektronen- 

temperaturen durch die Inte- 

gration weitgehend ausgeglichen 

werden. Ein nichtmonotoner Ver- 

lauf von U,(x) geht in einen 


ll monotonen Verlauf der A-Funk- 
KG 2 18 x5 30 3642 48 960 66 72 2 = tion über. Wie im folgenden Ab- 


schnitt noch gezeigt wird, be- 
6 hnitt h gezeigt wird, b 

 Abb.8. Quecksilberkonzentrationsverlauf nach (deutet dies. daß der prinzipielle 

= den Elektronentemperaturmessungen ($) und Verlauf d K trati 

(A) fir i = 3A. Entsprechender theoretischer YeTlaul der Konzentration gegen- 

Verlauf gestrichelt über solchen Elektronentempe- 
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raturschwankungen unempfindlich ist. Der Intensitätsverlauf reagiert dar- 
auf stärker. 

Mit Hilfe der A-Funktionen läßt sich nun ohne weiteres n(2) (x) nach (1.12) 
bestimmen, wenn man beachtet, daß die Quecksilberkonzentration an der 
Stelle x, bekannt ist. Das Ergebnis der Rechnungen ist in Abb. 8 wieder- 
gegeben. Infolge der ausgleichenden Wirkung der Integralfunktion unter- 
scheiden sich die Konzentrationsverläufe trotz der beträchtlichen Abwei- 
chungen der Elektronentemperatur- 499 
kurven nicht sehr stark voneinander. oL 
Ob die durchschnittliche Steilheit des 
KonzentrationsverlaufesderWirklich- dei 
keit entspricht, kann nur der Ver- Intensität 


gleich mit den daraus ermittelten 
Intensitätsabhängigkeiten zeigen. ar ’ 


Die mit der A-Funktion errech- aut ite 
neten Intensitätskurven sind zu- „oL 
sammen mit dem für ©=3A ge- 
messenen Intensitätsverlauf in Abb. 9 


aufgetragen. Es zeigt sich, daB die i 


mit Hilfe des Elektronentemperatur- sj 

verlaufes A ermittelte Intensitäts- Qe 
kurve gut mit der gemessenen über- rein theoretisch 
einstimmt. Bis zur Stelle x = 78 Boe tie 

findet man sogar in Feinheiten Uber- 


einstimmung. Man beachte z.B. die  08- 
Ausdellungen an den Stellen 2 = 36 06F 
und 72. Sie spiegeln das Verhalten 


diesen Stellen deutlich wider. Dar- K 6 12 18 x 30 36 42 48 54 60 66 72 78 BUA 
aus geht hervor, daB der nach Argon- [em] 


intensitätsmessungen ermittelte Ver- Abb. 9. Quecksilberintensität nach den 
lauf der Elektronentemperatur der Elektronentemperaturen (S) und (A). Zum 
Vergleich sind der rein theoretische sowie 
W irklichkeit besser entspricht, als der gemessene Intensitätsverlauf einge- 
die groben Elektronentemperatur- zeichnet. Die gemessene Intensität wurde 
bestimmungen durch Sondenmes- an der Stelle x = 6 willkürlich = 100 ge- 
setzt. Die berechneten Intensitäten wurden 

sungen. 

der Stelle r= der dort gemessenen 

& Intensität angeglichen 
theoretisch bestimmten Intensitäts- 
verlaufes mit dem gemessenen ist ebenfalls bedriedigend. Damit sind prinzipiell 
die Formeln (1.12) und (2.7) bestätigt. Ferner ist der Zusammenhang zwischen 
Konzentration und Intensität richtig (3.2). Insbesondere ist dann auch der 
nach A, Abb. 7, mit Hilfe von (1.12) ermittelte Konzentrationsverlauf A, 
Abb. 8, als der Wirklichkeit entsprechend anzusehen. 

Die Abweichungen von den rein theoretischen Kurven in Kathoden- und 
Anodennähe des Rohres sind auf Unsicherheiten im experimentellen Elek- 
tronentemperaturverauf in diesen Bereichen zurückzuführen. Ferner ist zu 
bedenken, daß die Elektronenbeweglichkeit als Konstante behandelt wurde. 
In Kathodennähe z. B. wird b, wegen der dort herrschenden größeren Queck- 
silberkonzentration kleiner sein, wodurch ein etwas steilerer Anfangsabfall 
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erforderlich. 


a 8.3 Hg-Konzentration und Intensität der Hg-Linien bei verschie- 
denen Stromstärken 


d “ auftritt, gewinnt man folgendermaßen: 

u 3 

! 

2 


4 oe) 
2 


K 36 U2 48 5h 60 66 72 78 84 A 
(@ 
Abb. 10. Konzentrationsverlauf in Abhängig- 
rope von der Stromstirke. An der Stelle 
- 84 beträgt bei 5A die Konzentration 
10%. Der Konzentrationsunterschied 
zwischen Kathoden- und Anodennähe beträgt 
damit mehr als 3 Größenordnungen 


der Konzentration und damit der Intensität bewirkt wird. 
ist die Quecksilberkonzentration so klein, daß die Beweglichkeit nur noch 
von der Elektronentemperatur abhängt. 
abschnitt offenbar schwankt (nach den A-Kurven um wenigstens 1000°), 
wären auch in der Konzentrations- und Intensitätsberechnung Korrekturen 


In Anodennähe 


Da diese in dem genannten Rohr- 


Um einen Überblick über die Abhängigkeit der Konzentrationskurven 
vom Strom zu gewinnen genügt es, sich auf einen Elektronentemperatur- 
verlauf zu beschränken. Eine Abschätzung des Fehlers, der durch die Nicht- 
berücksichtigung der Stromabhängigkeit der Elektronentemperaturkurven 


Nach der theoretisch ermittelten 
Kurvenschar Abb. 4 fällt U, links 
von x, und steigt rechts von x,, 
wenn der Entladungsstrom zu- 
nimmt. Die Integration der ein- 
zelnen Kurven einer solchen 
Kurvenschar hebt die genannte 
Tendenz jedoch auf. In Anoden- 
nähe ist daher A (x) weitgehend 
unabhängig von i. Der Konzen- 
trationsverlauf wird dort durch 
einen beliebigen Verlauf von’ A (x), 
etwa durch den für «=3A, und 
durch die neben A(x) im Expo- 
nenten von (1.12) noch auftretende 
Stromstärke bestimmt. In der 
Kathodenhäliie macht sich der 
Einfluß der Abhängigkeit von A (x) 
von i etwas stärker bemerkbar. 
Höhere Ströme haben hier einen 
steileren Anstieg von U,(x) zur 
Folge. Da die ausgleichende Wir- 
kung der Integration für kleine x 
fehlt, wird mit einer für alle Ströme 
einheitlich angenommenen Funk- 
tion U, (x) bei höheren Strömen 
ein etwas zu schwacher Konzen- 
trationsabfall errechnet. Der Feh- 
ler ist jedoch wegen der nach 
Theorie und Experiment offenbar 


geringen Änderung von U,(x) mit dem Strom auch hier vernachlässigbar 


klein. Vergleichende Rechnungen haben gezeigt, 
_ ponenten von (1.12) im Falle einer geringen Gesamtänderung der Elektronen- 
temperatur über die Rohrlänge stets den stärksten Einfluß besitzt. Da- 
, daß die Unterschiede des Elek- 


daß der Strom im Ex- 


oo psp 


| 
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»nnähe | tronentemperaturverlaufs bei verschiedenen Strömen nicht eindeutig er- 
r noch | kennbar waren. In Abb. 10 sind nach (1.12) mit Hilfe des fiir i= 3A ge- 
Rohr- | messenen Elektronentemperaturverlaufes berechneten Konzentrationsab- 
1000°), | hängigkeiten für i = 3, 4 und 5A dargestellt. Nach der angeführten Fehler- 
kturen | abschätzung dürften sie der Wirklichkeit gut entsprechen. Sie sind den im 
folgenden angeführten Intensitätsberechnungen zugrunde gelegt. 

Die Intensitäten für 3, 4 und 5 A ergeben sich aus Abb. 10 im wesentlichen 


U, v\ 


schie- | durch Multiplikation mit dem Faktor exp|- uf Wegen dieses Faktors 


reagiert die Intensität wesentlich stärker auf Elektronentemperaturschwan- 
curven | kungen, als die Konzentration. Unterscheiden sich die Elektronentem- 
ratur- | peraturkurven für die einzelnen Stromstärken so, wie dies aus Abb. 4 hervor- 
Nicht- | geht, oder stärker, so muß der entsprechende Einfluß auf die Intensitäts- 
curven | abhängigkeiten bemerkbar sein. Die Abbn. 11 und 12 zeigen einmal die ge- 


ren | ar koi 
olchen 60F \\, Intensität 

ehend 20 NN \ 
onzen- 

% © 8 
A(x), 10 |+- * 
, und 8 1 \. 
Expo- 6 


nm 


‘tende 4b \ 
ı der x= 3 Amp AN 
h der 4= 4Amp N 4 


n A (x) o= 08 

rkbar. \ 06 + 

c) zur 08 04 

K6 36 4248 54 60 66 72 78 B4A K6 1218 36 42 48 60.66 72 79 

[cm] cm] 

a Abb. 11. Sp2ktrographisch ermittelter In- Abb. 12. Intensitätsverlauf für 1 = 3, 4 und 

a tensitätsverlauf der Quecksilberlinie 5461A 5A, berechnet nach (3.2) mit n) nach 

‘omen für 3, 4 und 5A. Die Intensität an der Abb. 10 und U, nach (A). Die freie Kon- 

nzen- Stelle 2=6 ist für jede Stromstärke stante ist jeweils so gewählt, daß die be- 

. Feh- willkürlich = 100 gesetzt rechnete Intensität im Fixpunkt x, mit 
h der dort gemessenen Intensität überein- 

nac stimmt. (Vgl. Abb. 11) 

enbar 

sigbar messenen Intensitätskurven und zum anderen die Intensitätskurven, die sich 

| Ex- aus den Konzentrationsabhängigkeiten nach Abb. 10 unter der Annahme eines 

onen- einzigen Elektronentemperaturverlaufes ergeben. Die Übereinstimmung 

. Da- ist gut. Es zeigt sich jedoch, daß die berechneten Kurven stärker gespreizt 

Elek- sind, als die gemessenen. Tatsächlich muß der einem höheren Strom zu- 
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kommende steilere Verlauf der Elektronentemperaturkurve einen flacheren 
Verlauf der Intensität bewirken, wodurch die Parallelität der gemessenen 
4- und 5-A-Kurven bezüglich der 3-A-Kurve hervorgerufen wird. Da die 
_ Bteomabhängigkeit der Elektronentemperaturkurven nicht eindeutig be- 
stimmt werden konnte, muß auf eine eingehendere quantitative Diskussion 
verzichtet werden. 


9. Resultate bei T, = 50° 


9.1 Elektronentemperaturen 


Die aus Sondenmessungen gewonnenen Elektronentemperaturen sind in 
_ Abb. 13 wiedergegeben. Es fällt zunächst auf, daß U, (x,) vom Strom ab- 
hängt. Dies läßt sich durch das Auftreten von Stufenprozessen bei der Ioni- 

sation des Quecksilbers erklären. Der dafür maßgebliche Faktor y läßt sich 
wie folgt ermitteln: Zunächst kann 


| U, (x,) fiir y = 1 nach 
a) 2 
os 5 Amp ? nit) 2820 e fol . bed % 
20000 &= 4 Amp db \m 
r x = 3 Amp 
18000} += 2 Amp mi 
4 o= 1 Amp 9 +n =] 
4 
oo 
12000 H- — (10.11) 
5 + mit 
10000 £ a) 40) fico R R 
9000 F — = 2,42 . 
(1) ’ > 
K 6 12 18 xz 30 3642 48 54 60 66 72 7884 A n (2,) AD, 
Abb. 13. Elektronentemperatur nach Sonden- My b, 


 messungen für 7, = 50°. Die Abweichungen 
an Sonde V konnten nicht erklärt werden bestimmt werden. Man erhält 


U, = 1,065 eV = 12350° K. 
In der Nähe dieses Wertes liegen nach Abb. 13 nur die Kurven U, (x) für 
 i=05A und 1A. Die Kurven füri >1 A passieren die Stelle x, bei 10000 
und 11000°. Man hat also für diese Elektronentemperaturen (10.11) zu 
lösen, wobei sich neue k,-Werte ergeben, die dann in die Faktoren k, - y 
_ aufgespalten werden können. Man erhält 


ig @et.- U, (°K) U, (x,) (eV) k, 
10500 | 0,905 | 8,36 - 10 | 5,22 (10.12) 
11000 0,948 | 5,33 - 104 
12350 1,065 1,6 - 108 1 


Die Annahme solcher Stufenprozesse fiir Argon ist unwahrscheinlich, da 
Argonatome eine für metastabile Niveaus besitzen, die 
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höher liegt, als die des Quecksilbers. Ferner ist uf 
(v metastabil), klein gegen Ue. » so daB durch Stufenionisation keine wesent- 
lich günstigeren Verhältnisse entstehen würden. Dies läßt sich zeigen, indem 
man die Argonionisationsenergie gleich der Anregungsenergie metastabiler 
Stufen des Argons setzt (~ 12V). Nimmt man noch gleichen Stoßquer- 
schnitt an (ungünstigste Annahme), so wird (10.11) gelöst durch 


Die für i = 3, 4 und 5A gemessenen Elektronentemperaturen liegen aber 
noch darunter. 

Ebensowenig können Stöße zwischen Argonatomen im metastabilen 
Zustand und Quecksilberatomen zu einer Herabsetzung der Elektronen- 
temperatur führen, da dies eine hohe, mit der Elektronenkonzentration ver- 
gleichbare Konzentration metastabiler Argonatome erfordert. Zur Erklärung 
der gemessenen Elektronentemperaturen bleibt daher nur noch die Annahme 
von Stufenionisationen der Quecksilberatome. 

Außer an Sonde V entspricht im übrigen die Kurvenschar U,(x) quali- 
tativ den Erwartungen. Auf Grund des höheren Quecksilbergehaltes im Gegen- 
satz zu 7’, = 30° ist hier deutlich eine Abhängigkeit vom Strom erkennbar. 

Die außerordentlich hohen Elektronentemperaturen an Sonde V haben 
keine Erklärung gefunden. Sie stehen im Widerspruch zu den aus Argon 
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Abb. 15. Intensitätsverlauf der Linie 5461 A 
fiir einige Stromstirken. An der Stelle 
x = 6 ist die Intensität fiir i = 5 A willkür- 
lich gleich 100 gesetzt. Die eingezeichneten 
Intensitätswerte für i= 4, 3 und 2A stehen 
im richtigen Verhältnis zu denen bei 5A 


K 6 2 8 36 42 Sh 60 6672 78 


Abb. 14. Elektronentemperaturverlauf 
fir «= 5A aus Sondenmessungen (S) 
und Argonintensitä De (4) 


intensitätsmessungen gewonnenen Werten. Dazu Abb. 14. Eine quantitative 
Diskussion an Hand von (2.7) ist wegen der erwähnten Stufenprozesse sowie 
wegen der zunehmenden Abhängigkeit der Elektronenbeweglichkeit von der 
Quecksilberkonzentration in diesem Rahmen nicht — 


q i 
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— q 
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9.2 Intensität der Hg-Linien 
R ‘Die Intensitätskurven der Abb. 15 sind mit Hilfe einer Photozelle aufge- 
nommen. Die wiedergegebenen Intensitätswerte für verschiedene Ströme an 
gleichen Stellen des Rohres stehen daher im richtigen Verhältnis zueinander. 


Die Messungen entsprechen qualitativ den theoretischen Erwartungen. 

Mit wachsendem Strom nimmt der Entmischungsgrad zu. Gleichzeitig nimmt 

£ U, im Durchschnitt zu, so daß die Intensitätskurven in Anodennähe parallel 
verlaufen. Bei Strömen kleiner als 2A konnte keine Intensitätsänderung 
mehr werden. Daraus darf aber nicht geschlossen werden, daß 


n ar = (ist. Auch für i = 0,5 A geht aus Abb. 13 noch eindeutig U, + 0 


wo "Damit ist dann auch die Konzentrationsänderung des Quecksilbers 
von Null verschieden. 


Man erhält den ungefähren Verlauf der Konzentration für 7, = 50° 
durch Verschiebung der entsprechenden Kurven für 30° (Abb. 10) derart, 
daß diese an der Stelle «= x, der Wert n® (x,) = 3,6 - 101% durchlaufen. 
ae Bei den gemessenen Intensitätskurven der Abb. 11 und 15 ist ebenfalls nur 
eine geeignete Verschiebung erforderlich, um entsprechende Kurven zur 

Deckung zu bringen. Die der Teile der Intensitäts- 


IV. Elektronenbeweglichkeit 


In einer Arbeit von Wojaczek und Rademacher”) werden verschiedene 
Formeln für die Elektronenbeweglichkeit durch Messungen an reinem Argon 
überprüft. Die Meßwerte genügen für 1Torr Argon mit guter Überein- 

stimmung der dort angeführten Formel (8): 


Für zwei Gase mit den Konzentrationen n() und n(2 und den Querschnitten 
co und of) läßt sich daher schreiben 


bar 
€ 1/b, m f In + n® v F(v) dv, 


wobei - 


= 1/0 + 1/6 
1/6, 1/b + 1/b, yer os (IV. 2) 


13) K. Wojaczek u. K. Rademacher, Ann. Physik 18, 237 (1956). bh ch 


_ gesetzt wurde. Daraus folgt unmittelbar 
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d.h., b, kann in Einzelbeweglichkeiten zerlegt werden derart, daß sich die 
reziproken Beweglichkeiten bezüglich der einzelnen Gase addieren. Nun 
kann man die Funktion 1/b{” (U,, n®) nach den Ergebnissen von Wojaczek 
und Rademacher durch eine passende Funktion approximieren und die 
gefundene Abhängigkeit verallgemeinern. Eine hinreichend anpassungs- 
ist 


Nimmt man den Verlauf von b,” als bekannt an, so wird damit u 


1/b, = = t;. 


und zur Bestimmung von « und f sind nur noch zwei Gleichungen mit zwei 
MeBwerten fiir 6, im Gasgemisch erforderlich. Diese können nun aus den 
vorliegenden Intensitäts- und Elektronentemperaturmessungen gewonnen 


werden. 
Aus der Gl. (1.9) folgt Es 


b,=—- — 2 — 46-10%, 
én D- u? 
6x 


Mit (3.3) und (3,4) erhält man 


(2) 


I, Isa, U, und U, sind den Messungen zu entnehmen. Es darf natürlich 
nicht gleichzeitig 7, und U; = 0 sein. In diesem Falle verliert die Diffusions- 
gleichung ihre Gültigkeit. 


Die gewählten Meßwerte sind in der folgenden Tabelle angeführt: Phe 
x U, U, y n® (x) b, 
x, | 0,892 1,35 - 10-? 8 90 0,893 3,6 - 10 0,88 - 108 
60 | 1,292 0 1 12,3 0,773 1,47 - 1013 3,88 - 10% 


U, und U; sind den Messungen Abb. 14 entnommen. Die Ubereinstimmung 
der Elektronentemperaturkurven aus Sondenmessungen und aus Argon- 
intensitätsmessungen erhöht hier die Sicherheit für die vier angeführten 
Werte. J und J’ entnimmt man der 5-A-Kurve der Abb. 15. n(2 (x) für x = 60 
ist richtig, da in (3.3) die Ankniipfungskonstante 6, „nicht von der Elektronen- 


beweglichkeit abhängt. 
Ann. Physik. 7. Folge, Bd. 3 
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a Werden diese Werte in (IV.4) eingesetzt, so ist damit fiir einen gewissen be- 
Anwendungsbereich (0 < n® < 1015 em-%, 10000° << U,<16000°K, 
nl) — 3,3. 101%) die Elektronenbeweglichkeit im Gemisch Argon— Queck- 
silber als Funktion der Quecksilberkonzentration und der Elektronentem- 
peratur darstellbar. Das Ergebnis ist in Abb. 16 wiedergegeben. Fa 


a: T.=10000 °K 
6:T.- 
d: T.=15000 °K 
e:7.=16000°K 


Abb. 16. Elektronenbeweglichkeit als 
4 and der Quecksilberkonzentration bei verschiedenen 
Elektronentemperaturen. Die Angaben beziehen 
ee 7 sich auf eine Neutralgastemperatur von 70° C 
. und einen Argonfülldruck von 1 mm bei Zimmer- 
temperatur 


Es zeigt sich, daß eine Konzentration der Zusatzkomponente < 1011 cm? 
noch keine Änderung der Beweglichkeit bewirkt. Für n@ > 1012 cm? fällt 
b, dann stark ab, und zwar um so steiler, je niedriger die Elektronentemperatur 
j ist. Zwischen n®) — 6. 1013 und 5 - 10 em”? han: t 6, nur schwach von U, 
ab. Der Einfluß der Quecksilberkonzentration überwiegt hier den der Elek- 
tronentemperatur. Aus dem gleichen. Grund kehrt sich die Abhängigkeit von 
_ U, um, d.h., daß der höheren Elektronentemperatur auch die höhere Beweg- 
lichkeit im Gemisch entspricht. 


An den Stellen 7,85 - 10% und 6 - 101% cm? liest man etwa die in Ab- 
schnitt 8 angeführten, im wesentlichen durch Probieren gefundenen Durch- 
schnittswerte b,, und b, ab, die den Berechnungen des Abschnittes 8 zu- 
grunde gelegt sind. Desgleichen bestätigen einige weitere, mit Hilfe von 
(IV.5) gewonnene Beweglichkeitswerte im Rahmen der Meßgenauigkeit die 
Abhängigkeiten der Abb. 16. Eine endgültige Bestätigung muß jedoch wei- 
teren Messungen, eventuell auf anderer Grundlage überlassen bleiben. 


© 


| 
I 
| \ 
= 
3 8 \ J 
AS 


sen be- 
900° K, 
Queck- 
entem- 


1 em-3 
3 fällt 
eratur 
ron U, 
- Elek- 
it von 
Beweg- 


in Ab- 
Durch- 
8 zu- 
fe von 
eit die 
'h wei- 


H. Morgenroth: Über die Elektrophorese in Argon—Quecksilber-Niederdruckentladungen 395 


Ich danke Herrn Prof. Dr. R. Rompe, Direktor des Instituts für technische 
Physik der Akademie der Wissenschaften zu Berlin, für sein Interesse 
an der Arbeit. Für anregende und fördernde Diskussionen bin ich Herrn 
Dr. K. Rademacher und Herrn K. Wojaczek zu Dank verpflichtet. 


Berlin, Physikalisch-Technisches Institut der Deutschen Akademie der 


f 
Bei der Redaktion eingegangen am 16. rw 1988. ar 


4 


r 
| 
> 


» 


j Oberflächenleitfähigkeit unter Berücksichtigung 
u = der Streuung der Ladungsträger an. der Oberfläche 


Erz 


bei kleinen Bandverbiegungen 


Mit 15 Abbildungen 


Inhaltsiibersicht 


Fir Halbleiter wird die Leitfähigkeit von Raumladungsrandschichten 
parallel zur Oberfläche unter Berücksichtigung der Streuung der Ladungs- 
träger an der Oberfläche berechnet. Die Berechnungen wurden für kleine 
Bandverbiegungen, wie sie im Feldeffekt auftreten, durchgeführt. Die hierbei 
notwendige Integration der Poisson-Gleichung lieferte ein Diagramm, aus 
dem sich in übersichtlicher Weise der Potentialverlauf in Raumladungsrand- 
schichten bei beliebiger Größe der Bandverbiegung und beliebiger Lage des 
Fermi-Niveaus im Halbleiterinneren ergibt. Numerische Auswertungen 
wurden für Germanium durchgeführt. 


at Einleitung 


An der Obfl.!) eines HL!) entsteht durch verschiedene Obfl.-Einflüsse 
eine dünne Schicht, in der die Bänder verbogen sind. In dieser Raumladungs- 
randschicht ist die Leitfähigkeit für Stromfluß parallel zur Obfl. gegenüber 
der Leitfähigkeit im HL-Inneren verändert. Für verschiedene physikalische 
Probleme ist die Kenntnis der Leitwertänderung dieser Schicht in Abhängig- 


En keit von der Größe der Bandverbiegung oder der auf der Obfl. induzierten 
e 3 _ Ladung notwendig. Die Leitwertänderung der Randschicht, die sich unter 
e“ einem quadratischen Obfl.-Stück der Abmessung 1 x 1 cm befindet, gegenüber 


dem Fall unverbogener Bänder für Stromfluß parallel zu einer Kante dieses 
Quadrates, nennt man Obfl.-Leitfahigkeit. Kann man in der Randschicht 
für die Ladungsträger die gleiche Beweglichkeit voraussetzen wie im Inneren, 
an wird die Obfl.-Leitfähigkeit gegeben durch 


o, =e (u, ON + mw, OP), (1) 

wobei ON und 6P die zusätzlich zum Fall der neutralen Obfl. vorhandenen 
Elektronen bzw. Löcher sind. 

Durch Streuung der Ladungsträger an der Obfl. kann aber die Leitfähig- 

keit einer dünnen Obfl.-Schicht verändert werden, nämlich dann, wenn die 

Anzahl der Stöße an der Obfl. vergleichbar wird mit der Anzahl der Stöße 


pS 1) Die Termini ,,Oberfliche“ und ,,Halbleiter‘‘ und ihre Ableitungen sind im folgenden 
7 ihrer Häufigkeit wegen stets mit ,,Obfl.“‘ bzw. „HL‘ abgekürzt. 
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untereinander. Ein solcher Fall tritt z. B. im Feldeffekt auf, wobei induzierte 
Ladung parallel zur Obfl. in einer solch dünnen Schicht fließt. In der Folge 
soll der Fall betrachtet werden, daß die Ladungsträger eine statistische Streu- 
ung an der Obfl. erfahren. Es ist dies der entgegengesetzte Grenzfall zur 
Reflexion, für den keine Beweglichkeitsänderung eintritt. Das tatsächliche 
Verhalten wird zwischen diesen beiden Möglichkeiten liegen. Man kann ver- 
suchen, dieser Tatsache dadurch Rechnung zu tragen, daß man einen Bruch- 
teil als statistisch gestreut, den Rest als reflektiert annimmt?). In der Folge 
soll jedoch nur die reine statistische Streuung betrachtet werden. Ein Energie- 
austausch beim Stoß kann hier, genau wie bei der Wechselwirkung mit dem 
Gitter, unberücksichtigt bleiben. Der Einfluß der Stöße soll durch eine 
effektive Beweglichkeit beschrieben werden, die durch 


0, = € (ünett ON + Uperr OP) (2) 


definiert wird. Effektive Beweglichkeiten in Randschichten unter der An- 
nahme statistischer Streuung wurden zuerst von Schrieffer berechnet?). 
Inzwischen konnte die Annahme der statistischen Streuung an der Obfl. 
auch experimentell gestützt werden**). Die von Schrieffer berechneten Werte 
gelten jedoch nur, falls die Gesamtbandverbiegung groß gegen k T ist. Im 
Feldeffekt-Experiment jedoch liegen die Bandverbiegungen zum großen 
Teil in Bereichen, die außerhalb der Gültigkeit der Schriefferschen Berech- 
nungen liegen. Für dieses Gebiet der kleinen Bandverbiegungen sollen die 
effektiven Beweglichkeiten berechnet werden. 

Für diese Berechnungen war eine Integration der Poisson-Gleichung 
notwendig, um den Gesamtverlauf des Potentials in der Randschicht bei be- 
liebiger Größe der Bandverbiegung und 


beliebiger Lage des Fermi-Niveaus T% ] | 
im HL-Inneren in die Hand zu be- 100 | m. = 


/ 
kommen. Auch die Gesamtzahl der La- | obernaib™ | / 
dungsträger und die Gesamtladung in [| gemischte leitung abernalo72 / 
der Raumladungsschicht wurde be- / 


rechnet, um alle notwendigen Daten 
aus einheitlicher Quelle zu beziehen. 
Diese letztgenannten Daten stimmen 
mit denen aus Kingston und Neu- 
stadters Tabellen*) entnehmbaren [0 Zimmertemperatu_ | 
Werten gut überein. 


|] 


1. Bemerkungen zu den Volumendaten | 


Mp = INg-Nal 


Die zur Auswertung notwendigen 
theoretischen Berechnungen wurden 


10” 10" 10” 0? we 70" 8 09 


bisher nur fiir Erschöpfungs-HL durch- INd-Nal 
geführt. Dies ist aber für Germanium app. 1. Überblick über den Leitungs- 
keine große Einschränkung. Abb. 1 charakter von Germanium 


2) A.H. Wilson, Theory of Metals, Cambridge 1953, S. 242ff. 

3) J.R. Schrieffer, Physic. Rev. 97, 641 (1955). 

3a) J. N. Zemel u. R.L. Petritz, J. Phys. Chem. Solids 8, 102 (1959); R. Miss- 
man u. P. Handler, J. Phys. Chem. Solids 8, 109 (1959). 


> 
4 
3 
oe 
= 
a} 
/ 
= 
| = | 
') R. H. Kingst J. Appl. Phys. 26, 718 (1955) om 


Annalen der Physik. 7. Folge. Band 3. 1959 — 


gibt einen Überblick darüber, welcher Spielraum in Dotierung und Tempe- 
ratur durch diese Berechnungen erfaßt wird. 
Unter Zugrundelegung der Neutralitätsbedingung im HL-Inneren 
1 1 

1+e ** 1+e*? 
‚urde der Temperaturbereich zwischen 7, und 7, berechnet, in dem die 
Elektronen- bzw. Löcherkonzentration auf 5% genau gleich | NV, — N,| ist. 
Oberhalb 7, ist dann gemischte Leitung. 

Fir 7 >T, gilt für Elektronen und Löcher die Boltzmann-Verteilung, 
d.h. 


_I-Ei 

p=n,e (4) 

Für die Berechnungen wurde ein Bandabstand von 0,72 eV und eine Stör- 

stellenaktivierungsenergie von 0,01 eV zugrunde gelegt, die für Störelemente 
| der III. und V. Gruppe des periodischen 


Systems zutrifft5). Es ergibt sich fiir Ger- 
1% manium die Dichte der Ladungsträger im 

eigenleitenden Zustand bei Zimmertempe- 
10% ratur zu n, = 2,4-10%cm®). Der Zu- 
sammenhang zwischen Störkonzentration 
und Lage des Fermi-Niveaus ergibt sich 
im Erschöpfungs-HL zu 


10% Na — Na = Ny — Py= 20, sinh 9, (5) 
wobei g der Wert von 1/kT(£—E,) im 
HL-Inneren ist. 

Da vom Ausgangsmaterial die Leitfähig- 

0” keit gemessen wird, und für die theore- 

tischen Berechnungen die Lage des Fermi- 
Niveaus im neutralen HL-Inneren 9 be- 
kannt sein muß, ist die Abhängigkeit in 

Abb. 2 angegeben. Zugrunde gelegt wurden 

für dieses Diagramm die Arbeiten von 

Prince und van Roosbroeck’)®). Es 

wurde Ausgangsmaterial ohne Kompen- 

sation vorausgesetzt (N, > WN, bei n- und 

N,>N, bei p-Material). 


03 
GF 
Abb. 2. Zusammenhang zwischen 
Leitfähigkeit, Störstellenkonzentra- 
tion und Lage des Fermi-Niveaus 
für Germanium ohne Störstellen- 
kompensation 


2. Der Potentialverlauf in der Randschicht 


Für die Berechnung der Ladungsträgerzahlen und effektiven Beweglich- 
keiten in der Randschicht ist es notwendig, den Potentialverlauf in der Rand- 
schicht zu kennen. Ausgangspunkt hierfür ist die Poisson-Gleichung 


div? = 0 (6) 


5) J. A. Burton, Proc. Intern. Conf. on Semiconductors Amsterdam, Physica 20, 
801 (1954). 
8%) W. van Roosbroeck, Physic. Rev. 91, 287 (1953). 
oar 7) N. B. Prince, Physic. Rev. 92, 681 (1953); 98, 1204 (1954). 
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mit 
für vollständig dissoziierte Störstellen. 


Im neutralen HL-Inneren ist «= 0. Für u > 0 nimmt die Anzahl der 
Löcher zu und die der Elektronen ab, also ist o > 0; fürau< Oist es umgekehrt. 


wird die Poisson-Gleichung im eindimensionalen Fall 
- (Ne - N, — My + Py ®*). (9) 


Die erste Integration wird mit der Bedingung u’ = 0 fiir u = 0 durchgeführt 
und ergibt 


w= (Ne —N,)u—ny (l—e™) + py (e*—1). (10) 


Für u > 0 gilt das negative Vorzeichen, für u< 0 das positive Vorzeichen 
Der Ausdruck unter der Wurzel ist für beliebige u stets > 0. 

Die weitere Integration ist geschlossen nicht durchführbar. Zur näherungs- 
weisen Integration werden für den Radikanden in Gl. (10) folgende > 


ausdrücke verwendet: u ve 
a) (my + py) fiir kleine Werte von |u| 

b) (N,—N,) u im Erschöpfungsgebiet 


c) bzw. p,e* in Elektronen-Anreicherungsschichten auf n-Ma- 
terial bzw. in Löcheranreicherungsschichten auf 
p-Material und in Inversionsschichten. 


Die Berechnungen haben gezeigt, daß man am besten die Näherung a) 
für |w| < 2 benutzt und daran je nach dem Wert von N,— N, entweder b) 
oder c) anschließt. Der Übergang von der Erschöpfungsschicht zur Inver- 
sionsschicht wird für den Wert u; vorgenommen, für den „Näherung b) = 
Näherung ist. 


a) Die Näherungslösung für den Fall |u| < 2. % 
Gl. (10) geht in der Näherung a) über in “ae 

mit he 
kT [ek T 
= 1 — 12 
le Ny + Pa 2n, 
Als Lösung hiervon ergibt sich einheitlich für positive und negative u 
u= u, e~*lh, (13) 


u, ist der Wert von wu an der Obfl. Die exakte Lösung im Bereich |u| < 2 
kann man sich durch Gl. (13) mit ortsabhängigem J, gegeben denken. Ein Ver- 
gleich von Gl. (11) mit Gl. (10) zeigt, daß dieses exakte /, für |u| = 2 um 
an: 30% von dem in Gl. (12) gegebenen Wert abweicht. 
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b) Die Näherungslösung für |u| > 2im Erschöpfungsgebiet. 

Ein Erschöpfungsgebiet tritt für n-HL bei Elektronen-Verarmung, für 
p-HL bei Löcher-Verarmung dann auf, wenn die Raumladung nur durch die 
Störkonzentration bestimmt ist. Die Bedingung hierfür ist 


N,—N,| >n,ele#+2, (14) 


Wegen Gl. (5) erhält man für |p| die Forderung |p| > 1,1. Gl. (10) geht in 
in dieser Näherung über in 


2(Na— - N,) 9 n, sinh @ 
u + 2e// u, tfirus 0. (15) 


Die Lösung für positive und negative u zusammengefaßt ist 


IN. — N,| n, sinh || 
Die Schichtdicke, in der das Potential von seinem Wert an der Obfl. auf den 
Wert + 2 abfällt, ist 


d, =I) 2ctgh|p| (Viuj -V2). 
Da |p| > 1,1, ist etgh |p| ~ 1 und es wird van 


(8) 


Gegenüber dem Fall a) ergibt sich hier eine merkliche Änderung der Rand- 
schichtdicke mit der Größe der Bandverbiegung. 
Die Lösungen (13) und (16) gehen mit stetiger Tangente ineinander über. 
Der Übergang zur Inversionsschicht wird für den Wert von |w,| vorgenommen, 
für den 

[Ne mer 


ist, oder unter Zuhilfenahme von Gl. (5) el, 


c) Näherungslösung für | «| >2 bei Anreicherungs- und Inversionsschichten. 
Die folgende Näherung findet Anwendung für Elektronen-Anreicherung 
auf n-Material und Löcher-Anreicherung auf p-Material sowie für Inversions- 
schichten; also für alle in b) nicht behandelten Fälle. Bestimmen die Löcher 
die Raumladung, dann wird aus Gl. (10) 


ekT 

Bestimmen die Elektronen die Raumladung, dann ist 
=e 2ny -u/2 — ] -(u-g)/2 
== 2 = u< 0. (22) 

Die Lösung kann für beide Fälle wieder zusammengefaßt werden zu 

erlull2 — e-lusl/2 — etPl2.2/2 L + fiir 0. (23) 


Dabei gilt das Pluszeichen fiir Anreicherungs- und Inversionsschichten der 
Elektronen, das Minuszeichen für solche der Löcher. 


elwil = 2 (e2lel — 1). (20), 
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Schließt die Anreicherungs- bzw. Inversionsschicht direkt an den Wert 
|u| = 2 an, dann ist die Dicke dieser Schicht nach Gl. (23) unter Berück- 
sichtigung von Gl. (12) 

d, = 21, V cosh — er \wll2) + fürru (24) 


Beginnt die Inversionsschicht erst für den Wert u, (Anschluß an Erschöp- 
fungsschicht), dann ist ihre Dicke 


d, = cosh p e*#/2 (e-!mil2 — e-lusl/2) fir uS 0. (25) 


Aus den in a), b) und c) gegebenen Lösungen läßt sich der gesamte Verlauf 
des Randschichtpotentials ermitteln. Abb. 3 gibt eine solche Schar von 


u 


ata 10 05 1 


LL 


0 05 10 15 20 25 30 4 


Abb. 3. Der Potentialverlauf in der Randschicht für _ 
verschiedene Werte des Parameters gy. Die Kurve fiir | 
y = — 2 gilt gleichzeitig für alle p < — 2 


Potentialverläufen für u > 0 (Parameter 9). Der Nullpunkt der Abszisse ist 
so gewählt, daß man für 2/l, > 1 direkt die Randschichtdicke entnehmen kann: 


l,=1,+ 4, bei Erschépfungsschichten 
1=4+4, bei Anreicherungsschichten (26) 


l, =1,+d,-+-d; bei Inversionsschichten. 


Der Zusammenhang zwischen der Ortskoordinate z und der im Diagramm ver- 
wendeten Koordinate = ist 


(27) 


Es ist zu beachten, daß /, und damit der Maßstab der Abszisse nach GI. (12) 
von p abhängt. Für «< 0 erhält man die entsprechende Kurvenschar unter 
Berücksichtigung der aus Symmetriegriinden unmittelbar einleuchtenden 
Beziehung z (— u, 9) = z (u, —¢). 


3. Bestimmung der Ladungsträgerzahl in der Randschicht 


Nunmehr läßt sich die Anzahl der Ladungsträger in der Randschicht in 
Abhängigkeit von der Bandverbiegung berechnen. Es ist hier zweckmäßig, 
die Abweichung von der Anzahl bei ungeladener Obfl. pro Flächeneinheit der 


+) 


= 


% 
für 3 
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Obfl. anzugeben. Die Anzahl der zusätzlichen Löcher ist also tlie j 


= 6P = Pu) dz = Py 1) dz = | (e*—1)dz (28) 


und die Anzahl der zusätzlichen Elektronen 


= de = my ef (e-"—1)dz. (29) 


| Zwischen dN und öP besteht die Beziehung 
B ON (u, 9) = 6P (—u, —p). (30) 
f Es genügt also, dP zu berechnen, und 6N dann mit Hilfe von Gl. (30) zu be- 
stimmen. 
27 Für |u| < 2 wird u aus Gl. (13) entnommen und es wird mit Berück- 
As sichtigung von Gl. 
31) 


Für wu > 2 sind zwei Fälle, p< 1,1 und g >1,1, zu unterscheiden. Im 
Fall 9< 1,1 entnimmt man u der Gl. (23) und es ergibt sich unter Be- 
rücksichtigung der Gl. (24) und (12) 


| 
. 


(32) 
=n,L4e 1) u,>2,9< 1,1 


dP, ist der Wert von 6P für u = 2. 

Fir gy >1,1 wird das Potential durch die Erschöpfungsschicht der 
Elektronen bestimmt und wu der Gl. (16) entnommen. Damit wird unter Be- 
rücksichtigung von Gl. (18) und (12) 


= 


d, 
OP — py S (e“ _ 1) dz, > Us > 2,9 >11. 
e"dy+Y 2 Vs) 


_ Das hier auftretende Integral wurde in Abhängigkeit von der oberen Inte- 
grationsgrenze numerisch integriert. Das Ergebnis ist in Abb. 4 angegeben. 


Be Für u, > 4 ist gut brauchbar 


E Für den Wert u, , geht die Erschöpfungsschicht in eine Inversionsschicht über 
und wird aus Gl. entnommen. Unter von Gl. (25) 
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und (12) ergibt 


L4ee2 (cosh — cosh u >u,p>Li. 
(9) 


7 
103 AA 
(30) 
u be- = 
x 
| 
0° 


4. Ergebnis der numerischen 
Ve 
Integration von f e” dt 
. Im Ya 
r Be- 


Für u< — 2 sind wieder die beiden Fälle — 1,1 und g > — 1,1 zu unter- 
scheiden. Im Fallp > — 1,1 entnimmt man u der Gl. (23). Mit Berücksichtigung 


der Gl. (24) und (12) ergibt sich 
da 
6P — 6P_, = Py | [GEL dz 
0 


— e- 1/3 ( e-3 — e-3lwll2)] 
u, — 23,9 — 11. 


von Gl. (18) und (12) 


e 


— 6P_, = Dy f (e _ 1) dz 
A 0 
3 Vius| ( 
Inte- | Für den Wert u, ; geht die Erschöpfungsschicht in eine über, 


eben. [|% wird aus Gl. (23) entnommen. Es ergibt sich — Berücksichtigung von 
Gl. (25) und (12). 


di 
cgay | OP — OP = Pu dz 
0 (38) 


=—n,L2 3¢/2 [e «1/2 — e-Iusli2 — 1/3 (e 3 lusl/2 e~ 3 |us|/2)] 


U< Uyi,9 >—1,1. 
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Fir g< — 1,1 wird das Potential durch die Erschopfungsschicht der Löcher 
(16) entnommen. Damit wird unter Berücksichtigung 
q 
it über 
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Die bei den Berechnungen stets wieder auftretenden Funktionen en 


t #2 3 
(t) = cosh — cosh 1 


/ 


Vt 

Br sind in Abb. 5 angegeben. 
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Abb. 5. Die Funktionen g, (é), ..., 9; (t) nach Gl. (39) 
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Abb. 6. Gesamtzahl der Löcher in der Randschicht in Abhängigkeit von der Band- 
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öP in Abhängigkeit von u, für gebräuchliche Werte des Parameters ¢ ist 
in Abb. 6 für Germanium angegeben. Den Berechnungen wurde zugrunde 
gelegt 


cm? c m? 


34-10 cm „= = 38007 = = 18207 (40) 


Sind 6P und 6N bekannt, dann läßt sich auch sofort die in der Randschicht 
befindliche Ladung Q, bestimmen. Es ist 


(41) 


Ebenso läßt sich die Obfl.-Leitfähigkeit ohne Berücksiehtigung der Streuung 

der Ladungstriger an der Obfl. nach Gl. (1) nunmehr berechnen. 

4. Berechnung der Verteilungsfunktion 

Zur Berechnung der effektiven Beweglichkeit benötigt man weiter die 

Verteilungsfunktion. Ausgangspunkt für die Berechnung der Verteilungs- 
funktion f ist die 


a, 4 


z ist wieder die Koordinate senkrecht zur Obfl., z = 0 die Gleichung der 
Obfl. In x-Richtung fließt der Strom. qg ist die Ladung der Elektronen bzw. B 
Löcher. t ist die Relaxationszeit, die im allgemeinen noch von der Geschwin- Be; 
digkeit abhängt. Im folgenden soll jedoch t = const angenommen werden, i 
die Feldabhängigkeit der Beweglichkeit also unberücksichtigt bleiben. 
fo ist die Verteilung im Gleichgewicht, also für E, = 0: 


+... Elektronen 
o= Ce _ für Löcher (43) 
C bestimmt sich aus der Dichte der Ladungsträger im neutrale on H HL zu ; Pr = 
m 3/2 + Elektronen 
et? ii 
= k 7) für scher. sida 


Fön: 
Die folgenden Rechnungen sollen für Löcher durchgeführt RA Die Er- 
gebnisse für Elektronen erhält man dann auf Grund der bereits erwähnten 
Symmetrie. 

Mit dem Ansatz f = f, + f, und unter Berücksichtigung der Koordinaten- 
transformation 


—e E,x(E,= const) ,= sv2+kT(u—u,) (4) 
geht (42) über in 


é Bah 1 e ; 


Die statistische Streuung an der Obfl. fordert, daß f, = 0 wird an der ts 
Obfl.; die Verteilung geht in die reine statistische Verteilung f, über. Da man a 
sich im vorliegenden Falle auf eine Lösung beschränken kann, die linear in 
E,ist, kann in der (46) das Produkt E, v, Ver- 
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nachlässigt werden. Die so vereinfachte Differentialgleichung lautet 
pth she 


Für den vorliegenden Fall genügt es weiter, diese Differentialgleichung unter 
der Bedingung E, = const zu integrieren. Damit wird eine Geschwindigkeits- 
 verteilung in einer idealisierten Randschicht (Abb. 7) betrachtet. Die Lösung 
unter der vorgegebenen Bedingung an der Ar ist Lit 


ux 406 


hi ~ 
= 


Der Unterschied in den Verteilungsfunktionen für zu > bzw. <eE,z 


hat folgende Ursache. Die Ladungsträger, für die Zu <eE,z ist, können 


u a gr wegen zu kleiner kinetischer Energie nicht zur 
ba Obfl. gelangen und darum nicht am Streuprozeß 
teilnehmen. Ihre Gleichgewichtsverteilung wird 
‘ nur durch das elektrische Feld E, gestört. Der 
rd he Exponentialterm in der oberen Gl. (48) gibt 
rae he also den Einfluß der Streuung an der Obfl. 


wieder. 


Im Grenzfall E, gegen Null ergibt sich 
ATus_ - 
ef, a ( ) 
Abb. 7. Idealisierter Poten- O14 ziel 
tialverlauf fi fo\l—e m). (49) 


KERN Hierin bedeutet f§ den Wert von f, für «= 0. Durch f, wird wegen 
DS f f, dv, =0 keine Konzentrationsinderung der Ladungsträger bewirkt 


es gilt darum weiterhin 


1 


5. Einfiihrung der effektiven Beweglichkeit 
Die Stromdichte i im Abstand z von der Obfl. ist ey hie = 
5 
tne (z) =e So }, dv, dv, dv, =. 
 =eu,P(e) E, (1- f dt 
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benutzt. Der erste Term auf der rechten Seite von Gl. (51) gibt die ungestörte 
Leitfähigkeit, während der Exponentialterm den Einfluß der Stöße mit der 
Obfl., eine Verminderung der ungestörten Leitfähigkeit, angibt. In der Nähe 
der Obfl. ist also die Beweglichkeit ortsabhängig. Unter Berücksichtigung der 
Beziehung 

ing (2) = em (2) p(z) By (53) 
und Gl. (51) erhält man sein 


-00 


Für z = 0 ist die Beweglichkeit u (0) = 0. Sie wächst mitgrößer werdendem Ab- 
stand von der Obfl. an, und es wird für große z: u (2) = u 
Für den vorliegenden Fall ist es zweckmäßiger, eine effektive Beweglichkeit 

zu definieren, die auf den gesamten zusätzlichen Strom in der Randschicht 
bezogen ist und definiert wird durch nu 

=f (i (2) — (00)) dz = pert OP E,. (55) 

Schrieffer berechnet effektive Beweglichkeiten, wie sie durch Gl. (55) 
definiert werden, für zwei Fälle. Einmal für den Fall der Erschöpfungsschicht, 
deren Potential wie in Abb. 7 gezeigt idealisiert wird, und weiter für größere 
Bandverbiegungen, bei denen noch ein Inversionsgebiet auftritt. Reine An- 
reicherungsschichten werden nicht betrachtet. Für das Gebiet kleiner Band- 
verbiegungen können Schrieffers Rechnungen nicht ohne weiteres über- 
nommen werden. Der hierfür erforderliche Grenzübergang u, > 0, « > 0 ist 
nicht durchführbar, da die Lösung ein singuläres Verhalten zeigt und das 
Modell dann physikalisch auch nicht mehr zulässig ist. 
Es soll darum hier ein Weg beschritten werden, der unter Verwendung der 
Verteilungsfunktion Gl. (49) zu Werten der effektiven Beweglichkeit bei 
kleinen Bandverbiegungen führt. Zunächst wird ein einfaches Modell betachtet. 


6. Vereinfachtes Modell zur Berechnung der effektiven Beweglichkeit in einer 
dünnen Schicht 


Es sei in einer Obfl.-Schicht 0 < z</ die D:chtu der zusätzlichen Löcher 
konstant, z. B. = pp. Der zugehörige Strom ist dann 


I» fidv,dv,dvu,dz. 166) 
Hierin ist ff aus Gl. (49) mit 
m \-3/2 
(5; k 7) _ 
einzusetzen. Gl. (49) und (67) in (56) epee ergibt 


ve = er) f vz e "ar “(4 du, dv, dv, dz. (58) 
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Dieser Strom wird aufgeteilt in 

pp k 7) J wer du, de, dv,dz 65 
0 
und 

m -3/2 3; - 
I,.= Pol k r) J dy, dv, dv, dz. (60) 


I, ist der Strom ohne Berücksichtigung der Stöße an der Obfl. Unter 
Br horton der Gl. (52) und der Tatsache, daß p,! = 6P ist, wird 


In = ep E,: (61) 
Aus I,, wird mit Hilfe von 3 37 
[ze dx =5-|/— (62) 
und nach Ausführung der Integration über z 
° 7 2yn oo ( ) 
_ Fiithrt man hier die freie Weglänge 4 = 3/2 | dk t ein, dann ergibt sich 


i 


v? als neuer Integrationsveränderlicher. Der Gesamtstrom ist 


ET 


= 152 — (65) 

Nach der Definition (55) ergibt sich unter Berücksichtigung von Gl. (61) und 
(64) eine effektive Beweglichkeit in dieser Schicht zu cf rae 
3Ya 1 1 3 

dt . (66) 


oo ı 
4 Aye um 
fe dt~1 a7 (67) 


und damit sowohl Mer als auch duer/dl = 0 für 1/A = 0. Für 1 > co geht 
das Integral gegen Null und damit zer gegen u. Dies sind Ergebnisse, die man 
auch physikalisch erwartet. Mit der neuen Veränderlichen 2 = a - + 


wird Gl. (66) zu 
té 


Mett =—]— f e-t- z/Vt dt q (68) 
0 


Da diese vereinfachte Obfl.-Schicht als Modell fiir diinne Schichten ihre Be- 
ie de »utung hat und aus den hier gewonnenen Daten die Näherungsergebnisse 
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für die realen Randschichten abgeleitet werden, soll Gl. (68) eeannuies aus- 
gewertet werden. 
Die numerische nun. des Integrals - 


als Funktion ‘en FREE x wurde für Parameterwerte im Bereich 
0,025 < x < 10 auf drei Stellen 


genau durchgeführt. Für «>10 Fix) (delt 
ist praktisch F(x2)=(0, und "od 
damit 07 
06 
Va x 03 
02 7 
Für «< 0,025 ist praktisch a 4 


Hert = 90. Das Ergebnis dieser Om? 512 502 5 100 


Integration ist in Abb. 8 ange- Abb. 8. Ergebnis der numerischen Integration 
geben. Die Kreuze geben die be- m. 
= 


rechneten Werte an. Die Berech- f ur 
nung der effektiven Beweglichkeit (x)= e ‘dt und die effektive 


nach Gl. (68) ist ebenfalls in Abb. 8 Beweglichkeit des vereinfachten Modells für 
dargestellt. dünne Schichten nach GI. (68) 


7. Anwendung des einfachen Modells auf reale Randschichten Er ae 8 


Dieses einfache Modell wird nun in folgender Weise für die Berechnung der 
effektiven Beweglichkeiten in HL-Randschichten für kleine Bandverbiegungen 
verwendet. 

Auch für die allgemeine Randschicht wird die Aufteilung (65) des Stromes 
vorgenommen. Der Anteil /,, ohne Berücksichtigung der Streuung ist wieder 
Ing =euöPE,. Für MR läßt sich eine Abschätzung angeben. Diese Ab- 
schätzung muß für wu, > 0 und u,< 0 getrennt durchgeführt werden. Zu- 
nächst wird der Fall «, > 0 betrachtet, für den an der Obfl. die größte Konzen- 
tration herrscht. Man wird sicher / px ZU grow ansetzen, wenn es nach Gl. (64) 
berechnet wird mit / als Randschichtdicke und py = p, — py, also der größten 
Überschußkonzentration. Da die Störung durch Stöße mit der Obfl. im Ab- 
stand A praktisch abgeklungen ist, wird FAR sicher zu klein angegeben durch 
Gl. (64), wenn man für p, den Wert p, — p,, wählt, wobei p, der Konzentra- 
tionswert im Abstand A von der Obfl. ist. Der ne Wert w — also 
zwischen diesen beiden Grenzen liegen: BR, 2 


Der obere Wert würde eine zu kleine, der untere Wert eine zu große effektive 
Beweglichkeit liefern. Wird für Ei , das arithmetische Mittel aus oberem und 
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unterem Ww ert verwendet, dann wird die effektive Beweglichkeit 


® 
= 
= 


| 


Hierin ist p, = 1/, (p, + pa). Spezialisiert man diese Gleichung auf das 
von Schrieffer verwendete Modell, dann sieht man, daß Gl. (72) für größere 
Bandverbiegungen in befriedigender Weise in die dort gültige Schrieffer- 
Formel einmündet. Der Übergang zur Schrieffer-Formel findet dann statt, 
wenn der Abfall des Potentials am Rande auf einer freien Weglänge größer 
als k T wird. 

An einem Modell, für das die effektive Beweglichkeit exakt berechnet 
werden kann, soll die Brauchbarkeit des Vorgehens gezeigt und Wege zur 
eventuellen Erweiterung des Gültigkeitsbereichs der Formeln aufgezeigt 
werden. 

Gl. (49) legt nahe, in grober Näherung ein Verhalten der Beweglichkeit wie 


anzusetzen, wobei e~*/* den Einfluß der Streuung angibt. — 
Mit der weiteren Annahme 


P — Py = Poe 


— dies entspricht dem linearen Potential der Abb. 7 bei sehr ont u, und 
geeigneter Normierung — wird 


Die effektive Beweglichkeit wird damit 


Wie zu erwarten, wird 


Met = 0 für 2 = 0 und 
Her = u für 2, = oo. 
Berechnet man für 
dieses Modell die effek- 
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7, 
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Abb. 9. Vergleich der Näherungen für ein Modell, bei dem Ten Möglichkeiten ver- 
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hiervon zeigt sich für dieses Modell das arithmetische Mittel aus (77) und der 
unteren Grenze. 

Dieser Vergleich gibt einen Hinweis auf eine Möglichkeit der weiteren Ver- 
besserung der Annäherung der effektiven Beweglichkeit durch das arithme- 
tische Mittel. 

Für u,< 0 muß berücksichtigt werden, daß ein gewisser Teil der Löcher 
durch das Potential daran gehindert wird, an die Obfl. zu gelangen. Es soll 
angenommen werden, daß die Löcher im Mittel das Potential |u, — u;|/2 
zu überwinden haben, und die zum Stoß gelangende Anzahl damit 
(Ps — Py) e~'“s-“al/2 wird. Der Näherungswert für die effektive Beweglich- 
keit wird damit 


f Ps— P a 2 4 
Ime-wall2 A(1—F(z)) 0. (78) 
Für große Werte von u, geht per gegen u, für u, = 0 schließt sich (78) 
stetig an (72) an. Es ist 


im far _ (79) 
2X 


u 
wobei x, der Wert von x für u, = 0 ist. 
Nach den Gl. (72), (78) und (79) werden die effektiven Beweglichkeiten 
für kleine Bandverbiegungen berechnet. Die angeschriebenen Formeln gelten 
ihrer Herleitung nach für Löcher. Genau genommen hätte man die Größen 
u, A und I noch mit dem Index p versehen müssen. Ersetzt man diese Größen 
durch die entsprechenden der Elektronen und weiter 6P durch dN, p, — py 
durch n, — n,,, dann gelten die Gl. (72), (78) und (79) für Elektronen, und zwar 
Gl. (72) für u,< 0, Gl. (78) für u, > 0. 


— Ceétto 
— 
09 tLöcher 
etronen 025m 19 ser 


05 

Abb. 10. zu, 0/12 für Löcher und Elektronen Abb. 11. Die Temperaturabhängigkeit von 
an der Obfl. von Germanium Her o/f¢ für eigenleitendes Germaniuia 


Als Ausgangsdaten für die Berechnung im Falle des Germaniums werden 
wieder die Werte (40) benutzt. Hieraus ergeben sich die freien Weglängen 


für Löcher bzw. Elektronen zu A, = 1,32- 10-5] = em und A, = 


‘em. Hierbei sind m, und m, die effektiven Massen der Löcher 


2,75 - 105 | — 

7 
und Elektronen, die als Skalar betrachtet werden sollen mit dem Wert m, = 
Q7* 


Bar 
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m, — 0,25m in Übereinstimmung mit Schrieffer und anderen Autoren. 
Effekte, die durch den Tensorcharakter der effektiven Masse hervorgerufen 
werden®), bleiben hier unberücksichtigt. 

Einen ersten Überblick über die Größe des Effektes gewinnt man durch 
Betrachtung von werro/u nach Gl. (79) in Abhängigkeit von g und der Tem- 


en peratur. In Abb. 10 ist wero/u 
eo in Abhängigkeit von p bei Zimmer- 
| temperatur aufgetragen, und 
N) zwar für die beiden Fälle m, = 

m, = 0,25 m und m, = m, = m. 

NN Die Werte, die sich für ", = 
08 an ~ m, =m ergeben, sind die für 


u,=O0 kleinstmöglichen Werte 


0? N a | der effektiven Beweglichkeit. Für 
N \ We | sehr reine Materialien ist also der 


Effekt bei Zimmertemperatur nicht 


\ N Q EN sehr groß. 
05 Y=) AA NV Dies ändert sich aber bei höhe- 


2 0 4 2 -6 ~-t#us-@ ren Temperaturen. Um einen 
Abb. 12. Effektive Beweglichkeit der Elek- Überblick hierüber zu geben, soll 
tronen in Abhängigkeit von der Bandverbiegung die Temperaturabhängigkeit der 
an der Obfl. von Germanium effektiven Beweglichkeit von Lö- 
chern und Elektronen für u, = 0 
in eigenleitendem Material angegeben werden. Eigenleitendes Material wird 
deshalb gewählt, um die Temperaturabhängigkeit von p zu eliminieren und 
weiter auch, weil es eine obere Grenze für die effektiven Beweglichkeiten 
angibt. 
Unter Berücksichtigung der Temperaturabhängigkeit von Z nach Gl. (12) 
und derjenigen der freien Weglänge im Gebiet der Streuung an Schallquanten, 
ergibt sich eine Tempe- 


Cet raturabhängigkeit von 2 
zu 
> = x» (T) = 2, (300) 


x, Die mit diesem Wert 

07 N SN Ww auf Grund der Gl. (79) er- 
\\N rechnete Temperaturab- 

06 \ hängigkeit ist in Abb. 11 


05 Y= A A 2 Die effektive Beweglich- 
Abb. 13. Effektive Beweglichkeit der Löcher in Ab- 


_ hangigkeit von der Bandverbiegung an der Obfl. von sich nun ebenfalls bestim- 
hee Germanium men. 
‘ae = 8) J. R. Schrieffer, Semiconductor Surface Physics, Hrsg. R. H. Kingston, 


Fr 


Philadelphia 1957. 
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Abb. 14. Obfl.-Leitfahigkeit unter Berücksichtigung 
der Streuung der Ladungsträger an der Obfl. o, in Ab- — 
hängigkeit von der Bandverbiegung für Germanium 


In den Abb. 12 und 13 sind die nach Gl. (72) und (78) bereehneten Werte 
für Elektronen und Löcher für einige gebräuchliche Werte des Parameters ~ 
angegeben. 

Nach Gl. (2) ergibt sich nun auch 


die Obfl.-Leitfähigkeit unter der An- N \\G use” 
nahme statistischer Streuung der La- Gy \\ 15 
dungsträger an der Obfl.- Diese Obfl.- 

Leitfähigkeit wurde in Abb. 14 darge- \W 49. 

stellt. Abb. 15 gibt fiir eigenleitendes AM 10 
Material einen Vergleich der Obfl.- \\ 
Leitfähigkeitskurven ohne Modifikation \ 5 

der Beweglichkeit an der Obfl. (Kurve N 

ö,), Schrieffers Näherung und die N Vi 
hier berechnete Näherung für kleine 0 lf 
Bandverbiegungen (Kurve o,). 16 2 0 2 4 Gus 


Abb. 15. Vergleich der Obfl.-Leitfahig- 
1 0. Hachenberg keitskurven mit und ohne Beriicksichti- 
danke ich fiir das ördernde Interesse, gung der Streuung und Schrieffers 
das er der Arbeit entgegenbrachte. Lösung für eigenleitendes Germanium 


Berlin-Adlershof, Heinrich-Hertz-Institut der Deutschen Akademie 


der Wissenschaften zu Berlin. 
Kine Bei der Redaktion eingegangen am 22. Januar 1959. Fs) 


47 ° 

= 
= NYY / 
der N 
nicht 
J = 


Feldeffekt und 
an der Oberfläche von Germanium _ 


Von Hans Flietner 


Mit 11 Abbildungen 


Inhaltsübersicht 


Mit Hilfe des Feldeffektes wurden Reaktionen verfolgt, die an der Ober- 
fläche von Germanium-Proben verschiedener Leitfähigkeit stattfanden. Es 
wurde festgestellt, daß die Einwirkung von Ozon eine Vergrößerung der An- 
zahl der ,,schnellen‘‘ Terme mit sich bringt. Im Fall der nicht mit Ozon vor- 
behandelten Oberfläche war dieser Zuwachs an ‚schnellen‘ Termen begleitet 
von einem Verschwinden „langsamer‘‘ Akzeptor-Terme. Diese „langsamen“ 
Akzeptor-Terme, die durch Sauerstoff-Adsorption zu verstehen sind, werden 
durch die Ozonbehandlung in ‚schnelle‘ Akzeptor-Terme übergeführt, die 
Bindungsmöglichkeiten von freien Elektronen an freien Valenzen von Ober- 
_ flächen-Atomen entsprechen. Durch eine feuchte Atmosphäre wird die Wir- 
all kung auf die „schnellen“ Terme wieder vollständig rückgängig gemacht. 
Es bleibt jedoch die Ladung in den „langsamen“ Termen erhalten. Dies 
bedeutet, daß durch die Vermittlung von H,O ein Bindungstyp an der 
Oberfläche entsteht, der keine Ladungsänderung mit sich bringt. Ein Bindungs- 
typ mit diesen Eigenschaften ist die GeO,-Bindung. Bei einer mit Ozon und 
Feuchtigkeit vorbehandelten Oberfläche bewirkt eine erneute Ozonbehand- 
lung ebenfalls eine Vergrößerung der Anzahl der ‚schnellen‘ Terme, jedoch 
ohne gleichzeitige Ladungsänderung in den ‚langsamen‘ Termen. Man ver- 
steht dies so, daß Ozon ebenfalls die GeO,-Bindung auftrennt, um Bindungs- 
möglichkeiten für freie Elektronen an freien Valenzen zu schaffen.  ~ 


Einleitung 


Wird auf der Obfl.!) eines HL!) Ladung induziert, so verbiegen sich die 

re Bänder in einer dünnen Obfl.-Schicht und die Leitfähigkeit ändert sich. 

a coe Diesen Effekt nennt man Feldeffekt. Mit Hilfe dieses Effektes ist es möglich, 

die in Obfl.-Termen befindliche Ladung in Abhängigkeit von der Lage des 

ee Fermi-Niveaus zu bestimmen?) Diese Kenntnis ermöglicht Aussagen über 


De Lage, Dichte und Charakter der Obfl.-Terme. Insbesondere ist es möglich, 


!) Die Termini ,,Oberflache“ und ‚Halbleiter‘ und ihre Ableitungen sind im folgenden 
ihrer wegen stets mit „„Obfl.‘“ und abgekürzt. 
2) W. L. Brown, Physic. Rev. 100, 590 (1955). 
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Das Modell, das heute fiir Obfl.-Untersuchungen herangezogen wird, ist 
in Abb. 1 gezeigt?)®). 
Auf der Obfl. des Kristalls befinden sich zusätzliche Terme, deren Energie- 
niveaus im Gebiet der verbotenen Zone liegen. Diese Terme nehmen dadurch 
Einfluß auf die elektrischen Eigenschaften, daß sie Elektronen mit dem HL 


austauschen. Hierdurch 
entsteht an der Obfl. eine 
Raumladungsschicht. Die neutraler naibleiter 
die Gesamtladung auf der Se it 
Obfl. und die Obfl.-Terme ee 


festgelegt. Diese  Obfl.- 


> 
Terme kann man nach & = 
ihren Reaktionszeiten in t F „langsame Terme 
zwei Klassen einteilen: In 
Fermi-Niveou 
die sogenannten ,,lang- & -*=- 


samen‘ Terme, deren Reak- 
tionszeiten zwischen 0,1 sec 
und einigen Minuten liegen, 
und in die „schnellen“ sanatt (p-Typ) 
Terme mit Reaktionszeiten Abb. 1. Das Termschema de: freien HL-Obfl. 

in der Größenordnung von 

10 usec. Die Anzahl der ‚langsamen‘ Terme ist 1013—10!4 cm?. Über ihre ener- 
getische Lage ist noch wenig bekannt°), Ref. Bardeen. Die Anzahl der 
„schnellen‘‘ Terme ist viel geringer, doch gelang es, genauere Angaben über 
ihre Lage und Dichte zu machen. Es gibt im wesentlichen eine stetige Ver- 
teilung um die Bandmitte herum mit einer Dichte von 101—1012 em”? eV-1 
und weitere Terme in größerer Dichte in einem Abstand von der Bandmitte, 
der größer ist als etwa 0,1 eV), Ref. Garrett. 

Anlaß zu den vorliegenden Untersuchungen war die Feststellung des 
Effektes, daß Ozon eine Veränderung der Obfl.-Leitfähigkeitskurven im Feld- 
effekt bewirkte, die auf eine Änderung des Spektrums der ,,schnellen‘‘ Terme 
schließen ließ. Da sich einerseits im Feldeffekt eine Änderung der Konzen- 
tration und Lage der ,,schnellen‘‘ Terme gut verfolgen läßt, andererseits eine 
chemische Reaktion die Terme verändert, versprachen diese Untersuchungen 
einige Aufschlüsse über die Bindungsverhältnisse. 

Man kennt schon länger Einflüsse, die die ‚schnellen‘ Terme verändern, 
z.B. Schmirgeln und Ätzen. Wang und Wallis®) fanden weiter eine Än- 
derung der „schnellen“ Terme kurz nach dem Ätzen. Erst einige Zeit nach dem 
Ätzen ist also die Termstruktur auf der Obfl. stabil, d.h. zeigt keine vom 
Ätzprozeß mehr herrührende Veränderungen. Die hier beschriebenen Än- 
derungen der ‚schnellen‘ Terme wurden an stabilen Obfl. herbeigeführt. 
Nach der Veränderung der Termstruktur durch Ozon konnte durch bloßen 
Wechsel der umgebenden Atmosphäre die Ausgangssituation wieder hergestellt 


3) R. H. Kingston, Review of Germanium Surface Phenomena, J. Appl. Phys. 27, 
101 (1956). 
4) Semiconductor Surface Physics, Hrsg. R. H. Kingston, Philadelphia 1957. ee ta 
5) Intern. Koll. üb. Halbleiter u. Phosphore, Garmisch-Partenkirchen 1956. ni 
6) S. Wang, G. Wallis, Physic. Rev. 107, 947 (1957). en Kiss a 
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werden, d.h. die an der Obfl. ablaufenden Reaktionen wurden reproduzierbar 
beherrscht. 

Bei der Ermittlung der Lage der Fermi-Niveaus an der Obfl. relativ 
zu den Bändern aus den gemessenen Obfl.-Leitfähigkeitskurven wurde die 
Streuung der Ladungsträger an der Obfl. berücksichtigt”). 


1. Vorbereitung der Proben 


Als Ausgangsmaterial wurden Germanium-Einkristallscheiben bekannter 
Kristallorientierung von etwa 0,5mm Stärke verwendet. Diese Scheiben 
wurden rechteckig auf etwa 6x 12 mm zurechtgeschnitten, mit feinem Schmir- 
gel auf etwa 0,3mm Dicke abgeschmirgelt und anschließend poliert. Als 
letzte Obfl. „Behandlung folgte eine Ätzung mit CP 4 und gründliche Spülung 
in destilliertem Wasser. Eine Verfolgung des Ätzprozesses mit Mikroauf- 
nahmen ergab, daß eine Ätzzeit von 1,5—2 Minuten am günstigsten ist. 
Nach Durchführung einer Anzahl Messungen wurden die Proben oft von neuem 
fein geschmirgelt und geätzt, so daß im Mittel die Dicke der Proben 0,15 bis 
0,2mm war. 

Abb. 2 und 3 zeigen eine Obfl.-Aufnahme in etwa achtzigfacher Ver- 
_ größerung und eine interferenzmikroskopische Aufnahme einer so behandelten 
Obfl. Zwischen je zwei Linien der interferenzmikroskopischen Aufnahme ist 


Abb. 2. cP 4-geätzte Germanium- Abb. 3. Interferenzmikroskopische 
Obfl. (Vergr. etwa 80fach) Aufnahme einer CP 4-geätzten Ger- 
manium-Obfl. (Vergr. etwa 80fach) 


Obfl. nur einen Flächeninhalt hat als die sie ideali- 
sierende Ebene, während das Verhältnis bei einer geschmirgelten Obfl. einen 
Faktor 5—10 ausmachen kann. Die geätzten Obfl. hatten im Durchschnitt 
eine Obfl.-Rekombinationsgeschwindigkeit von 100 cm/sec (gemessen nach 
der Methode von van Roosbroeck?). 

‘ Die fertig bearbeitete Probe wurde mit Ohmschen Kontakten versehen. 


u H. Flietner, Ann. Physik 3, 396 (1959). 
8) W. van Roosbroeck, Physic. Rev. 101, 1713 (1956). 
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2. Apparatur zur Messung des Feldeffektes 


Die so vorbereitete Probe wurde in den im Meßgefäß befindlichen Feld- 
kondensator eingebracht. Hier wurde sie gegenüber einer Metallelektrode 
angebracht, so daß Metallelektrode und Probe einen Kondensator bildeten. 
Als Dielektrikum wurde, je nach den Absichten, entweder Luft (bzw. die 
entsprechende Atmosphäre) oder SrTiO, (er = 300) verwendet. Die wirk- 
same Kondensatorfläche war !/, bzw. 1 cm?, je nachdem sich nur auf einer bzw. 
auf beiden Seiten der Probe Elektroden befanden. In Luft ergab sich bei Aus- 
nutzung von nur einer Probenseite eine Kapazität von einigen pF, bei Ver- 
wendung von SrTiO, etwa 50 pF. Der Faktor 300 konnte sich nicht voll aus- 
wirken, da die SrTiO,-Folien nicht so dünn gemacht werden konnten wie der 
Luftspalt (0,1 mm). 

Das Meßgefäß war abgedunkelt, um inneren und Obfl.-Photoeffekt zu 
vermeiden. Durch dieses Gefäß wurde das Gas geleitet, dessen Wirkung auf 
die Obfl. untersucht werden sollte. Das Ausgangsgas war reines N, und Q,. 
Für die Trocknung wurde das Gas über CaCl,, Silikagel und P,O, geleitet. 
Ein Haarhygrometer kontrollierte die Trocknung. Dann lief das Gas durch 
ein Glasgefäß, in dem eine Quecksilber-Niederdrucklampe zur Erzeugung 
von Ozon eingebaut war. Feuchtes Gas wurde erzeugt, indem das Ausgangs- 
gas durch Wasser hindurchgeleitet wurde. Um eine Kondensation im Meßgefäß 
weitgehend zu verhindern, war eine Flasche mit kleinen Glasperlen vorge- 
schaltet. 

Die Gasgeschwindigkeiten lagen je nach Absicht zwischen 2und 10 Litern 
pro Stunde. Eine Abhängigkeit des Effektes von der Gasgeschwindigkeit 
und eventuell damit verbundenen Temperaturschwankungen wurde nicht 
festgestellt. Die Messungen fanden bei Zimmertemperatur statt. 


Das trockene Gas hatte eine relative Feuchte unter 5%, das feuchte über 
95%. Damit ist man sicher, die hierdurch erzielbare größtmögliche Band- 
verbiegung erreicht zu haben?). 


Der elektrische Teil der Apparatur registrierte oszillographisch in einer 
Brückenschaltung die ‚Leitwertänderung der Probe in Abhängigkeit von der 
auf der Proben-Obfl. induzierten Ladung. Die Spannung am Feldkondensator 
war sinusförmig, ihr Effektivwert etwa 40—60 V, ihre u an 
200 und 400 Hz. eS aw 


it 
2 


Im folgenden soll über Messungen an vier verschiedenen Germanium- 
Proben berichtet werden. Die Proben hatten folgende Daten: 


Probe 11: p-Typ, 17 Ohm cm; Obfl. etwa 4° gegen 111-Kristallebene 


geneigt, 
Probe 12: p-Typ, 34 Ohm cm; Obfl. etwa 5,5° gegen 111-Kristallebene 
geneigt, 
Probe 13: p-Typ, 40 Ohm cm; Obfl. etwa 5° gegen 111-Kristallebene 
geneigt, 
Probe 15: n-Typ, 17 Ohm cm; Obfl. etwa 8° gegen 111-Kristallebene 
%) ALR. atest, Physic. Rev. 102, 381 
Ann. Physik. 7. Folge, Bd. 3 Ne : 27b 


2 
> 
as 
r 
Als 
ilung 
oauf- 
1 ist. 
euem 
(5 bis = N 
Ver- 
lelten 
ne ist 
tea 
4 
= 3 
7 


418 Annalen der Physik. 7. Folge. Band 3. 1959 


An Probe 11 wurde der in der Einleitung genannte Ozoneffekt festgestellt, 
_ der es gestattete, wie an Probe 12 und 15 gezeigt wird, das Spektrum der 
„schnellen“ Terme in reproduzierbarer Weise zu beeinflussen. Dieser Effekt 
trat an stabilen Obfl. auf, d.h. an Obfl., die nach dem Ätzen on. ov 
auch Wochen an Luft gelagert waren. 


Messungen an Probe 11 


An dieser Probe wurde in verschieden zusammengesetzten Zyklen der 
Einfluß von O, trocken und feucht, N, trocken und feucht und Ozon auf die 
mit CP 4 geätzte Obfl. untersucht. Die Probe wurde vor der Messung so 
_ lange gelagert, bis der Feldeffekt in einem O,-trocken-feucht-Zyklus stabil war. 
7 Es ergab sich, daß sich bei Zyklen, die nur O, trocken und feucht N, 

trocken und feucht enthielten, die Teilstücke der Obfl.-Leitfähigkeitekurven 
_ mit nur geringer Streuung zu einer gemeinsamen Kurve zusammensetzen ließen. 
Da durch diese Kurve die 


\ \ 160°2") 10, Coulfan? „schnellen“ Obfl.-Terme 
\ | 14 | fe 210°? 7 5 | festgelegt sind, finden bei 
diesen Zyklen keine oder 
\0s \és 12 | nur geringe chemische 
\ \ * | 2 Reaktionen statt, die die 
LA „schnellen“ Terme beein- 

\ 7 flussen. 


\ 


S 
4+ 


\ 3 WA 2 ust Anders wird dies, wenn 
der Zyklus Ozon enthilt. 

A 


Dann findet eine merkbare 


29 7 Anderung der Obfl.-Leit- 
\ YO, fähigkeitskurve statt. Am 
N NE YA deutlichsten wird diese 
= Veränderung bei einem 


543 74 Mo 7 2 3 4 5 60 reinen O,-trocken-Ozon- 

Vor Gu/em“ Zyklus (Abb. 4). Als Aus- 
Abb. 4. Ergebnis der Obfl.-Leitfähigkeitsmessungen an gangskurve der nicht mit 
Probe 11 und Bestimmung der Ladung in den ,,schnel- Ozon behandelten Obfl. er- 


] [77 T 
en‘‘ Termen gab sich die Kurve A. 


Obfl. „Leitfähigkeitakurve suschends flacher und weenie sich nach rechts. 
Im Endzustand nach etwa zwei Minuten erhält man die Kurve C. Nun wird 
der UV-Strahler wieder ausgeschaltet und die Obfl. der Probe ist von nun an 
nur trockenem O, ausgesetzt. Man bemerkt, daß selbst nach einigen Stunden 
die Obfl.-Leitfähigkeitskurve nicht auf die Ausgangskurve A zurückgeht, 
sondern in der Lage der Kurve B verbleibt. Nachdem die Probe dann wieder 
mehrere Stunden normaler Zimmerluft ausgesetzt wurde, konnte die Ausgangs- 
kurve A wieder sehr genau reproduziert werden. Es war lediglich festzustellen, 
daß der Punkt, an dem die Obfl. neutral ist, (?), etwas nach links verschoben 
war. Auf Grund dieser Tatsache lag die Vermutung nahe, daß der Einfluß 
der Feuchtigkeit entscheidend für den Rückgang auf die Ausgangskurve A 
. = Die Untersuchung dieses Sachverhaltes wurde an den Proben 12 und 15 
durchgeführt. 
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Messungen an Probe 15 


Zur Untersuchung des an Probe 11 beschriebenen Effektes wurden fol- 
fende Messungen durchgefiihrt: Zuerst wurde die Ausgangskurve A der nicht 
mit Ozon behandelten Obfl. bestimmt (Abb. 5). Dann wurde die Obfl. etwa 
2 Minuten ozonhaltiger At- 


mosphäre ausgesetzt. Hier- \ | Pa 

auf folgte O, trocken. 7 3:98 — 
Dieser Ozon-O,-trocken- 

Zyklus wurde einige Male C 

durchgeführt und zwar so \ \6s 


lange, bis an der O,-trocken- 
Kurve keine merkbaren N B \ 


Änderungen mehr auf- 
traten. Die sich dann er- ies \ 


\ 
gebende O,-trocken-Kurve Q \ 


B war flacher als die Aus- £) 
gangskurve, und der Neu- 
tralpunkt (?,) war nach 

links verschoben (?,). Das 


-6 + ewe ow ? 4109 
gleichzeitige Flacherwer- 


= der Kurve = ao Abb. 5. Ergebnis der Obfl.-Leitfähigkeitsmessungen 
dern des Neutralpunktes an Probe 15 im O,-trocken-Ozon-Zyklus 
konnte man an den Über- 


gangskurven, die zwischen A und B lagen, deutlich beobachten. Die Kurve B 
blieb in trockener O,-Atmosphäre stabil. Während 15—20 Stunden war keine 
Änderung festzustellen. Wurde die Obfl. dann einige Minuten mit feuchtem O, 
behandelt, dann ergab sich die Kurve C, die fast mit der Ausgangsk: rve 
übereinstimmte, jedoch war der Neutralpunkt nicht wieder in die alte Lage 
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Abb. 6. Gesamtkurve der Obfl.-Leitfähigkeit für Probe 15 und Bestimmung der Ladung 
in den „schnellen‘‘ Termen 


zurückgekehrt (1.). Nachdem die Obfl. dann etwa 15—20 Stunden der 
Zimmerluft ausgesetzt war (Meßgefäß geöffnet), ging auch der Neutral- 
punkt in die Ausgangslage zurück. 
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Von dieser Probe wurde dann noch die Obfl.-Leitfähigkeit in einem O,- 
_ trocken-naß- N,-trocken-naß- Zyklus ohne Ozonvorbehandlung bestimmt, um 
die Obfl. Leitfähigkeit für einen größeren Bereich zu erhalten. An diese 
Kurve wurde dann noch der Ozon-Ast angehängt. Diese Gesamtkurve ist 
in Abb. 6 gezeigt. 
Messungen an Probe 12 
An Probe 12 wurde festgestellt, daß die Verflachung der Obfl.-Leitfähig- 
keitskurve auch dann eintritt, wenn die O,-trocken-Ausgangskurve A sich 
schon sehr stark auf der Seite der n-Typ-Obfl. befindet (Abb. 7). Durch einen 
O,-trocken-Ozon-Zyklus, in diesem Falle dreimal hintereinander 2 Minuten 
Ozon, 10—15 Minuten O,-trocken, tritt wieder eine Verflachung ein (Kurve B). 
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Nach Anwendung einer feuchten Atmosphäre wird diese Wirkung wieder 
rückgängig gemacht. Eine merkbare Verschiebung des Neutralpunktes trat 
in diesem Fall nicht auf. Kurve C in Abb. 7 ist die Ozonkurve. 


Messungen an Probe 13 


An Probe 13 traten keine prinzipiell neuen Effekte auf. Es soll davon nur 
eine möglichst weit ausgesteuerte Obfl.-Leitfähigkeitskurve angegeben 
werden, um eine weitere Vergleichsmöglichkeit bis der Bestimmung aaa Obfl.- 


Terme zu haben (Abb. 8). | 


4. Bestimmung der Ladung in aon te a in | Abhängigkeit von 
der Bandverbiegung 


Die im voraufgegangenen Abschnitt angegebenen Obfl.-Leitfähigkeits- 
kurven werden nun in folgender von Brown?) beschriebenen Weise zur Be- 
stimmung der Ladung in den Obfl.-Termen verwendet. Die gemessene und 
die berechnete Obfl.-Leitfähigkeitskurven, 

9,(Q) und 0, (Q,) Q ist die Gesamt- 

ladung auf der Obfl., Q, der in der Raum- gu OF 
ladungsrandschicht befindliche Teil —, wer- ih 65 
den so in ein gemeinsames Diagramm ein- ur 

getragen, daß ihre Minima im Koordinaten- 
ursprung liegen (Abb. 9). 

Da die Obfl.-Leitfähigkeit durch die Ver- 
hältnisse in der Raumladungsrendschicht 
eindeutig bestimmt ist, gehören die Punkte 
gleicher Leitfähigkeit auf o, und g, zur glei- 
chen Bandverbiegung, falls nur noch das Vor- 
zeichen der Bandverbiegung wu, übereinstimmt, 
die Punkte also im vorliegenden Diagramm 
im gleichen Quadranten liegen. Die gesamte auf der Obfl. befindliche Ladung 
teilt sich auf in die Ladung der Raumladungsrandschicht Q, und den Anteil 
Q,, der sich in den Obfl.-Zuständen befindet. Da im vorliegenden Falle mit 
einer Frequenz von 200—400 Hz gearbeitet wird, ändert sich nur die Ladung 
in den „schnellen‘‘ Obfl.-Termen, die Ladung in den „langsamen“ Termen 
bleibt konstant. Trägt man das bei dieser Frequenz gemessene g, in das 
Diagramm ein, so ist das daraus entnommene Q, bis auf eine Konstante die 
Ladung in den ‚schnellen‘ Termen. Bei der Berechnung der Obfl.-Leit- 
fähigkeit o, wurde die Streuung der Ladungsträger berücksichtigt”). Die in 
dieser Weise ausgewerteten Messungen sind in den Abb. 4 bis 8 dargestellt. 

Im nächsten Abschnitt wird die Bestimmung der ‚schnellen‘ Obfl.-Terme 
aus der so erhaltenen Abhängigkeit Q, (u,) vorgenommen. 


5. Das Modell fiir die „schnellen“ Obfl.-Terme 


Der Verlauf der Q, (w,)-Kurven legt es nahe, die Beschreibung mit fol- 
gendem einfachen Modell zu versuchen: Die ‚schnellen‘‘ Obfl.-Terme sollen 
aus einem Akzeptor-Term E, mit der Konzentration N, pro em? Obfl. und 
einem Donator-Term E, mit der EEE N, bestehen (vel. Abb. 1). 
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Abb. 9. Bestimmung der Ladung 
in den Obfl.-Termen 
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Man sieht, daß durch Änderung der Bandverbiegung u, die relative Lage des 
 Fermi-Niveaus zu den Obfl.-Termen und damit deren Besetzung geändert 
wird. Die Beziehung zwischen [, und u, ist 


Die Ladung in den Obfl.-Termen Q, wird 

1 
= N, + — N, (ea + 1, (2) 

. wobei ¢, =E,/kT und e, = E,/k T ist. Hieraus entnimmt man, daß 
im Q,=Q,,=eN, lim Q, = Q,, = —eN,. (3) 
n—-00 n—>00 


» _ Aus der Kenntnis dieser beiden Grenzwerte kann man also die Konzentration 
_ der Terme entnehmen. Da Q, aber nur bis auf eine Konstante bekannt ist, 
ist nur die Differenz der beiden Grenzwerte zugänglich. Es ist 


Qua — Qua = (Ng + Na). (4) 
5 pe i Solange n< e,, wird Q, praktisch konstant bleiben. Uberstreicht 7 den Wert 


€,, dann werden die Akzeptoren umgeladen. Da sie im unbesetzten Zustand 
neutral sind, und bei Besetzung negativ geladen werden, wird Q, in der Um 
Qs gebung von e, um etwa e WN, ab- 
eNy 2 nehmen, falls Akzeptoren und 
Donatoren nicht zu dicht liegen und 
die Umladung der Donatoren in 
diesem Gebiet ohne Einfluß ist. Im 
Gebiet um e, wird dann der rest- 
liche Abfall um eN, stattfinden, 
und Q, für 7 >e, wieder praktisch 
konstant sein. Man erwartet also 
j einen Verlauf Q,(n), wie er schema- 
tisch in Abb. 10 gezeigt ist. Der ausgezogene Abschnitt ist der Ausschnitt aus 
dem Gesamtverlauf, der etwa durch das Experiment erfaßt wird. In dieser 
Unvollständigkeit liegt natürlich ein Unsicherheitsfaktor bei der Bestimmung 
der Termdaten. So läßt sich z. B. Gl. (4) nur dahingehend auswerten, daß 
man eine untere Grenze für die Konzentration der Terme erhält. 
sees Zur Bestimmung der Termdaten wird eine genauere Analyse des Kurven- 
_ verlaufs durchgeführt. Hierfür werden noch die beiden ersten Ableitungen 
von Q, (n) benötigt: 


0 


ENy4 
Abb. 10. Schematischer Gesamtverlauf von 
7 Q, (n ) 


= N, cosh— 25 (n — &4) + N, cosh-* — &) (9) 
4 dQ, 1 
e = N, cosh-* (n — tgh (n — 


+ N, cosh > — €,) tgh (n — 
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age des | d. h., die Kurve geht durch Spiegelung am Punkt (3 (e, + &,), 0) in sich über, 
»ändert | falls Akzeptoren- und Donatoren-Konzentrationen einander gleich sind. Au 
Umgekehrt kann man auf diese Gleichheit schließen, falls die Kurve eine x 
solche Symmetrie hat. Da die aus den Messungen ermittelten Kurven diese : 
(1) | Symmetrie aufweisen, braucht im folgenden nur der Fall N, = N, betrachtet u’. 
zu werden. a 


Fir die weitere Auswertung ist die Lage der Wendepunkte wichtig. Sie 
(2) | werden nach Gl. (6) wegen Te = 0 bestimmt aus 


(3) cosh? 4 (4 — €,) tgh 4 (n — &,) + cosh? (ny — e,) tgh 4 (n — &,) = 0. (8) 


tration ws 2 hat als Funktion von 7 für 2,—e,< 2,634 eine Nullstelle, für 


unt ist, — &, > 2,634 stets drei Nullstellen. Aus Abb. 10 und den Q, (w,)-Kurven 
di Abb. 4 bis 8 geht hervor, daß im vorliegenden Falle stets drei Wendepunkte 
(4) auftreten, also bestimmt ¢, — e, > 2,634 ist. Ist e, — e, sehr groß, so ent- 

nimmt man Gl. (8), daß die beiden äußeren Wendepunkte praktisch bei 
n Wert | 7 = &,, &, liegen. Aus den oben angeführten Symmetriegründen liegt der 
ustand | dritte Wendepunkt dann bei 7 = } (e, + &,). Diese Situation ist bei den aus 


sr Um- | den Messungen hervorgegangenen @ (n)-Kurven gegeben. öe 

V, ab- Die Steigung in den Wendepunkten wird nach Gl. (5) 
und 
en und 4 dQ 
ren in| e = == N, fürn = =H, (9) 
rrest- | - 2 fir te) (10) 
finden, 2 
ktisch Aus Gl. (9) kann man also in guter Näherung die Konzentration der dis- f 


t also | kreten Terme ermitteln. Es liegt allerdings hier wieder die Schwierigkeit vor, 

'hema- daß die Wendepunkte experimentell noch nicht erfaßt wurden. Man kann also 
itt aus | Lave und Konzentration der diskreten Terme nur abschätzen (vgl. die ent- 
dieser | sprechende Situation bei). 


nmun 
. re Nun kommt man aber zur Erklärung der gefundenen Q, (n)-Kurven nicht ? 

; mit zwei diskreten Termen aus. Nach Gl. (10) ergibt sich eine Steigung am 7 
urven- | mittleren Wendepunkt, die gegenüber dem Wert aus den Messungen viel zu ei 
Jungen klein ist. Es müssen also noch zusätzlich Terme vorhanden sein. Der Verlauf = 

um den mittleren Wendepunkt herum kann befriedigend wiedergegeben 
werden, wenn man eine stetige Verteilung mit konstanter Termdichte den zwei Be 


(5) | diskreten Termen überlagert. Wenn man die Steigung, die von den diskreten 
Termen herrührt, am mittleren Wendepunkt vernachlässigen kann, dann er- 
gibt sich die Termdichte der stetigen Verteilung N, aus dem Wert der Steigung . 

an diesem Wendepunkt zu 


6 
- N, hat die Dimension em”? eV-!, een 


10) W.L. Brown, H. C. Montgomery, Physic. Rev. 108, 865 (1956). oe 
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( 6. Bestimmung der Termdaten und Diskussion der Ergebnisse 


u I Mit Hilfe der im vorigen Abschnitt gewonnenen Formeln kann nun das 
Spektrum der ‚schnellen‘ Terme bestimmt werden. Das Ergebnis ist in 


Tab. 1 wiedergegeben. 


Auf Grund der Messungen ermittelte Termdaten 
N,cm-? eV-1 | 3 (eu + €4) | N,em? 
Probe 11: A 2,4 . 1011 0,8 >6 >8,5 - 1010 
B 3,4 - 104 
Probe 12: A 4,3 - 10" 0,8 8 ~1,2- 104 
B 6,0 - 10" 
Probe 13: 6,2. 10" 0,6 z1,1. 10" 
Probe 15: A 3,1. 100 0,0 >7 > 7,5 . 1010 
B 4,5. 104 


Betrachtet man zunächst die stetige Verteilung, so fällt auf, daß die Term- 
- konzentration von allen Proben in derselben Größenordnung liegt. Im Mittel 
_ ergibt sich für die nicht mit Ozon behandelte CP 4-geätzte Obfl. (Fall A) 
eine Dichte von N, = 4-10" cm-?eV-. Das Verhältnis der Konzentration 
nach Ozonbehandlung (Fall B) zu der Konzentration vor Ozonbehandlung 
_ (Fall A) ist bei allen Proben recht genau 1,4. Diese Tatsachen weisen darauf 
hin, daß bei allen Proben die gleiche Art der Obfl.-Terme vorliegt und, da die 
Proben in ihrem Leitfähigkeitscharakter recht unterschiedlich sind, die Vo- 
_ lumendaten keinen Einfluß auf die Obfl.-Terme haben. 

Auch bei den diskreten Termen fällt esauf, daß sie in ihren Daten im wesent- 
. 7 lichen übereinstimmen. Aus diesen Daten entnimmt man, daß bei Probe 12 
und 13 etwa die Wendepunkte gerade erfaßt wurden oder nicht weitab liegen, 
während bei Probe 11 und 15 die Wendepunkte noch etwas weiter außerhalb 
liegen. Sehr viel weiter auseinander als bei Probe 12 und 13 können die 
Wendepunkte aber nicht liegen, da sich sonst in dem gemessenen Bereich 
noch keine Krümmung in der Q, (n)-Kurve bemerkbar machen dürfte. Nimmt 
man nun extrapolierend an, daß die Wendepunkte für die Proben 11 und 15 
an den gleichen Stellen liegen, wie für die Proben 12 und 13, dann muß, um 
die gemessenen Kurven zu erklären, N, auf einen Wert erhöht werden, der 
re-ht gut mit den N ,-Werten für Probe 12 und 13 übereinstimmt. Die Mes- 
sungen können also an allen Proben mit dem gleichen diskreten Spektrum 
erklärt werden. Dieses Spektrum wird gegeben durch 


+e)=055 e,—-e,=85 N,=N,=12- 101 cm. 


Während beim kontinuierlichen Teil des Spektrums eine gewisse Streuung 
in den N,-Werten auftritt, ändert sich das diskrete Spektrum von Probe zu 
Probe nicnt merklich. Über die Änderung des diskreten Spektrums durch die 
Ozonbehandlung ließ sich keine nähere Aussage machen, da die gemessenen 
Kurven einen zu kleinen Abschnitt in der Umgebung des mittleren Wende- 
punktes lieferten. Nur soviel ließ sich sagen, daß durch die Ozonbehandlung 
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sicher keine zusätzlichen diskreten Terme erzeugt wurden, die näher am mitt- 
leren Wendepunkt liegen. 

Dann ist noch die Verschiebung des Neutralpunktes (f, nach ?, 
in Abb. 5) zu beachten. Diese Verschiebung gibt an, um welchen Betrag sich 
die Ladung in den ‚langsamen‘ Obfl.-Termen geändert hat. Eine solche 
Ladungsänderung kann zwei Ursachen haben; einmal eine Lageänderung des 
Fermi-Niveaus, zum andern eine Änderung der Termdichte. Die Annahme. 
daß es sich um eine bloße Lageänderung des Fermi-Niveaus handelt, führt 
aber zu einem Widerspruch. Hat sich nämlich auf Grund der Änderung der 
Termdichte der ‚schnellen‘ Terme im Fall der neutralen Obfl. das Fermi- 
Niveau relativ zu den ‚langsamen‘ Termen verschoben, so müßte es und damit 
die Ladung in den „langsamen“ Termen auf den Ausgangswert (1) zurück- 
gehen, sobald die Termdichte der ‚schnellen‘ Terme auf den Ausgangswert 
zurückgeht. Dies geschieht aber nicht, wie man an Hand der Lage von (f,) 
in Abb. 5 feststellt. Es liegt also hier eine Änderung der Termdichte der ,,lang- 
samen‘‘ Terme vor, die durch Ozon bewirkt wird. 

Die Ladungsänderung ist bei Probe 15 (Abb. 5) 2,6 - 10-° Coul/em?. Es 
sind also 1,6 - 10! Bindungen pro cm? aufgelöst worden, die vorher mit Elek- 
tronen besetzt waren. 

Bei den Messungen an Probe 12 brachte die Einwirkung von Ozon keine 
merkliche Verschiebung des Neutralpunktes mehr mit sich (Abb. 7). Die 


Ladung in den „langsamen“ Obfl.-Termen hat sich bei dieser Probe also nicht 
geändert; es sind keine Bindungen ,,langsamer‘‘ Terme aufgelöst worden. Die 
Ausgangskurve stellt aber hier bereits den Zustand einer ozonbehandelten 
Obfl. dar, deren Vergrößerung der Anzahl der ‚schnellen‘ Terme bereits 
wieder durch Feuchtigkeit rückgängig gemacht worden ist. Der Ausgangs- 


zustand bei dieser Probe ist also 


vergleichbar mit dem bei Probe 15 b fins 


durch die Kurve C dargestellten _ 
Zustand. Der Neutralpunkt t, der 
Ausgangskurve A ist etwa um 


der ,,langsamen‘‘ Terme mit Ak- 
zeptorcharakter gelöst. Diese , lang- 
samen‘ Terme kann man sich 
durch Adsorption von Sauerstoff 
entstanden denken, etwa in der 


in Abb. 11b gezeigten Weise. 
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Abb. 11. In Sauerstoffatmosphäre auftretende 

Bindungstypen an der Obfl.*), Ref. Hand- 

ler. a) Freie Valenz als Akzeptor; b) Ak- 

zeptortyp der Sauerstoffadsorption; c) GeO;- 
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Nimmt man an, daß diese Terme alle die gleiche energetische Lage haben, 
dann geht, wegen ihrer großen Anzahl, das Fermi-Niveau durch sie 


6 - 109 Coul/em? gegenüber dem 

Neutralpunkt bei nicht mit Ozon mat 

behandelter Obfl. verschoben. 

Die dargestellten Ergebnisse 

lassen sich unter Zuhilfenahme 1. 

folgender Vorstellungen beschrei- pair fie 

ben. Bei Probe 15 (Abb. 5, 6) wird Germanium co 

durch die Einwirkung von Ozon SEE 

eine gewisse Anzahl Bindungen O Sauerstoff 
C e Elektron 
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oben abgeschätzt, sich 1,6 - 1010 besetzte Bindungen pro cm? lösen, werden 
sich also insgesamt 3,2 - 1010 Bindungen pro cm? lésen. Ist es nun möglich, daß 
diese Terme sich in die zusätzlich auftretenden ,,schnellen‘‘ Terme verwandelt 
haben? Die positive Beantwortung dieser Frage hängt davon ab, ob der 
Bereich ÖE, über den sich die stetige Verteilung um den mittleren Wende- 
punkt herum erstreckt, groß genug ist. öE darf nicht wesentlich kleiner sein 
als die Energiedifferenz, die sich aus dem Abstand der beiden äußeren Wende- 
punkte ergibt, da sich sonst die Steigung der Q, (n)-Kurven in dem Bereich 
außerhalb der stetigen Verteilung als zu klein ergeben würde. 

_ Die Konzentrationsänderung der „schnellen“ Terme bei Probe 15 ist nach 
Tab. 1: 1,4 101! cm? eV-1, Soll die Änderung ihrer Gesamtzahl gleich der 
Anzahl der gelösten Bindungen werden, dann muß sich die stetige Verteilung 
über den Bereich 6# erstrecken, der sich aus 1,4 - 1011. öE = 3,2. 1010 zu 
öÖE = 0,2283eV = 9k T ergibt. Auch bei Probe 12 wird man für die Lage des 
Neutralpunktes die gleiche Reaktion verantwortlich machen. Mit den für 


ner a Da die Anzahl der „schnellen“ Some von allen Autoren für die gleiche 
Obfl.-Behandlung in der gleichen Größenordnung angegebe n wird und bei 
 Atmosphärenwechsel auch keine Änderung um eine Größenordnung oder 
‘mehr eintritt, ist für ihr Entstehen eine Ursache verantwortlich, die unabhängig 
von leicht desorbierbaren Adsorbenten ist*), Ref. Garrett u.a. Solche 
Ursachen sind Gitterbaufehler, Ecken, Kanten. Die Vermutung, daß diese 
Obfl.-Eigenschaften die ‚schnellen‘ Terme bestimmen, wird weiter begün- 
 stigt dadurch, daß beim Schmirgeln eine Obfl.-Vergrößerung entsteht, die 
‘bel weitem nicht ausreicht, die Vergrößerung der Obfl.-Rekombinations- 
geschwindigkeit zu erklären. Während die Vergrößerung der Obfl. durch 
_Schmirgeln höchstens etwa den Faktor 10 ausmachen kann (vgl. Abschn. 1), 
nimmt die Obfl.-Rekombinationsgeschwindigkeit um den Faktor 100 bis 1000 
zu. Eine dem Flächeninhalt der realen Obfl. proportionale Termdichie kann 
also diesen Effekt nicht erklären. Die Vermutung, daß Ecken und Kanten 
entscheidend zur Vermehrung der Terme beitragen, ist also naheliegend. 

Diese Terme zeigen Akzeptorcharakter, indem sie an einer freien Valenz 
eines Obfl.-Atoms ein Elektron aus dem Innern binden können, etwa in der 
in Abb. lla angegebenen Weise. Die Ozonbehandlung der Probe 15 bewirkt 
also bei einem Teil der Obfl.-Bindungen die Verwandlung des Bindungstyps 
der Abb. 11b in den der Abb. lla. An der Obfl. findet also die Reaktion 
O, + O 20, statt. 

Die Verminderung der ,,schnellen‘‘ Terme ohne wesentliche Ladungs- 
änderung in den ,,langsamen‘‘ Termen durch die Einwirkung der feuchten 
Atmosphäre läßt sich folgendermaßen erklären: Die aufgelöste Bindung muß 
einen neutralen Term ergeben. Eine solche Möglichkeit bietet die Adsorption 
von Sauerstoff in der in Abb. 1le angegebenen Weise als GeO,. Die freien 
Valenzen der Obfl.-Atome werden dabei so gebunden, daß pro adsorbiertes 
Sauerstoffatom ein Obfl.-Zustand beseitigt wird. 

Diese Terme sind sehr stabil, wie man daran erkennt, daß der Neutral- 
punkt ft, nicht ohne weiteres auf den Ausgangswert 1, zurückgeht. 


Mehrfache Behandlung in feuchter Atmosphäre (>95°/, relative Feuchte für 
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10—20 Minuten) änderte die Lage des Neutralpunktes nicht. In trockener 
Sauerstoffatmosphäre war während der Dauer von 20 Stunden auch keine 
Änderung zu bemerken. Ein Rückgang zum Ausgangs-Neutralpunkt ?, er- 
folgte, wenn die Probe etwa 15 Stunden der Normalatmosphäre (Zimmerluft) 
bei geöffnetem Meßgerät ausgesetzt war. 

Dieses Verhalten — Bildung des Zustandes innerhalb einiger Minuten 
und dann eine relative Stabilität — ist ein Merkmal der GeO,-Bindung, das 
von Dillon bemerkt und von Handler‘) zitiert wird und für das Vorliegen 
dieses Bindungstyps spricht. H,O wirkt bei dieser Reaktion als Katalysator. 

Die Entstehung „schneller“ Seems durch die Ozonbehandlung bei Probe 12 
ohne gleichzeitige Änderung der Ladung in den „langsamen“ Termen wird 
man nach dem Vorangegangenen nun so verstehen, daß die GeO,-Bindung 
durch Ozon gelöst und in die „schnelle“ Akzeptorbindung der Abb. 11a über- 
führt wird. 


Herrn Prof. Dr. ©. Hachenberg danke ich für die Förderung der Arbeit. 


Berlin-Adlershof, Heinrich-Hertz-Institut der Deutschen Akademie 
der Wissenschaften zu Berlin. 
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Berichtigung zur Arbeit: 
»Die Depolarisation der Natrium-Resonanzfluoreszenz“ 


Ann. Physik (7) 2, 393 (1959) 


Von W.Ermisch und R.Seiwert 


Bei der Ermittlung der Teilchendichte N der Natrium-Atome mit Hilfe 
der Formel N(T) = p(T)/kT sind, wie sich bei weiteren Berechnungen zum 
Thema der Arbeit herausgestellt hat, im Nenner 7-Werte aus der falschen 
Spalte des Rechenschemas eingesetzt worden. Um den dadurch entstandenen 
Fehler zu korrigieren, muß man sämtliche N-Werte in den Abb. 1 bis 3 mit 
dem Faktor T’°C/T°K multiplizieren. Dementsprechend ändern sich die 
Maßzahlen der unteren Skala der Abb. 1 und 2, und die Linearität des Maßstabes 
geht verloren, was hier jedoch ohne Bedeutung ist. In Abb. 3 ist die Teilchen- 
dichte als Parameter der Meßkurven angegeben. Entsprechend der erforder- 
lichen Korrektur muß es heißen 


statt 1- 1010 2,73. 109, statt 8.101 299.100 
statt 1. 1012 3,20. 1010, statt 1-10 3,79.10 
statt 2-101 6,81 - 1010, statt 2-10 7,90-10%, 
an statt 4-10" 1,43 . 104, statt 4 1012 1,65 - 1012, 

u; statt € - 1011 2,20. 104, statt 6 - 1012 2,53 - 1012, 


An der Deutung des Verlaufs der Meßkurven und den Folgerungen aus dem 
Meßresultat ändert sich durch die Berichtigung der N-Werte nichts. 


Berlin, II. Physikalisches Institut der Humboldt-Universität und In- 
stitut für Optik und Spektroskopie der Deutschen Akademie der Wissenschaften 


zu Berlin. 
Bei der Redaktion eingegangen am 26. Februar 1959. i im 
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DGL-PRESSLER LEIPZIG 


Lexikon der Kern- und Reaktortechnik 


Herausgegeben von K. H. HöckEr und K. WEIMER 
Zwei Bände Hütte-Format, 


zusammen 1685 Seiten mit 310 Abbildungen und Tabellen 
ee In Leinen gebunden DM 125.— ee 


Dieses Lexikon faBt zum ersten Male auf der Welt die Atomwissenschaft enzyklo- 
pädisch zusammen. Bei der Definition und Erläuterung der 3000 Stichworte sind 
die Ergebnisse der 2. Genfer Atomkonferenz „Atome für den Frieden‘ berück- 
sichtigt. 
Eine sehr ausführliche Isotopentabelle gibt der Forschung und Industrie die 
notwendigen Daten (Wirkungsquerschnitte, Halbwertzeiten usw.). Der umfang- 
reiche Tabellenteil ist zweckmäßig für die praktische Arbeit zusammengestellt und 
enthält Zahlen, die im allgemeinen nur nach mühsamer Kleinarbeit zugänglich sind. 


Die bereits arbeitenden und auch die geplanten Reaktoren aller Länder sind hier 

zum ersten Male tabellarisch zusammengefaßt und nach Aufbau, Arbeitsweise und 

Leistung aufgeschlüsselt. Rasches Auswerten der internationalen Fachliteratur 

gewährleisten die englischen und französischen Fachbezeichnungen nach jedem 
Stichwort und das englisch-französisch-deutsche Vokabular. 


In der Deutschen Demokratischen Republik beziehbar im Rahmen des Kontingents 


FRANCKH’SCHE VERLAGSHANDLUNG STUTTGART 


Quantentheorie des Atoms 
Von Prof. Dr. A. RUBINOWICZ, Warschau 


Übersetzung nach der 2. polnischen Auflage des Werkes „Kwantowa Teoria 
Atomu“. VII, 486 Seiten mit 38 Abb. im Text 1959. Gr. 8°. 
DM 29.80, Leinen DM 31.80 


Dieses Lehrbuch ist aus langjährigen Vorlesungen hervorgegangen Im ersten Teil des 
Buches wird die Entwicklung der älteren Quantentheorie dargestellt. Im zweiten Teil 
werden die Grundlagen und Anwendungen der neuen Quantentheorie gebracht. Dem 
Kapitel II ist ein mathematischer Anhang beigegeben, der die Sommerfeldsche Polynom- 
methode zur Lösung der Eigenwertprobleme der Quantentheorie in ihrer allgemeinen 
Fassung darstellt. Außerdem enthält sie einige spezielle im Buchtext nicht berücksich- 
tigte Tatsachen und Überlegungen. 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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Annalen der Physik. 7. Folge. Band 3. Heft 7-8. 1959 


Einfihrung in die Theorie der Elektronenoptik 
Von Dr. phil. JOHANNES PICHT 


o. Prof. für theoretische Physik und Optik, Pädagogische Hochschule Potsdam-Sanssouci 
(früher Techn. Hochschule Berlin-Charlottenburg) 
2., erweiterte Auflage 
VIII, 274 Seiten mit 70 Abbildungen im Text. 1957. Gr. 8° 
DM 29.20, Leinen DM 31.— 
Über die erste im Jahre 1939 erschienene Auflage des Buches schrieb Professor Dr. Brüche 
in der „Elektrotechnischen Zeitschrift‘; 
„Eine mit großer Sachkenntnis geschriebene Einführung in die Theorie, die durch ihr 
hohes Niveau und durch die klare Darstellung dem Theoretiker ein schnelles und 
tiefes Eindringen in die Grundlagen und den mathematischen Apparat der geome- 
trischen Elektronenoptik ermöglicht. Das Buch kann jedem, der sich ernstlich mit 
der Elektronenoptik befassen will, aufs wärmste empfohlen werden.“ 
Die neue Auflage ist auf den neuesten Stand der Wissenschaft gebracht worden, wodurch 
sich eine Erweiterung des Umfanges ergab. Der Verfasser war bemüht, den Charakter 
einer Einführung in das Gebiet zu wahren, insbesondere die mathematischen und ge- 
dauklichen Überlegungen soweit durchzuführen, wie dies zum ausreiehenden Verständnis 
der theoretischen Zusammenhänge erforderlich schien. 


Arbeitstagung Festkörperphysik I 


7s 


vom 8. bis 11. April 1954 in Dresden 
Redaktion Dr. K. W. Böer, Berlin 
219 Seiten mit 205 Abbildungen im Text. 1955. Gr. 8°. DM 21.— 


Neue Hütte: Der dargebotene Stoff hat keineswegs nur für den theoretischen Physiker 
oder Physiko-Chemiker, sondern auch für den physikalisch-chemisch orientierten Metall- 
kundler besonderes Interesse. Man kann sagen, daß die Niederlegung der Tagungsberichte 
in Buchform sehr wertvoll ist, da sich in den Vorträgen eine große Zahl fruchtbarer Ge- 
danken und Ansätze zur Theorienbildung zeigt. Die Beiträge sind von hohem wissen- 
schaftlichen Niveau. Thümniler 


Klassische Arbeiten deutscher Physiker 


Herausgegeben von der Physikalischen Gesellschaft in der 
Deutschen Demokratischen Republik 


Heft 1: W. C. RÖNTGEN, Grundlegende Abhandlungen über die X-Strahlen 
44 Seiten und 1 Porträt. 1954. Gr. 8°. DM 2.70 


Die Technik: Es ist sehr zu begrüßen, daß die Originalarbeiten Röntgens über seine 
Entdeckung der X-Strahlen neu herausgegeben worden sind. Das Studium der kleinen 
Schrift führt uns die Größe der Entdeckung Röntgens vor Augen und zeigt, daß er dank 
seiner Gründlichkeit bereits die wichtigsten Eigenschaften der X-Strahlen feststellen 
konnte. Das Studium der Abhandlungen kann sehr empfohlen werden. Freyer 


Heft 2: W. FRIEDRICH, P. KNIPPING, M.v. LAUE, Interferenzerseheinungen 
bei Röntgenstrahlen 


36 Seiten mit 3 Abb. im Text und 4 Tafeln. 1955. Gr. 8%, DM 3.— 


Röntgen- und Laboratoriumspraxis: Nachdem im ersten Heft die grundlegenden 
Arbeiten Röntgens abgedruckt worden waren, folgt jetzt im zweiten Heft der Serie die 
originalgetreue Wiedergabe der Veröffentlichung von Laue, Friedrich und Knipping aus 
dem Jahrgang 1913 der „Annalen der Physik‘ über Interferenzerscheinungen bei Röntgen- 
strahlen. Es ist außerordentlich reizvoll, den grundlegenden theoretischen Ausführungen 
zu folgen, die Laue über die vermutliche Wellennatur der Röntgenstrahlen und die Mög- 
lichkeit, sie mit Hilfe von Interferenzerscheinungen an dem regelmäßigen Raumgitter von 
Kristallen zu beweisen, in dieser Arbeit gemacht hat. Gj. 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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